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

烈酒焚心 提交于 2020-01-07 04:30:39

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

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!