The macro (%newsurv) will do them for you if you want to have them with a macro call like:
data BMT;
input DIAGNOSIS Ftime Status Gender@@;
label Ftime="Days";
format Diagnosis grpLabel.;
datalines;
1 2081 0 1 1 1602 0 1
1 1496 0 1 1 1462 0 0
1 1433 0 1 1 1377 0 1
1 1330 0 1 1 996 0 1
1 226 0 0 1 1199 0 1
1 1111 0 1 1 530 0 1
1 1182 0 0 1 1167 0 0
1 418 2 1 1 383 1 1
1 276 2 0 1 104 1 1
1 609 1 1 1 172 2 0
1 487 2 1 1 662 1 1
1 194 2 0 1 230 1 0
1 526 2 1 1 122 2 1
1 129 1 0 1 74 1 1
1 122 1 0 1 86 2 1
1 466 2 1 1 192 1 1
1 109 1 1 1 55 1 0
1 1 2 1 1 107 2 1
1 110 1 0 1 332 2 1
2 2569 0 1 2 2506 0 1
2 2409 0 1 2 2218 0 1
2 1857 0 0 2 1829 0 1
2 1562 0 1 2 1470 0 1
2 1363 0 1 2 1030 0 0
2 860 0 0 2 1258 0 0
2 2246 0 0 2 1870 0 0
2 1799 0 1 2 1709 0 0
2 1674 0 1 2 1568 0 1
2 1527 0 0 2 1324 0 1
2 957 0 1 2 932 0 0
2 847 0 1 2 848 0 1
2 1850 0 0 2 1843 0 0
2 1535 0 0 2 1447 0 0
2 1384 0 0 2 414 2 1
2 2204 2 0 2 1063 2 1
2 481 2 1 2 105 2 1
2 641 2 1 2 390 2 1
2 288 2 1 2 421 1 1
2 79 2 0 2 748 1 1
2 486 1 0 2 48 2 0
2 272 1 0 2 1074 2 1
2 381 1 0 2 10 2 1
2 53 2 0 2 80 2 0
2 35 2 0 2 248 1 1
2 704 2 0 2 211 1 1
2 219 1 1 2 606 1 1
3 2640 0 1 3 2430 0 1
3 2252 0 1 3 2140 0 1
3 2133 0 0 3 1238 0 1
3 1631 0 1 3 2024 0 0
3 1345 0 1 3 1136 0 1
3 845 0 0 3 422 1 0
3 162 2 1 3 84 1 0
3 100 1 1 3 2 2 1
3 47 1 1 3 242 1 1
3 456 1 1 3 268 1 0
3 318 2 0 3 32 1 1
3 467 1 0 3 47 1 1
3 390 1 1 3 183 2 0
3 105 2 1 3 115 1 0
3 164 2 0 3 93 1 0
3 120 1 0 3 80 2 1
3 677 2 1 3 64 1 0
3 168 2 0 3 74 2 0
3 16 2 0 3 157 1 0
3 625 1 0 3 48 1 0
3 273 1 1 3 63 2 1
3 76 1 1 3 113 1 0
3 363 2 1
;
run;
%newsurv(data=bmt,time=ftime,cens=status,class=diagnosis, method=cif,ev_vl=1, summary=0, gpath=~/ibm/,
risklist=0 to 2500 by 500);
The code that runs in the background would look something like this (it calculates the N at risk manually):
proc lifetest data=bmt alpha=0.05 alphaqt=0.05 cifvar error=aalen
outcif=_surv (where=((event>. or t1=0 or c1=1))
rename=(stratum=stratumnum ftime=t1 cif=s1 censored=c1 cif_lcl=lcl1 cif_ucl=ucl1 variance=vcif)
keep=cif ftime censored event atrisk cif_lcl cif_ucl variance stratum)
conftype=LOG;
strata diagnosis;
time ftime * status(0) / failcode=1;
run;
data _splot;
set _surv;
by stratumnum;
array _risktimes_ {6} (0, 500, 1000, 1500, 2000, 2500);
retain _ncens _nevents _lagncens _lagnevents;
if first.stratumnum then do;
_ncens=0;
_nevents=0;
_lagncens=0;
_lagnevents=0;
_count_=1;
end;
if _count_ le dim(_risktimes_) then do;
if t1 <=_risktimes_{_count_} then do;
if event>. then _nevents=_nevents+event;
if c1>. then _ncens=_ncens+c1;
end;
else do;
if event>. then _lagnevents=_lagnevents+event;
if c1>. then _lagncens=_lagncens+c1;
end;
end;
if _count_ le dim(_risktimes_) then do;
if t1 >=_risktimes_(_count_) then do until(x=1);
tatrisk=_risktimes_(_count_);
time=tatrisk;
call missing(event, c1);
output;
_nevents=_nevents+_lagnevents;
_ncens=_lagncens+_ncens;
_lagnevents=0;
_lagncens=0;
_count_+1;
if _count_>dim(_risktimes_) then x=1;
else if t1 < _risktimes_(_count_) then x=1;
else x=0;
end;
if last.stratumnum then do i=_count_ to dim(_risktimes_);
tatrisk=_risktimes_(_count_);
time=tatrisk;
atrisk=0;
call missing(event, c1);
output;
_nevents=_nevents+_lagnevents;
_ncens=_lagncens+_ncens;
_lagnevents=0;
_lagncens=0;
end;
end;
keep tatrisk time stratumnum atrisk event c1 t1 _count_ _nevents _ncens _lagnevents _lagncens;
rename c1=censored;
run;
How to Concatenate Values
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.
-
3770 views
-
2 likes
-
-
