The syntax would be the same as for the above design, so I will repeat that code, and then insert your factors.
/* previous poster */
PROC GLIMMIX data=data;
class site origin block cut tree;
model recovery = site|origin|cut/ddfm=kr;
random block block*site block*origin block*cut block*site*origin block*site*cut block*origin*cut block*origin*site*cut;;
/* change variable names to organic, system, nitrogen */
PROC GLIMMIX data=data;
class block organic system nitrogen;
model recovery = organic|system|nitrogen/ddfm=kr2;
random block block*organic block*system block*nitrogen block*organic*system block*organic*nitrogen block*organic*system*nitrogen;
Again, I want to point out that you probably do not have enough data to estimate the three-way and four-way effects, and they probably should be considered residual error in any case. To improve the ability to converge, you may then want this for the RANDOM statement:
random intercept organic system nitrogen/subject=block;
SteveDenham
... View more