/* PubH 7450 Exampel 4.1c: Pointwise 95% CI of a survival curve */ options ls=80 center pagesize=60 number label; data cancer; infile '/home/merganser/weip/public_html/course/data/leukemia.dat' firstobs=2; input group $ weeks relaps ; data grp6MP; set cancer; if group='6-MP'; /* (1) obtain point-wise 95% (linear) CI: proc lifetest method=km data=grp6MP outsurv=a; by group; time weeks*relaps(0) ; proc print data=a; run; /* (2) to obtain other types of pointwise CI based on KM estimator, use SAS macro %surv available from http://mayoresearch.mayo.edu/mayo/research/biostat/sasmacros.cfm The key parameter is cl: cl=1: linear transformation cl=3: log trans cl=5: log(-log) trans cl=7: logit trans */ %include 'surv.sas'; %surv(time=weeks, event=relaps, cen_vl=0, data=grp6MP, cl=5, out=survout) proc print data=survout; run; ********************************Output: **************(1): Obs group weeks _CENSOR_ SURVIVAL SDF_LCL SDF_UCL 1 6-MP 0 0 1.00000 1.00000 1.00000 2 6-MP 6 0 0.85714 0.70748 1.00000 3 6-MP 6 1 0.85714 . . 4 6-MP 7 0 0.80672 0.63633 0.97711 5 6-MP 9 1 0.80672 . . 6 6-MP 10 0 0.75294 0.56410 0.94178 7 6-MP 10 1 0.75294 . . 8 6-MP 11 1 0.75294 . . 9 6-MP 13 0 0.69020 0.48084 0.89955 10 6-MP 16 0 0.62745 0.40391 0.85099 11 6-MP 17 1 0.62745 . . 12 6-MP 19 1 0.62745 . . 13 6-MP 20 1 0.62745 . . 14 6-MP 22 0 0.53782 0.28648 0.78915 15 6-MP 23 0 0.44818 0.18439 0.71197 16 6-MP 25 1 . . . 17 6-MP 32 1 . . . 18 6-MP 32 1 . . . 19 6-MP 34 1 . . . 20 6-MP 35 1 . . . **************************(2): p l u n c o o p n c c u i w p w e e n u m n e e e v n r m _ t r r O e e s i _ c f _ _ b k p n o s e e l p s c c s s t t r k v n g e e l l 1 0 1.00000 0 0 21 0 0 . 0.00000 0.00000 1.00000 1.00000 2 6 0.85714 3 1 21 3 1 . 0.14286 0.07636 0.61972 0.95155 3 7 0.80672 1 0 17 4 1 . 0.19328 0.08694 0.56315 0.92281 4 9 0.80672 0 1 16 4 2 . 0.19328 0.08694 0.56315 0.92281 5 10 0.75294 1 1 15 5 3 . 0.24706 0.09635 0.50320 0.88936 6 11 0.75294 0 1 13 5 4 . 0.24706 0.09635 0.50320 0.88936 7 13 0.69020 1 0 12 6 4 . 0.30980 0.10681 0.43161 0.84907 8 16 0.62745 1 0 11 7 4 . 0.37255 0.11405 0.36751 0.80491 9 17 0.62745 0 1 10 7 5 . 0.37255 0.11405 0.36751 0.80491 10 19 0.62745 0 1 9 7 6 . 0.37255 0.11405 0.36751 0.80491 11 20 0.62745 0 1 8 7 7 . 0.37255 0.11405 0.36751 0.80491 12 22 0.53782 1 0 7 8 7 . 0.46218 0.12823 0.26778 0.74679 13 23 0.44818 1 0 6 9 7 . 0.55182 0.13459 0.18805 0.68014 14 25 0.44818 0 1 5 9 8 . 0.55182 0.13459 0.18805 0.68014 15 32 0.44818 0 2 4 9 10 . 0.55182 0.13459 0.18805 0.68014 16 34 0.44818 0 1 2 9 11 . 0.55182 0.13459 0.18805 0.68014 17 35 0.44818 0 1 1 9 12 . 0.55182 0.13459 0.18805 0.68014