JOB '4 x 5 Rectangular lattice '
"4x5 rectangular lattice - refer Cochran and Cox; Page 435"
UNITS[60]
Factor[levels=5;nvalues=60] PF[1...3]
FACTOR[LEVELS=3] Rep
FACTOR[LEVELS=4] Plot
& [LEVELS=5] Block
& [LEVELS=5] PF[1...3]
& [LEVELS=20;values=10,12,11,9,2,4,1,3,20,18,17,19,5,6,8,7,15,14,13,16,\
3,11,20,7,4,16,12,8,10,1,18,14,17,3,13,5,6,15,2,19,8,11,15,18,19,1,5,12,\
6,3,13,10,2,9,16,20,7,4,14,17] Treat
"READ"
GENERATE Rep,Block,Plot
& [TREATMENTS=Treat;REPLICATES=Rep;BLOCKS=Block] PF[1...3]
BLOCKS Rep/Block/Plot
TREATMENTS Treat//PF[]
ANOVA
"
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]
"
ENDJOB