/* 
Chapter 9 
++++++++++++++++++++++++++++++++++++++++++++
page 209  Program 9.1 */ 


data epilong ;
  infile 'c:\book\RepeatedMeasure\epilepsy_long.txt' missover;
  input no  subject treatment  visit seizure age;
  if visit=0 then post=0; else post=1;
  if visit=0 then logt=log(4); else logt=log(1);
  run;

  ods graphics on;

 /* Model I */
 proc glimmix data=epilong method=quad (qpoints=7)
   plots=studentpanel( conditional blup);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  visit treatment*visit 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept /  subject= subject g gcorr ;
   estimate 'progabide at 8w' visit 0 0 0 1 -1
                       treatment*visit 0 0 0 1 -1 /divisor=1  cl  alpha=0.05;
run;

/* Model IIa */
 proc glimmix data=epilong method=quad (qpoints=7)
   plots=studentpanel( conditional blup);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  visit treatment*visit 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept post/  type=simple subject= subject g gcorr ;
   estimate 'progabide at 8w' visit 0 0 0 1 -1
                       treatment*visit 0 0 0 1 -1 /divisor=1  cl  alpha=0.05;
run;

/* Model IIb */
 proc glimmix data=epilong method=quad (qpoints=7)
   plots=studentpanel( conditional blup);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  visit treatment*visit 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept post / type=un subject= subject g gcorr ;
   estimate 'progabide at 8w' visit 0 0 0 1 -1
                       treatment*visit 0 0 0 1 -1 /divisor=1  cl  alpha=0.05;
run;


/* 
++++++++++++++++++++++++++++++++++++++++++++
page 216  Program 9.3 */ 

/*  Model IV */
proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  post treatment*post 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept /  subject= subject g gcorr ;
 run ;


/*  Model Va */
proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  post treatment*post 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept post / type=simple subject= subject g gcorr ;
 run ;


/*  Model Vb */
proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  post treatment*post 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept post / type=un subject= subject g gcorr ;
 run ;


/* 
++++++++++++++++++++++++++++++++++++++++++++
page 220  Program 9.4 */ 


data epilong ;
  infile 'c:\book\RepeatedMeasure\epilepsy_long.txt' missover;
  input no  subject treatment  visit seizure age;
  if visit=0 then post=0; else post=1;
  if visit=0 then logt=log(4); else logt=log(1);
  xvisit=visit;
  run;

  /* Model VI */
  proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  xvisit treatment*xvisit 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept / subject= subject g gcorr ;
run;


  /* Model VIIa */
  proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  xvisit treatment*xvisit 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept xvisit / type=simple subject= subject g gcorr ;
run;


  /* Model VIIb */
  proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  xvisit treatment*xvisit 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept xvisit / type=un subject= subject g gcorr ;
run;



/* 
++++++++++++++++++++++++++++++++++++++++++++
page 224  program 9.5  */ 


data epilong ;
  infile 'c:\book\RepeatedMeasure\epilepsy_long.txt' missover;
  input no  subject treatment  visit seizure age;
  if visit=0 then post=0; else post=1;
  if visit=0 then logt=log(4); else logt=log(1);
  xvisit=visit;
  xvisit2=xvisit*xvisit;
  run;

  /* Model VIQ */
  proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  xvisit xvisit2 treatment*xvisit treatment*xvisit2 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept / subject= subject g gcorr ;
run;

 /* Model VIIQ */
  proc glimmix data=epilong method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  age treatment  xvisit xvisit2 treatment*xvisit treatment*xvisit2 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept xvisit xvisit2/ type=un subject= subject g gcorr ;
run;



/* 
++++++++++++++++++++++++++++++++++++++++++++
page 230  Program 9.6  */ 


data epilongB ;
  infile 'c:\book\RepeatedMeasure\epilepsyB_long.txt' missover;
  input no  subject treatment  visit seizure age base;
  logt=log(1);
  logbase=log(base/4);
  run;

 /* Model VIII */
proc glimmix data=epilongB method=quad (qpoints=7);
   class subject visit/ref=first  ;
   model seizure  =  logbase  age treatment  visit treatment*visit 
                                                        / d=poisson link=log offset=logt s cl ;
   random intercept /  subject= subject g gcorr ;
   estimate 'treatment at 2w'
      treatment 1  treatment*visit  0 0 0 1 /divisor=1 cl alpha=0.05;
   estimate 'treatment at 4w'
      treatment 1  treatment*visit  1 0 0 0 /divisor=1 cl alpha=0.05;
   estimate 'atment at 6w'
      treatment 1  treatment*visit  0 1 0 0 /divisor=1 cl alpha=0.05;
   estimate 'treatment at 8w'
      treatment 1  treatment*visit  0 0 1 0 /divisor=1 cl alpha=0.05;
  run ;

/* Model  IX  */

proc glimmix data=epilongB method=quad (qpoints=7);
   class subject visit ;
   model seizure  =  logbase age  treatment
                                / d=poisson link=log offset=logt s cl ;
   random intercept /  subject= subject g gcorr ;
  run ;