5# fxf258
proc lifetest data= fxf outsurv=plotdata;*输出生存数据集;
time qtime*censor(0);
strata group;
run;
proc sql;
create table plotdata_ as
select *
, ifn(censor=0 & missing(survival), min(survival), survival) as survival, *如果最后一个观测是censored,那么survival=最小的survival;
1-survival as prob
from plotdata group by group order by group,qtime;
quit;
symbol i=stepj;
proc gplot data=plotdata_;
plot qtime*prob=group;
run;
*没有试运行,如果有错误再讨论;