/* Split Strip Plot Design */
options nocenter ps=64 ls=80;
data splstr;
   infile 'splstr.dat';
   input A S B C Y;
/* analyse the whole plot: trt = A */
proc sort; by A S;
proc means noprint; by A S;
   var Y;
   output out=WP mean=YWP;
proc glm;
   class A;
   model YWP = A / ss1;
/* analyse row split plot: trt = A,B */
data b; set splstr;
proc sort; by A S B;
proc means noprint; by A S B;
   var Y;
   output out=SPB mean=YSPB;
proc glm;
   class A S B;
   model YSPB = A S B A*B / ss1;
   test h = A e = S;
/* analyse column split plot: trt = A,C */
data c; set splstr;
proc sort; by A S C;
proc means noprint; by A S C;
   var Y;
   output out=SPC mean=YSPC;
proc glm;
   class A S C;
   model YSPC = A S C A*C / ss1;
   test h = A e = S;
/* analyse strip plot: trt = B,C */
data d; set splstr;
proc glm;
   class B C S;
   model Y = S B B*S C C*S B*C / ss1;
   test h = B e = B*S;
   test h = C e = C*S;
/* analyse split-strip plot: trt = A,B,C */
proc glm;
   class A B C S;
   model Y = A S B A*B B*S C A*C C*S B*C A*B*C / ss1;
   test h = A e = S;
   test h = B A*B e = B*S;
   test h = C A*C e = C*S;