***************************************************************;
* Generate a multistage experiment with 4 whole plot factors
* and 3 subplot factors in 16 runs.
* I=AwBwCw=DwAsBs=AwAsCs.
*
* Generate simulated random normal y-values, just to demonstrate
* how to find the aliasing structure using Proc GLM.
****************************************************************;
data multistage;
   do Aw=-1,1; do Bw=-1,1; do Cw=-1,1; do Dw=-1,1;
   do As=-1,1; do Bs=-1,1; do Cs=-1,1;
   y=rannor(10981);
   if Aw*Bw*Cw=1 and Dw*As*Bs=1 and Aw*As*Cs=1 then output;
   end; end; end; end;
   end; end; end;
proc print data=multistage;
proc glm data=multistage;
   model y = Aw | Bw | Cw | Dw | As | Bs| Cs @ 2 / aliasing;
run;