OPTIONS LS=80;
DATA ONE;
  INPUT ID Depress Intcpt Linear LinChang Order OrdLin OrdLCha;
  WEEK = LINEAR;

CARDS;
    1    1    1 -2.5 -0.5     0    0    0
    1    1    1 -1.5  0.0     0    0    0
    1    1    1 -0.5  0.5     0    0    0
    1    1    1  0.5  0.5     0    0    0
    1    1    1  1.5  0.0     0    0    0
    1    1    1  2.5 -0.5     0    0    0

   .... MORE DATA .....

   75    5    1 -2.5 -0.5     1 -2.5 -0.5
   75    5    1 -1.5    0     1 -1.5    0
   75    5    1 -0.5  0.5     1 -0.5  0.5
   75    5    1  0.5  0.5     1  0.5  0.5
   75    5    1  1.5    0     1  1.5    0
   75    4    1  2.5 -0.5     1  2.5 -0.5
;

/* ************************** */
/* some random effects models */
/* ************************** */

PROC MIXED METHOD=ML COVTEST;
  CLASS ID;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT /SUB=ID G;

PROC MIXED METHOD=ML COVTEST;
  CLASS ID;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT LINEAR /SUB=ID TYPE=UN G GCORR;

PROC MIXED METHOD=ML COVTEST;
  CLASS ID;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT LINEAR LINCHANG /SUB=ID TYPE=UN G GCORR;


/* ******************************** */
/* some covariance structure models */
/* ******************************** */

PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  REPEATED WEEK / SUB=ID TYPE=CS R RCORR;
 
PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  REPEATED WEEK / SUB=ID TYPE=AR(1) R RCORR;
 
PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  REPEATED WEEK / SUB=ID TYPE=TOEP R RCORR;
 
PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  REPEATED WEEK / SUB=ID TYPE=UN R RCORR;



/* ***************************************************** */
/* some random effects models with autocorrelated errors */
/* ***************************************************** */

PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT /SUB=ID G;
  REPEATED WEEK /SUB=ID TYPE=TOEP(4) R RCORR;

PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT /SUB=ID TYPE=UN G GCORR;
  REPEATED WEEK / SUB=ID TYPE=AR(1) R RCORR;

PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT LINEAR /SUB=ID TYPE=UN G GCORR;
  REPEATED WEEK /SUB=ID TYPE=TOEP(3) R RCORR;

PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT LINEAR /SUB=ID TYPE=UN G GCORR;
  REPEATED WEEK / SUB=ID TYPE=AR(1) R RCORR;

PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT LINEAR LINCHANG /SUB=ID TYPE=UN G GCORR;
  REPEATED WEEK / SUB=ID TYPE=TOEP(3) R RCORR;

PROC MIXED METHOD=ML COVTEST;
  CLASS ID WEEK;
  MODEL Depress = Linear LinChang Order OrdLin OrdLCha/SOLUTION;
  RANDOM INTERCEPT LINEAR LINCHANG /SUB=ID TYPE=UN G GCORR;
  REPEATED WEEK / SUB=ID TYPE=AR(1) R RCORR;