proc sort; by obs; *model in Table 13.1 and 13.2; proc mixed cl covtest ; class dyad obs sobvdru; model yt = trial sober drunk sober*ytm1o drunk*ytm1o sober*ytm1p drunk*ytm1p /s ddfm=SATTERTH noint; random sober drunk / sub=dyad type=un gcorr; random sober*ytm1o drunk*ytm1o / sub=dyad type=un gcorr; repeated sobvdru /type=csh subject=obs r; run;quit; *final growth curve model on page 370; proc mixed cl covtest ; class dyad obs m1m2; model yt = trial sovdru trial*sovdru /s ddfm=SATTERTH ; random intercept trial /type=un sub=dyad gcorr ; repeated m1m2 / sub=obs type=csh r; run;quit; * saturated growth model; proc mixed cl covtest ; class dyad obs m1m2; model yt = sober drunk trial*sober trial*drunk /s ddfm=SATTERTH noint; random sober drunk trial*sober trial*drunk / sub=dyad type=un; repeated m1m2 / sub=obs type=csh r; run;quit;