" Example C-C-418: 3x4 rectangular lattice"
UNITS [36]
FACTOR [LEVELS=3] Rep,Plot
& [LEVELS=4] Block
& [LEVELS=4] PF[1...3]
& [LEVELS=12] Treat
READ Treat,Y
10 7 12 12 11 11 2 9 3 4 1 16 7 16 9 15 8 23 4 0 5 3 6 11
3 11 6 20 9 17 1 17 11 8 8 19 12 9 2 10 5 6 10 6 4 5 7 14
8 20 6 15 12 10 9 16 10 9 2 15 11 6 3 3 4 1 5 11 1 22 7 17 :
GENERATE Rep,Block,Plot
& [TREATMENTS=Treat;REPLICATES=Rep;BLOCKS=Block] PF[1...3]
BLOCKS Rep/Block/Plot
TREATMENTS Treat//PF[]
ANOVA Y
"PF[1], PF[2] and PF[3] are mutually non-orthogonal, thus no means
are given for treatments. Correct means can be obtained from the
partial effects (accessed via AKEEP). N.B. The means are based on
estimates from the Rep.Block.Plot stratum only."
AKEEP TERMS=PF[]+Treat; PARTIALEFFECTS=T[1...4]
CALCULATE X[1...4]=NEWLEVELS(PF[],Treat;!(#T[1]),!(#T[2]),!(#T[3]),!(#T[4]))
& Y=MEAN(Y)+VTOTAL(X)
TABULATE Y; MEANS=T[4]
PRINT T[4]