filename grafout 'l13.ps';
goptions device=pslmono gsfname=grafout gsfmode=append;
/* Example of competing risks */
data bmt; infile 'bmt.dat';
   input g t1 t2 status relapse nodx tagvh agvhs tcgvh cgvhs
      tplat platstat patage donage patsex donsex patcmv 
      doncmv transwt fab hosp mtx; run;
data bmt; set bmt; keep time which; which=0;
   if agvhs=1 or cgvhs=1 then do;
      if tagvh<tcgvh then do;
         time=tagvh; which=1; end;
      else do;
         time=tcgvh; which=2; end; 
      end;
    else time=t1; run;
/***********************************************************/
/* I can't to the calculations below in one go using by or */
/* stratify, since they use different censoring  values.   */
/***********************************************************/
proc lifetest data=bmt outsurv=g1; time time*which(0,1); run;
proc lifetest data=bmt outsurv=g2; time time*which(0,2); run;
proc lifetest data=bmt outsurv=g3; time time*which(0); run;
proc print data=g1; run;
data g1; set g1; s1=survival; if survival=. then delete; 
   gp=1; run;
data g2; set g2; s2=survival; if survival=. then delete; 
   gp=2; run;
data g3; set g3; s3=survival; if survival=. then delete; 
   gp=3; run;
data g; merge g1 g2 g3; by time; retain sa1 sa2 sa3; 
   keep sa1 sa2 sa3 time gp;
   if s1=. then; else sa1=s1; 
   if s2=. then; else sa2=s2; 
   if s3=. then; else sa3=s3; run;
proc gplot data=g; symbol1 i=stepjl line=1; 
   symbol2 i=stepjl line=2; symbol3 i=stepjl line=3;
   plot sa1*time=1 sa2*time=2 sa3*time=3 /overlay; run;