macro
sim nrep

mcolumn resp acol bcol ccol blockcol ctrcol stdcol runcol CP SO RO noise
mconstant nrep sd

FFDesign 3 8;
CTPT ctrcol;
Blocks 1;
SOrder stdcol runcol;
Alias 3;
XMatrix blockcol acol bcol ccol.

let sd=sqrt(10/nrep)
let resp=120+.5*acol+2.0*acol*bcol+ccol
rand 8 noise;
norm 0 sd.
let resp=resp+noise

FFactorial;
Response resp;
Factors acol bcol ccol;
Terms acol bcol ccol acol*bcol acol*ccol bcol*ccol acol*bcol*ccol; 
InUnit 1;
Levels -1 1 -1 1 -1 1;
CTPT ctrcol;
GEffects.

endmacro