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;
Looks OK to me
I agree with Peter, the code formulation is perfect. I would hesitate on specifying a compound symmetry error structure, as it assumes equal correlations between days, no matter how far separated in time, and equal covariances at each time point. However, because of the small N for timepoints, and unequal spacing in time, this may be the only candidate that will be able to be fit.
Steve Denham
Thanks!
Don't miss out on SAS Innovate - Register now for the FREE Livestream!
Can't make it to Vegas? No problem! Watch our general sessions LIVE or on-demand starting April 17th. Hear from SAS execs, best-selling author Adam Grant, Hot Ones host Sean Evans, top tech journalist Kara Swisher, AI expert Cassie Kozyrkov, and the mind-blowing dance crew iLuminate! Plus, get access to over 20 breakout sessions.
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.