/*
Chapter 7
++++++++++++++++++++++++++++++++++++++++++++
page 104  Program 7.1  */


data dbb;
  infile 'c:\book\RepeatedMeasure\BBlong.txt' missover;
  input no subject drug length treatment visit bdi;
  if visit=0 then post=0; else post=1;
  if visit=2 then w1=1; else w1=0;
  if visit=3 then w2=1; else w2=0;
  if visit=5 then w3=1; else w3=0;
  if visit=8 then w4=1; else w4=0;
  run;

ods graphics on;

/* Model I */
proc mixed data=dbb method=reml covtest
  plots=studentpanel( conditional blup) ;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment visit treatment*visit  / s cl ddfm=sat ;
  random intercept / subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;

/* Model IIa */
proc mixed data=dbb method=reml covtest
  plots=studentpanel( conditional blup) ;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment visit treatment*visit  / s cl ddfm=sat ;
  random intercept post / type=simple subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;

/* Model IIb */
proc mixed data=dbb method=reml covtest
  plots=studentpanel( conditional blup) ;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment visit treatment*visit  / s cl ddfm=sat ;
  random intercept post / type=un subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;

/* Model III */
proc mixed data=dbb method=reml covtest
  plots=studentpanel( conditional blup) ;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment visit treatment*visit  / s cl ddfm=sat ;
  random intercept w1 w2 w3 w4 / type=un subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;


/*
++++++++++++++++++++++++++++++++++++++++++++
page 111  Program 7.2  */

/* Model IV */
proc mixed data=dbb method=reml covtes ;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment post treatment*post  / s cl ddfm=sat ;
  random intercept / subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;


/* Model Va */
proc mixed data=dbb method=reml covtes ;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment post treatment*post  / s cl ddfm=sat ;
  random intercept post / type=simple subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;

 
/* Model Vb */
proc mixed data=dbb method=reml covtes ;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment post treatment*post  / s cl ddfm=sat ;
  random intercept post / type=simple subject= subject g gcorr ;
  repeated visit / type = un subject = subject r  rcorr ;
 run ;



/*
++++++++++++++++++++++++++++++++++++++++++++
page 114  Program 7.3  */


proc mixed data=dbb method=reml covtes ;
  class subject  visit treatment /ref=first  ;
  model bdi = drug length treatment post treatment*post  / s cl ddfm=sat ;
  random intercept / type=simple subject= subject g gcorr ;
  random post / type=simple subject= subject g gcorr group=treatment ;
  repeated visit / type = simple subject = subject r  rcorr ;
run ;



/*
++++++++++++++++++++++++++++++++++++++++++++
page 117  Program 7.4  */


data dbb;
  infile 'c:\book\RepeatedMeasure\BBlong.txt' missover;
  input no subject drug length treatment visit bdi;
 if visit=0 then post=0; else post=1;
 xvisit=visit;

/* Model VI */
proc mixed data=dbb method=reml covtest;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment xvisit treatment*xvisit  / s cl ddfm=sat ;
  random intercept /  subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;


 /* Model VIIa */
proc mixed data=dbb method=reml covtest;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment xvisit treatment*xvisit  / s cl ddfm=sat ;
  random intercept xvisit/  subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;


 /* Model VIIb */
proc mixed data=dbb method=reml covtest;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment xvisit treatment*xvisit  / s cl ddfm=sat ;
  random intercept xvisit/ type=un  subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;

 
/*
++++++++++++++++++++++++++++++++++++++++++++
page 120  Program 7.5  */


data dbb;
  infile 'c:\book\RepeatedMeasure\BBlong.txt' missover;
  input no subject drug length treatment visit bdi;
 if visit=0 then post=0; else post=1;
 xvisit=visit;
 xvisit2=visit*visit/8;


/* Model VIIQ */
proc mixed data=dbb method=reml covtest;
  class subject  visit/ref=first  ;
  model bdi = drug length treatment xvisit xvisit2 treatment*xvisit treatment*xvisit2 / s cl ddfm=sat ;
  random intercept xvisit xvisit2 / type=un subject= subject g gcorr ;
  repeated visit / type = simple subject = subject r  rcorr ;
 run ;


/*
++++++++++++++++++++++++++++++++++++++++++++
page 126  Program 7.7  */


data dbb2;
  infile 'c:\book\RepeatedMeasure\BBlongAncova.txt' missover;
  input no subject drug length treatment visit bdi bdi0;
  if visit=3 then w2=1; else w2=0;
  if visit=5 then w3=1; else w3=0;
  if visit=8 then w4=1; else w4=0;
run;


/* Model VIIIa */
   proc mixed data=dbb2 method=reml covtest;
   class subject visit /ref=first  ;
   model bdi = drug length bdi0 visit treatment treatment*visit / s cl ddfm=sat ;
   random intercept /  subject= subject g gcorr ;
    repeated visit/ type =simple subject = subject r  rcorr ;
    estimate 'CFB:   2m   ' 
      treatment  1  treatment*visit  0 0 0 1/ divisor=1 cl alpha=0.05;
    estimate 'CFB:   3m   ' 
      treatment  1  treatment*visit  1 0 0 0 / divisor=1 cl alpha=0.05;
    estimate 'CFB:   5m   ' 
      treatment  1  treatment*visit  0 1 0 0 / divisor=1 cl alpha=0.05;
    estimate 'CFB:   8m   ' 
      treatment  1  treatment*visit  0 0 1 0 / divisor=1 cl alpha=0.05;
   run;


/* Model VIIIb */
   proc mixed data=dbb2 method=reml covtest;
   class subject visit /ref=first  ;
   model bdi = drug length bdi0 visit treatment treatment*visit / s cl ddfm=sat ;
   random intercept w2 w3 w4/  subject= subject g gcorr ;
    repeated visit/ type =simple subject = subject r  rcorr ;
    estimate 'CFB:   2m   ' 
      treatment  1  treatment*visit  0 0 0 1/ divisor=1 cl alpha=0.05;
    estimate 'CFB:   3m   ' 
      treatment  1  treatment*visit  1 0 0 0 / divisor=1 cl alpha=0.05;
    estimate 'CFB:   5m   ' 
      treatment  1  treatment*visit  0 1 0 0 / divisor=1 cl alpha=0.05;
    estimate 'CFB:   8m   ' 
      treatment  1  treatment*visit  0 0 1 0 / divisor=1 cl alpha=0.05;
   run;[

/* Model IX */
 proc mixed data=dbb2 method=reml covtest;
   class subject visit ;
   model bdi = drug length bdi0 treatment  / s cl ddfm=sat ;
   random intercept /  subject= subject g gcorr ;
   repeated visit / type = simple subject = subject r  rcorr ;
 run;