计算生存曲线四分位数以及中位数的置信区间

计算生存曲线四分位数以及中位数的置信区间:
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 累积发生函数估计