"REML Analysys of Cedara1.gsh - Rectangular Lattice"
set[outprint=dots]
VCOMPONENTS [FIXED=Cultivar] RANDOM=Rep/Block; INITIAL=1,1; CONSTRAINTS=positive,\
positive

For I=totDW[1...7],TOT_DW
REML [PRINT=model,components,means,waldTests; PSE=differences; MVINCLUDE=*; METHOD=AI] I
Endfor