12-12-2016 11:53 AM
how to analyze data from a split-split-plot design with multi years? Would the following code work?
wc is whole plot factor, n is the sub plot factor, and sr is the sub-sub plot factore.
proc mixed data=file2;
class year rep wc n sr;
model plants=rep year*rep year|wc|n|sr;
random year*rep*wc year*rep*wc*n;
12-12-2016 12:50 PM - edited 12-12-2016 12:55 PM
Same plots from year to year, so that these are repeated measures in time as well?
At least remove rep from the model statement, and move to a RANDOM statement. Maybe something like:
proc mixed data=file2; class year rep wc n sr; model plants= year|wc|n|sr; random intercept year year*wc year*wc*n/subject=rep; repeated year/type=un subject=rep; run;;