第一章 sas 生存 程序.docVIP

您现在的位置：?生物谷??生物??生物教程??SAS6.12统计教程??生物频道正文 SAS第十课：随访资料的生存分析－－非参数法与半参数Cox比例风险模型 信息来源：医学统计之星　更新时间：2004-7-12 22:20:00 ? 生存分析方法大体上可分为三类：非参数法、参数法和半参数方法，与之相对应，SAS提供了三个程序步用于生存分析，它们是： LIFETEST过程 提供非参数分析方法，用乘积极限法(Product limit method)和寿命表法(Life table method)估计生存率和中位生存时间等；用对数秩检验(Log-rank test)、Wilcoxon检验和似然比检验等做分组比较。该过程主要用于估计生存率及进行单因素分析。 LIFEREG过程 提供指数模型、Weibull模型、Gompertz模型等参数分析方法。 PHREG过程 提供半参数Cox比例风险模型分析。 本章只介绍常用的LIFETEST过程和PHREG过程。 §10.1 LIFETEST过程及其应用 10.1.1 语法格式 LIFETEST过程的语法格式如下： PROC LIFETEST [选项]; TIME 生存时间变量*截尾指示变量(数值) ; [TEST 分组变量名列; STRATA 分组变量名列; FREQ 变量名列; BY 变量名列;] 10.1.2 语法说明 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语句定义需检验的变量，即生存时间与该变量是否有关，如果它后面定义的变量为数值变量，则把该变量当作协变量检验与生存时间的关系。如果它定义的为分组变量，则分组比较生存时间有无差别。 10.1.3 应用实例 例10.1 观察两组卵巢腺癌患者的病程天数如下。请用对数秩检验比较两组的生存期差异有无统计学意义，并作生存率曲线。(医统P.321，7.2题) A组(低恶性高分化癌)：28 29 175 195 309 377＋ 393＋ 421＋ 447＋ 452 709＋ 744＋ 770＋ 1106＋ 1206 B组(高恶性低分化癌)：34 88 137 199 280 291 299＋ 300＋ 309 351 358 369 370 375 382 392 429＋ 451 1119＋ 解：程序如下： data a.yt7_2; 　input t @@; 　if t0 then censor=1; 　　else censor=0; 　if _N_16 then group=A; 　　else group=B; 　t=abs(t); 　cards; 　28 29 175 195 309 -377 ... 　... 382 392 -429 451 -1119 ; proc lifetest method=pl plots=(s,ls); 　time t * censor(1); 　strata group; proc lifetest method=lt plots=(s,h); 　time t * censor(1); 　strata group; run; 本例为不分组资料，生存时间有截尾值，应用乘积极限法和寿命表法分别作生存分析，结局变量只有一个，即生存时间，负值表示截尾，用if语句完