options linesize=80 pagesize=60;

data panss;
  infile 'panss-data-stacked.data';
  input id y0 score time week group;
run;

data panss; set panss;
  time8 = time - 8;
run;

**************************************;
***** Exchangeable Cov           *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	repeated / type=cs  subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Ran Int + Slope            *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	random intercept time / type=un subject=id g;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Ran Int + AR(1)            *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	random intercept / subject=id g;
	repeated / type=sp(pow)(time)  subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Ran Int + AR(1) + e        *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	random intercept / subject=id g;
	repeated / type=sp(pow)(time)  subject=id local;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Ran Int + Slope + AR(1)    *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	random intercept time / type=un subject=id g;
	repeated / type=sp(pow)(time)  subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Ran Int + Slope + AR(1) + e  ***;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	random intercept time / type=un subject=id g;
	repeated / type=sp(pow)(time)  subject=id local;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Ran Int + Slope + AR(1) + e  ***;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	random intercept time / type=un subject=id g;
	repeated / type=sp(exp)(time)  subject=id local;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Unstructured               *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time group time*group / s;
	repeated / type=un subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 
                          time*group -8 8 0  / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 
                          time*group -8 0 8 / cl alpha=0.05;
run;

**************************************;
***** Ran Int + Slope            *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time8 group time8*group / s;
	random intercept time8 / type=un subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 / cl alpha=0.05;
run;

**************************************;
***** Unstructured               *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time8 group time8*group / s;
	repeated / type=un subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 / cl alpha=0.05;
run;

**************************************;
***** Unstructured (y0*time8)    *****;
**************************************;
proc mixed data=panss method=ml;
        class id group;
	model score = y0 time8 y0*time8 group time8*group / s;
	repeated / type=un subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 / cl alpha=0.05;
run;

**************************************;
***** Unstructured (time8 factor) ****;
**************************************;
proc mixed data=panss method=ml;
        class id group time8;
	model score = y0 time8 group time8*group / s;
	repeated / type=un subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 / cl alpha=0.05;
run;

**************************************;
***** Unstructured (y0*time8 factor) *;
**************************************;
proc mixed data=panss method=ml;
        class id group time8;
	model score = y0 time8 y0*time8 group time8*group / s;
	repeated / type=un subject=id;
        estimate 'Group 2 versus Group 1 at 8 weeks'
                          group -1 1 0 / cl alpha=0.05;
        estimate 'Group 3 versus Group 1 at 8 weeks'
                          group -1 0 1 / cl alpha=0.05;
run;