KM and Life Table Estimation
/* Survival Analysis: K-M & Life-Table estimation */ data leukemia; input weeks status drug $; if drug='1' then drug='6-MP'; else drug='placebo'; datalines; 6 0 1 6 1 1 6 1 1 6 1 1 7 1 1 9 0 1 10 0 1 ... 1 1 0 1 1 0 2 1 0 2 1 0 .... ; run;
/* KM estimation */ ods pdf file='\\psf\Home\Desktop\km_lifetable.pdf' style=SASdocprinter; options nodate; title KM and Lifetable for Leukemia Pts: '6-MP' vs Placebo; proc lifetest data=leukemia method=km conftype=loglog plots=survival(cl) plots=(s,ls,lls) graphics outsurv=a; time weeks*status(0); strata drug; symbol1 v=none color=black line=1; symbol2 v=none color=black line=2; run
proc print data=a; run;
/* Life Table estimation */ proc lifetest data=leukemia method=act intervals=10,20,30,40 plots=(s,ls,lls) graphics outsurv=b; time weeks*status(0); strata drug; symbol1 v=none color=black line=1; symbol2 v=none color=black line=2; run;
proc print data=b; run;
title; ods pdf close;
Full data set:
data leukemia; input weeks status drug $; if drug='1' then drug='6-MP'; else drug='placebo'; datalines; 6 0 1 6 1 1 6 1 1 6 1 1 7 1 1 9 0 1 10 0 1 10 1 1 11 0 1 13 1 1 16 1 1 17 0 1 19 0 1 20 0 1 22 1 1 23 1 1 25 0 1 32 0 1 32 0 1 34 0 1 35 0 1 1 1 0 1 1 0 2 1 0 2 1 0 3 1 0 4 1 0 4 1 0 5 1 0 5 1 0 8 1 0 8 1 0 8 1 0 8 1 0 11 1 0 11 1 0 12 1 0 12 1 0 15 1 0 17 1 0 22 1 0 23 1 0 ; run;