It is often of interest to test whether the mean values of the dependent variable increases or decreases monotonically with certain factors. Hirotsu and Srivastava (2000) demonstrate one approach by using data (Moriguchi, 1976). The data consist of ferrite cores subjected to four increasing temperatures. The response variable is the magnetic force of each core.
data FerriteCores; do Temp = 1 to 4; do rep = 1 to 5; drop rep; input MagneticForce @@; output; end; end; datalines; 10.8 9.9 10.7 10.4 9.7 10.7 10.6 11.0 10.8 10.9 11.9 11.2 11.0 11.1 11.3 11.4 10.7 10.9 11.3 11.7 ;
The method presented by Hirotsu and Srivastava (2000) to test whether the magnetic force of the cores rises monotonically with temperature depends on the lower confidence limits of the isotonic contrasts of the force means at each temperature, adjusted for multiplicity. The corresponding isotonic contrast compares the average of a particular group and the preceding groups with the average of the succeeding groups. You can compute adjusted confidence intervals for isotonic contrasts by using the LSMESTIMATE statement.
The following statements analyze the FerriteCores
data as a one-way design and multiplicity-adjusted lower confidence limits for the isotonic contrasts. For the multiplicity adjustment, the LSMESTIMATE
statement employs simulation, which provides adjusted p-values and lower confidence limits that are exact up to Monte Carlo error.
proc mixed data=FerriteCores; class Temp; model MagneticForce = Temp; lsmestimate Temp 'avg(1:1)<avg(2:4)' -3 1 1 1 divisor=3, 'avg(1:2)<avg(3:4)' -1 -1 1 1 divisor=2, 'avg(1:3)<avg(4:4)' -1 -1 -1 3 divisor=3 / adjust=simulate(seed=1) cl upper; ods select LSMestimates; run;
The results are shown in Output 65.10.1.
Output 65.10.1: Analysis of LS-Means with Isotonic Contrasts
Least Squares Means Estimates Adjustment for Multiplicity: Simulated |
|||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Effect | Label | Estimate | Standard Error | DF | t Value | Tails | Pr > t | Adj P | Alpha | Lower | Upper | Adj Lower | Adj Upper |
Temp | avg(1:1)<avg(2:4) | 0.8000 | 0.1906 | 16 | 4.20 | Upper | 0.0003 | 0.0010 | 0.05 | 0.4672 | Infty | 0.3771 | Infty |
Temp | avg(1:2)<avg(3:4) | 0.7000 | 0.1651 | 16 | 4.24 | Upper | 0.0003 | 0.0009 | 0.05 | 0.4118 | Infty | 0.3337 | Infty |
Temp | avg(1:3)<avg(4:4) | 0.4000 | 0.1906 | 16 | 2.10 | Upper | 0.0260 | 0.0625 | 0.05 | 0.06721 | Infty | -0.02291 | Infty |
With an adjusted p-value of 0.001, the magnetic force at the first temperature is significantly less than the average of the other temperatures. Likewise, the average of the first two temperatures is significantly less than the average of the last two (p = 0.0009). However, the magnetic force at the last temperature is not significantly greater than the average magnetic force of the others (p = 0.0625). These results indicate a significant monotone increase over the first three temperatures, but not across all four temperatures.