文档介绍:SAS第十课:随访资料的生存分析--非参数法与半参数Cox比例风险模型
田晓燕分享| 收藏 
 
生存分析方法大体上可分为三类:非参数法、参数法和半参数方法,与之相对应,SAS提供了三个程序步用于生存分析,它们是:
LIFETEST过程提供非参数分析方法,用乘积极限法(Product limit method)和寿命表法(Life table method)估计生存率和中位生存时间等;用对数秩检验(Log-rank test)、Wilcoxon检验和似然比检验等做分组比较。该过程主要用于估计生存率及进行单因素分析。
LIFEREG过程提供指数模型、Weibull模型、Gompertz模型等参数分析方法。
PHREG过程提供半参数Cox比例风险模型分析。
本章只介绍常用的LIFETEST过程和PHREG过程。
§ LIFETEST过程及其应用
语法格式
LIFETEST过程的语法格式如下:
PROC LIFETEST [选项];
TIME <生存时间变量*截尾指示变量(数值)> ;
[TEST <分组变量名列>;
STRATA <分组变量名列>;
FREQ <变量名列>;
BY <变量名列>;]
语法说明
PROC 和TIME语句为必需的,其他语句都可以省略。TIME语句,为必需语句,定义生存时间和截尾指示变量
【过程选项】
METHOD=方法指定估计生存率所用的方法:
PL 要求用乘积极限法(即Kaplam-Meier法)估计生存率并计算中位生存时间等,为缺省方法。
LT 要求用寿命表法估计生存率等。
INTERVALS=(初值 TO 终值 BY 步长) 只能在指定分析方法为寿命表法时使用。用寿命表法分析时,程序会自动给定生存时间的区间。如果人为规定生存时间的分组区间,则需用该选项指定。步长的缺省值为1。
PLOTS=绘图类型要求输出生存分析图。可供输出的图形有:
S 对生存函数S(t)做图,横、纵坐标分别为t、S(t)。
LS 对-LOGS(t)做图,横、纵坐标分别为t、-LOGS(t)。
LLS 对LOG(-LOGS(t))做图,横、纵坐标分别为LOG(t)、LOG(-LOGS(t))。
H 对风险函数做图, 横、纵坐标分别为t、H(t)。
【TIME语句】
TIME语句用于定义生存时间和截尾指示变量。对截尾指示变量可以指定发生失效事件的数值,默认失效事件用0来表示,截尾事件用1来表示。
【STRATA语句和TEST语句】
STRATA语句定义生存率比较的分组变量,TEST语句定义生存率比较的分组变量或协变量。STRATA语句在这里的作用和BY语句类似,都是要求按分组变量名列进行分析,在计算生存率时各组分开计算;
TEST语句定义需检验的变量,即生存时间与该变量是否有关,如果它后面定义的变量为数值变量,则把该变量当作协变量检验与生存时间的关系。如果它定义的为分组变量,则分组比较生存时间有无差别。
应用实例
观察两组卵巢腺癌患者的病程天数如下。请用对数秩检验比较两组的生存期差异有无统计学意义,并作生存率曲线。(,)
A组(低恶性高分化癌):28 29 175 195 309 377+ 393