生存分析的SAS编程操作.ppt

生存分析的SAS编程操作 薛 富 波 2006/05/21 生存分析的有关概念 事件时间（time-to-event，TTE）：又称为生存时间（survival time）或失效时间（failure time），是指从研究开始到相应事件发生的时间。 删失值（censored value）：是指TTE在某些观测上其确切的取值是未知的，所能知道的仅仅是其测量值肯定大于或小于某个特定值（即其取值的下限或上限），或者位于某个特定的取值范围内（取值区间）。 右侧删失：反映真实值下限的数据。 左侧删失：反映真实值上限的数据。 区间删失：反映真实值所在区间的数据。 生存分析的方法分类 非参数法：乘积极限法、寿命表法等； 参数法：指数分布法、威布尔分布（Weibull Distribution）法、对数正态回归分析法以及对数logistic回归分析法等； 半参数法：Cox模型分析方法。 生存分析方法的SAS过程 非参数法：lifetest过程； 参数法：lifereg过程； 半参数法：phreg过程。 生存分析的非参数方法 非参数法生存分析示例1 45例乳腺癌患者中，免疫过氧化物酶检测结果显示9例阳性、36例阴性，比较其生存时间分布的SAS程序如下。 非参数法生存分析示例2 lifetest过程 用以对右侧删失类型的时间数据进行时间分布的描述以及多组间时间分布的比较（检验两个或多个样本是否来自具有同一个生存函数的总体）； 提供两种秩检验方法和一种似然比检验方法来检验多组生存函数的同一性，两种秩检验方法为包括广义Savage（指数得分）检验（即log-rank检验）和Wilcoxon检验。 lifetest过程可包含的语句 proc lifetest语句 proc lifetest语句为调用lifetest过程的开始，其主要语句选项及其可执行的功能见下表。 proc lifetest语句的主要选项 intervals=选项的设置方法 端点值间以空格分隔，如“intervals=10 30 50 70”； 端点之间以逗号分隔，如“intervals=10,30,50,70”； “a to b”的形式，如“intervals=1 to 7”； “a to b by c”的形式，如“intervals=1 to 7 by 2”； 混合方式，如“intervals=1,3 to 7, 10”。 plots=选项的绘图类型 strata语句 strata语句的作用是指定用以分层的变量； lifetest过程根据strata变量的非缺失值对输入数据进行分层； 如果要将strata变量的缺失值也作为一个合法的水平用以分层，可在strata语句中设置“missing” 选项； strata语句中可指定多个分层变量，分层情况将由各分层变量水平的组合来决定； 对于数值型分层变量，若为其设置分层区间端点列表，该变量的水平与形成的区间一一对应，设置方式与proc lifetest语句的“intervals=”选项中类似 ； 对于数值型分层变量，若为其设置分层区间端点列表，该变量的水平与形成的区间一一对应。所划分的区间均为半开半闭区间，第一个区间总是以-∞为下界，最后一个区间总是以∞为上界。 strata语句分层区间端点列表设置方式 ①端点值间以空格分隔，如： ? strata age(5 10 20 30); ②端点之间以逗号分隔，如： strata age(5,10,20,30); ③“a to b”的形式，如： strata age(5 to 10); ④“a to b by c”的形式，如： strata age(5 to 30 by 10); ⑤混合方式，如： strata age(5,10 to 50 by 10); 举例： strata wt(50 to 80 by 10); 此时变量wt将按照如下的5个区间划分成5个水平: test语句 test语句指定需要与生存时间进行关联性分析的协变量； test语句中的变量必须为数值型； 对于test变量，lifetest过程将通过两类秩检验统计量来检验它们与生存时间的关联性； 每一个test变量的单变量检验的统计量也将被给出； 一个有关协变量联合效应的统计量列表也将被给出，此统计量列表的顺序和各协变量对联合效应贡献的大小顺序相一致。 time语句 time语句为调用lifetest过程所必需，用以指定代表生存时间的变量。 time语句具有如下的形式： time variable *censor(list) ; 其中“variable”即为生存时间变量的名称，此变量名称必须指定； 随后的星号“*”及“censor(list)”项可有可无，