/* Chapter 5 ++++++++++++++++++++++++++++++++++++++++++++ pages 70-72 Programs 5.1-5.3 */ data d3; infile 'c:\book\RepeatedMeasure\ratlong.dat' missover; input no y id day; run; proc glm data=d3; class group; model y = day /solution clparm; proc mixed data = d3 method=reml covtest; class id ; model y = day / s cl ddfm=sat; random intercept / subject=id g gcorr ; repeated / type=simple subject=id r rcorr; run ; /* type = simple */ proc mixed data = d3 method=reml covtest; class id ; model y = day / s ddfm=sat cl ; random intercept day / type=simple subject=id g gcorr ; repeated / type=simple subject=id r rcorr; run ; /* type = un */ proc mixed data = d3 method=reml covtest; class id ; model y = day / s ddfm=sat cl ; random intercept day / type=un subject=id g gcorr ; repeated / type=simple subject=id r rcorr; run ; /* ++++++++++++++++++++++++++++++++++++++++++++ page 76 Program 5.4 */ data d1; infile 'c:\book\RepeatedMeasure\experimentRat.dat' missover; input id group week y; group=group-1; run; proc mixed data=d1 method=reml covtest; class id week / ref=first ; model y = group week group*week / s cl ddfm=sat ; random intercept / subject= id g gcorr ; repeated week/ type = simple subject = id r rcorr ; estimate 'mean CFB ' group*week 1 1 1 -3 / divisor=3 cl alpha=0.05; run;