计算生存曲线四分位数以及中位数的置信区间:
ods listing close;
ods results off;
proc lifetest data=adtte method=KM plots=(s);
time aval*cnsr(1);
strata trtan;
ods output Quartiles= part3 HomTests=HomTests;
run;
ods results on;
ods listing;
adtte是datain,part3是dataout,Quartiles=表示只输出包含四分位数以及置信区间的那部分;默认的是95%置信区间,alpha=。
strata分层分析,如果只有一层,则不需要使用HomTests=HomTests选项。
cnsr(1),表示在cnsr=1的时候,该条记录删失。
四分位数表示当有百分之25/50/75的人发生事件(cnsr=0)的时候,发生事件的时间(aval)的点估计;95%的置信区间上下限也是时间。四分位的分母是指所有人,并且随着时间的推移,不停的去掉删失的人。
ods listing close;
ods results off;
proc lifetest data=adtte plots=(survival(atrisk) logsurv) timelist=78 162;
time aval*cnsr(1)/eventcode=0;
ods output cif=part3;
run;
ods results on;
ods listing;
timelist=78 162表示在78天,162天的时候,event事件的发生率
/eventcode=0 除了算出时间点上事件的发生率之外,还可以得到发生率的95%置信区间
cif=part3 仅输出发生率和置信区间那块内容到dataset part3中,cif意为cumulative incidence function estimates 累积发生函数估计
来源:https://www.cnblogs.com/jwcutey/p/12159324.html