libname M236 'd:\m236 2003\cd4';

*PARALLEL PLOTS;
goptions cback=white;
axis1 minor=none label=( angle=90 'log2(CD4)' ) order=6 to 12 by 1;
axis2 minor=none label=('Binmonth');
symbol1 i=join v=none r=250  color=black;
title 'Parallel Plot of CD4 Data';
proc sort data=M236.CD4M236 out=sorted;
  by smoker;
run;
proc gplot data=SORTED;
  where smoker ne '';
  plot log2cd4*binmonth=id /vaxis=axis1 haxis=axis2 nolegend;
  by smoker;
run;
quit;
 
proc sort data=m236.cd4m236 out=sortformix;
  by id binmonth;
run;

*PROC MIXED;
proc mixed data=sortformix dfbw noclprint noitprint covtest;
  class id  smoker druguser conversion;
  model log2cd4 = age smoker druguser meancesd conversion
                  conversion*binmonth/s;
  random intercept binmonth / subject=id type=un;
  ods output solutionf=soln;
run;