options pagesize=47 linesize=64;
filename grafout 'l09.ps';
goptions device=ps gsfname=grafout gsfmode=append;
data mcac; infile 'mcac.dat';
   input helth life happy status entry time aaids; 
   if helth<0 then delete; one=1;run;
data bigmcac; set mcac; 
   if aaids<time then do; 
      aids=1; output; status=0;
   end; 
   aids=0; output; 
   run;
data bigmcac; set bigmcac; if aids=1 then entry=aaids; else
   time=aaids; run;
title1 'Fit to illustrate residuals';
%include 'l09m.sas';
%csplot(bigmcac,entry,one,one,time,status,
   'Ignoring AIDS status',trycov1=helth);
%csplot(bigmcac,entry,aids,one,time,status,
   'Adjusting for AIDS status',trycov1=helth);
%csplot(bigmcac,entry,one,aids,time,status,
   'Stratifying on AIDS status',trycov1=helth);
%csplot(bigmcac,entry,helth,aids,time,status,
   'Stratifying on AIDS status, adjusting for init. health',
   trycov1=helth,plotscore=1);