Chapter 16 SAS Analysis of Split-Plot Experiments

Read Data

proc import datafile=’C:\\Data\\FieldSplitPlotData.txt’
  dbms=TAB replace out=Field;
run;
proc print data=Field (obs=14);
run;

Fit Linear Mixed-Effects Model

proc mixed data=Field;
 class block geno fert;
 model y=geno fert geno*fert / ddfm=satterthwaite;
 random block block*geno;

Example Estimate Stateements

estimate ’geno 1’
  intercept 4 geno 4 0 0 fert 1 1 1 1
  geno*fert 1 1 1 1 0 0 0 0 0 0 0 0 / divisor=4 cl;

estimate ’geno 1 - geno 2’
  geno 4 -4 0
  geno*fert 1 1 1 1 -1 -1 -1 -1 0 0 0 0 / divisor=4 cl;

estimate ’geno 1 - geno 2 with no fertilizer’
  geno 1 -1 0 geno*fert 1 0 0 0 -1 0 0 0 0 0 0 0 / cl;
run;

Refit model with Fixed Block Effects

proc mixed data=Field;
  class block geno fert;
  model y=block geno fert geno*fert / ddfm=satterthwaite;
random block*geno;