I am working on sas code to test for treatment effects on soil pH through time. I have 4 plots (reps), two of which are fertilized, and two that are not fertilized (Treat) 1/2 of each plot was plowed (plow) I measured soil pH at four uneven time intervals (DOY), so end up with 32 data points. Here is the code I developed - am I missing something? Thanks! proc mixed data=soil; class treat plow reps doy; model pH = treat|plow|doy / Lcomponents; repeated doy / sub=treat type=cs; run;
... View more