文档介绍：生存时间分析的SAS过程一(非参数统计过程)
PROC LEFETEST <选择项>;
TIME 时间变量*结局变量(截尾值);
BY;
FREQ;
STRATA;
选择项说明:
METHOD=PL(或KM)| LT WIDTH=组距 PLOTS=(S)
STRATA语句为分组计算生存率并进行假设检验,统计量为LOG-RANK、WILCOXON
和-2LOG(LR)
小样本生存率的 Kaplan-Meier估计与生存率比较的Log-rank检验
此法直接用概率乘法原理估计生存率,故又称乘积极限法(Product-limit method PL),它由Kaplan-Meier于1958年提出,故又称Kaplan-Meier法(记为KM法)。这是一种非参数法,主要用于小样本资料的生存率估计。
例1 两组淋巴肉瘤患者治疗后复发时间(月数),对照组为“摘除十放疗”,处理组为“摘除十放疗十化疗”,问在“摘除十放疗”基础上附加“化疗”是否可延长缓解期(相当于生存期)?
两组淋巴瘤患者治疗后复发时间(月数,“+”表示未复发)
对照组
2
2
3
3
5+
9
9+
10
12+
处理组
2
5+
10
13
13+
14+
14+
19+
19+
对照组
19+
25
32+
36+
42+
处理组
25+
26+
30+
39+
41+
51+
57+
61+
data sur;
do r=1 to 2;
input n;
do id=1 to n;
input c t@@;
output;
end;
end;
cards;
14
1 2 1 2 1 3 1 3 0 5 1 9 0 9 1 10 0 12 0 19 1 25 0 32 0 36 0 42
17
1 2 0 5 1 10 1 13 0 13 0 14 0 14 0 19 0 19 0 25 0 26 0 30 0 39 0 41 0 51 0 57 0 61
;
Run;
proc lifetest method=pl plots=(s);
time t*c(0);
strata r;
run;
对照组的生存率
处理组的生存率
K-M生存曲线
生存时间分析SAS过程二(参数模型)
PROC LIFEREG <选择项>;
BY 分类变量;
CLASS 分组变量;
MODEL 时间变量*结局变量(截尾值)=自变量/<选择项>;
WEIGHT 加权变量;
关于MODEL语句的部分选择项说明:
| D=WEIBULL|EXPONENTIAL|NORMAL|LOGISTIC

=显著性水平