計算生存曲線四分位數以及中位數的置信區間:
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 累積發生函數估計