生存分析SAS程序示例

藤菜

1、单因素生存分析
data m;
input t d@@;
cards;
;
proc lifetest plots=(s);
time t*d(0);
run;



2、寿命表法生存分析
data m;
do c=0 to 1;
do t=0 to 15;
input f@@;
output;
end;end;
cards;
456  226 152  171  135  125  83  74 51 42  43  34  18  9  6  0
0 39  22  23  24  107  133  102 68  64  45  53  33  27  23  30
;
proc lifetest method=life width=1 plots=(s);
time t*c(1);
freq f;run;
[注]
频数表按是否截尾值分类，调用freq f;个体编号数据直接用；method=pl默认为乘积极限法





3、单因素两组生存曲线比较
data m;
input t d g@@;
cards;
;
proc lifetest plots=(s);
time t*d(0);
strata g;   
run;