寿命表示法在SAS统计系统中的应用

在SAS 统计系统中，PROC LIFETEST 过程提供有非参数分析方法，可以用乘积极限法和寿命表法估计生存率和中位生存时间等，用对数秩检验（Log-rank test）、Wilcoxon 检验和似然比检验等做分组比较。该过程主要用于估计生存率及进行单因素分析。

（一）调用格式

PROC LIFETEST 选项；

TIME 生存时间变量∗截尾指示变量（数字）；

TEST 分组变量名列表；

STRATA 分组变量名列表；

FREQ 变量名列表；

BY 变量名列表。

（二）语句说明

（1）PROC LIFETEST 语句

①DATA=SAS 数据集，规定PROC FREQ 语句使用的数据集。

②METHOD=SAS 方法，指定估计生存率所用的方法。

a.PL，要求用乘积极限法（即Kaplam-Meier 法）估计生存率并计算中位生存时间等，为缺省方法。

b.LT，要求用寿命表法估计生存率等。

c.INTERVALS=（初始值TO 终止值BY 步长）只能在指定分析方法为寿命表法时使用，用寿命表法分析时，程序会自动给定生存时间的区别。如果规定生存时间的分组区间，则需要指定该选项。步长的缺省值为1。

③WIDTH=宽度，指定用LT 法的生存时间区间的宽度。

④PLOTS=绘图类型，要求输出生存分析图。可供输出的图形有：

a.S，对生存函数S（t）作图，横、纵坐标分别为t、S（t）。

b.LS，对-LOG S（t）作图，横、纵坐标分别为t、-LOG S（t）。

c.LLS，对LOG（-LOG S（t））作图，横、纵坐标分别为-LOG S（t）和LOG（-LOG S（t））。

d.H，对风险函数作图，横、纵坐标分别为t、H（t）。

e.NOTABLE，指令不输出生存函数估计结果，只输出生存时间的截尾数据和完全数据的个数以及散点图和检验结果。

（2）TIME 语句

用于定义生存时间和截尾指示变量。对截尾指示变量可以指定发生失效时间的数值，默认失效事件用0 来表示，截尾事件用1 来表示。

（3）STRATA 语句

定义生存率比较的分组变量，TEST 语句定义生存率比较的分组变量或协变量。STRATA语句在这里的作用和BY 语句类似，都是要求按分量变量名列进行分析，在计算生存率时各组分开计算。

（4）TEST 语句

定义需检验的变量，即生存时间与该变量是否有关，如果它后面定义的变量为数值变量，则把该变量当作协变量检验与生存时间的关系。

（三）实例应用

例11.2　某兽医院检测某抗生素对犬某疾病进行治疗，治疗后收集35 只犬疾病治疗后的生存时间为t（周），培养观察后发现这些病犬中有一部分出现了白细胞（WBC）增加的现象。先将35 只犬按是否出现WBC 倍增分成两组如下（注：负值代表缺失数据），使用寿命表法分析白细胞有无倍增对其生存时间长短有无显著影响。(www.xing528.com)

A 组（倍增）：2 -2.5 3.5 4 4 -5 6 -6 7 -7 8 -9 10.5 12.5 19；

B 组（未增）：2.5 5 7 -8.5 9 -10 11 -11 12 13 -14 15 -16 17 -18 19 -20 21 24 32。

SAS 程序如图11.4 所示。

图11.4

输出结果为：

①本部分结果显示LT 法对第1 组资料进行统计描述。Life table 表给出了生存时间区间的下限和上限、死亡数、截尾数、有效样本大小、条件死亡概率、条件死亡概率的标准误、生存概率、死亡概率、死亡概率的标准误、中位数剩余生存寿命（即在时刻ti或者犬有一般可望生存到的时间）、中位数的标准误、区间中点概率密度函数的估计值PDF、PDF 的标准误、风险函数、风险函数的标准误等。

SAS 程序如图11.5 所示。

图11.5　例11.2SAS 程序图（1）

②本部分结果显示LT 法对第2 组资料进行统计描述的结果，同①，在此不再赘述。

SAS 程序如图11.6 所示。

图11.6　例11.2SAS 程序图（2）

③两组病犬的总数、死亡数、截尾数和截尾百分比。如倍增组总数为15，死亡数为10，截尾数为5，截尾百分比为33.33%。

SAS 程序如图11.7 所示。

图11.7　例11.2SAS 程序图（3）

④两组病犬的生存分布函数曲线。A 组（倍增）在第20 个周终止，B 组（未增）在第35 个周终止，A 组生存率下降速度显著快于B 组，说明B 组病犬比A 组病犬的生存时间长。

SAS 程序如图11.8 所示。

图11.8　例11.2SAS 程序图（4）

⑤给出生存时间T 对应的死亡风险函数。A 组死亡风险率较高，尤其是在0 ～18 周内，B组死亡风险率在0 ～22 周较为平缓，增长不显著，在23 ～28 周有下降的趋势，但在28 周以后死亡风险率直线上升。

SAS 程序如图11.9 所示。

图11.9　例11.2SAS 程序图（5）

⑥各组生存函数曲线齐性检验。依次给出秩次统计量、Log-rank 统计量的协方差矩阵、Wilcoxon 统计量的协方差矩阵、各组生存函数一致性检验结果等。本实验结果表明，此资料不服从指数分布，近似服从Weibull 分布，故选用Log-rank 法或Wilcoxon 法检验的结果，两条生存曲线分布有显著性差异（P<0.05）。因此，无WBC 倍增的病犬时间显著长于有WBC 倍增的病犬。

SAS 程序如图11.10 所示。

图11.10　例11.2SAS 程序图（6）