生存分析SPSS单因素和多因素对生存率的可能分析

1、生存分析生存分析spss过程过程（spss of survival analysis） 邹莉玲邹莉玲预防医学教研室预防医学教研室1. 1. 何为生存分析？何为生存分析？生存分析生存分析（survival analysis）是将是将事件的结果事件的结果（终点事（终点事件）和出现结果经历的件）和出现结果经历的时间时间结合起来分析的一种统计分析方法。结合起来分析的一种统计分析方法。2. 2. 生存分析的目的：生存分析的目的：(1)(1)描述生存过程描述生存过程：估计不同时间的总体生存率，计算中位生存期，估计不同时间的总体生存率，计算中位生存期，绘制生存函数曲线。统计方法包括绘制生存函数曲线。统计方法

2、包括kaplan-meierkaplan-meier（k-mk-m）法、寿命）法、寿命表法。表法。(2)(2)比较：比较：比较不同处理组的生存率，如比较不同疗法治疗脑瘤的比较不同处理组的生存率，如比较不同疗法治疗脑瘤的生存率，以了解哪种治疗方案较优。统计方法生存率，以了解哪种治疗方案较优。统计方法log-ranklog-rank检验等。检验等。(3)(3)影响因素分析：影响因素分析：研究某个或某些因素对生存率或生存时间的影研究某个或某些因素对生存率或生存时间的影响作用。如为改善脑瘤病人的预后，应了解影响病人预后的主要响作用。如为改善脑瘤病人的预后，应了解影响病人预后的主要因素，包括病人的年龄、

3、性别、病程、肿瘤分期、治疗方案等。因素，包括病人的年龄、性别、病程、肿瘤分期、治疗方案等。统计方法统计方法coxcox比例风险回归模型等。比例风险回归模型等。(4)(4)预测：预测：建立建立coxcox回归预测模型。回归预测模型。生存分析的理论复习生存分析的理论复习生存分析生存分析( (survival analysis) )菜单菜单寿命表寿命表（life tables）过程过程life tables 过程用于（小样本和大样本资料）：1.估计某生存时间的生存率，以及中位生存时间。2.绘制各种曲线：如生存函数、风险函数曲线等。3.对某一研究因素不同水平的生存时间分布的比较。4.控制另一个因素后对

4、研究因素不同水平的生存时间分布的比较。5.对多组生存时间分布进行两两比较。（比较总体生存时间分布采用wilcoxon检验）company logo实例分析实例分析例1：为了比较不同手术方法治疗肾上腺肿瘤的疗效，某研究者随机将43例病人分成两组，甲组23例、乙组20例的生存时间（月）如下所示：其中有“+”者是删失数据，表示病人仍生存或失访，括号内为死亡人数。（1）计算甲、乙两法术后10月的生存率和标准误。（2）估计两组的中位生存期。（3）绘制各组生存函数曲线。（4）比较两组的总体生存时间分布有无差别。company logo一、建立数据文件（一、建立数据文件（data-01.sav）定义定义5个

5、变量：个变量：生存时间变量：t，值标签“生存时间（月）”生存状态变量 ：status，取值“1=死亡，0=删失或存活”频数变量：freq，值标签“人数”分组变量：group，取值“1=甲组，2=乙组”生存时间序号变量(可无）：i二、操作过程二、操作过程主菜单：分析主菜单：分析analyze生存生存survival寿命表寿命表life tables对话框参数设置：对话框参数设置：1.时间time框：选入 “t”。2.显示时间间隔display time intervals框：步长by前面填入最大生存时间的上限（必须包括生存时间最大值），步长by后面填入生存时间的组距。本例上限填“60”，组距填“

6、1”。3.状态status框：选入“status”，击define events 钮，在single value框填入“1”4.因子factor框：选入“group”，定义最小值“1”，最大值“2”。5.单击选项option按钮，弹出对话框： 1）寿命表，系统默认。 2）图： 生存函数 3）比较第一个因子的水平： 整体比较三、主要输出结果三、主要输出结果1.10月生存率的估计：月生存率的估计： 甲法甲法 48%，标准误，标准误 0.1 乙法乙法 30%，标准误，标准误 0.1 2.两组的中位生存期估计：两组的中位生存期估计：company logo3. 绘制生存曲线：绘制生存曲线：compan

7、y logo4. 两组生存时间分布的比较：两组生存时间分布的比较：kaplan-meier 过程过程kaplan-meier过程用于(尤其小样本资料）：1.估计各生存时间的生存率以及中位生存时间。2.绘制各种曲线：如生存函数、风险函数曲线等。3.比较某研究因素不同水平的生存时间有无差异。4.控制某个分层因素后对研究因素不同水平的生存时间分布进行比较。5.对多组生存时间分布进行两两比较。（各总体分布比较采用log-rank等非参数方法）company logo实例分析实例分析例2：（数据同例1）为了比较不同手术方法治疗肾上腺肿瘤的疗效，某研究者随机将43例病人分成两组，甲组23例、乙组20例的生

8、存时间（月）如下所示：其中有“+”者是删失数据，表示病人仍生存或失访，括号内为死亡人数。（1）计算甲、乙两法各生存时间的生存率和标准误。（2）估计两组的中位生存期。（3）绘制各组生存函数曲线。（4）比较两组的总体生存时间分布有无差别。一、建立数据文件（同前）一、建立数据文件（同前）二、操作过程二、操作过程主菜单：分析主菜单：分析analyze生存生存survivalkaplan-meier对话框参数设置：对话框参数设置：1.时间time框：选入 “t”。2.状态status框：选入“status”，击define events 钮，在single value框填入“1”。3.因子factor框

9、：选入“group”。4.单击选项option按钮，弹出对话框： 1）统计量： 生存分析表，系统默认。 均值和中位生存时间，系统默认。 2）图： 生存函数5. 单击比较因子compare factor按钮，弹出对话框： 1）检验统计量test statistics： 都用于检验时间分布是否相同。 对数秩log-rank：各时间点的权重一样。 breslow：按各时间点的观察例数赋权。 tarone-ware：按各时间点观察例数的平方根赋权。二、操作过程二、操作过程2）水平间的两两比较。水平间的两两比较。6. 单击save按钮，弹出保存新变量save new variables 对话框：三、主要

10、输出结果三、主要输出结果1.生存表：生存表： 略略 2.两组的中位生存期估计：两组的中位生存期估计：3. 绘制生存曲线：绘制生存曲线：company logo4. 两组生存时间分布的比较：两组生存时间分布的比较：cox回归回归过程过程cox回归过程用于： 1. 多个因素对生存时间的影响作用分析和比较 2. 生存（或死亡）风险预测实例分析实例分析o例3：为探讨某恶性肿瘤的预后，某研究者收集了63例患者的生存时间、生存结局及影响因素。影响因素包括病人年龄、性别、组织学类型、治疗方式、淋巴结转移、肿瘤浸润程度，生存时间以月计算。变量的赋值和所收集的资料分别见表17-8和表17-9。试用cox回归模型

11、进行分析。no1234567891011121314151617181920212223242526272829303132x15457584348404436394242425155495248543840381967374349505332464344x200010000100110111101001010111011x301011110111011111011000000111100x410010000100110111101001100111011x510110001000100000100111111100011x600100001110001110101000000100100t5

12、251351037605829706766878582767463101100669324939015387120120120120120y00101001000000000001010011000000no33343536373839404142434445464748495051525354555657585960616263x162405033574828543547494348446040324448724263553944427461453862x20111110100101010001001000100100x30101100010011010100100100101010x401

13、00100100100010000001100100100x51010011111100101011111011011101x60100000010000000110000000010000t120 40 26120120120 3120 7 18120120 15 4120 16 24 19120 24 2120 12 5120120 7 40108 24 16y0110001111001101110110110011011表表17-9 63名某恶性肿瘤患者的生存时间（月）及影响因素名某恶性肿瘤患者的生存时间（月）及影响因素一、建立数据文件（一、建立数据文件（data-03.sav）二、操作

14、过程二、操作过程主菜单：分析主菜单：分析analyze生存生存survivalcox regressioncox回归主对话框回归主对话框二、操作过程二、操作过程主对话框参数设置：主对话框参数设置：1.时间time框：选入 “ t ”。2.状态status框：选入“ y ”，击define events 钮，在single value框填入“1”3.协变量covariates框：选入“x1x6”。4.方法method框：选择自变量进入cox模型的方法，spss提供以下7种方法。o1.分析例数描述案例处理摘要案例处理摘要2641.3%3758.7%63100.0%0.0%0.0%0.0%0.0%6

15、3100.0%事件a删失合计分析中可用的案例带有缺失值的案例带有负时间的案例层中的最早事件之前删失的案例合计删除的案例合计n百分比因变量: ta. 三、主要输出结果三、主要输出结果o2-1.模型检验（全变量模型）模模 型型 系系 数数 的的 综综 合合 测测 试试a a, ,b b180.05220.1746.00321.9426.00121.9426.001-2 倍对数似然值卡方dfsig.整体 (得分)卡方dfsig.从上一步骤开始更改卡方dfsig.从上一块开始更改起始块编号 0，最初的对数似然函数:-2 倍对数似然值: 201.994a. 起始块编号 1. 方法 = 输入b. 结果提示

16、：结果提示：（1）对模型总体检验有显著意义（）对模型总体检验有显著意义（p=0.003），即至少），即至少有有1个自变量的总体回归系数不为个自变量的总体回归系数不为0。o2-2.模型检验（逐步回归模型，method=向前法lr，自变量进入p0.05，剔除p0.10 ）模模 型型 系系 数数 的的 综综 合合 测测 试试a a, ,b b182.77717.5942.00019.2172.000步骤2-2 倍对数似然值卡方dfsig.整体 (得分)卡方dfsig.从上一块开始更改起始块编号 0，最初的对数似然函数:-2 倍对数似然值: 201.994a. 起始块编号 1. 方法 = 向前逐步（似

17、然比）b. （2）采用逐步回归法进行）采用逐步回归法进行cox模型分析的结果提示：模型拟合自变量进入和模型分析的结果提示：模型拟合自变量进入和剔除的检验水准分别为剔除的检验水准分别为0.05和和0.1时，筛选后的最佳模型包含两个协变量，时，筛选后的最佳模型包含两个协变量，即即x4（治疗方式）和（治疗方式）和x5（淋巴结是否转移），该拟合模型总体检验提示（淋巴结是否转移），该拟合模型总体检验提示具有统计学意义（整体卡方具有统计学意义（整体卡方=17.594，p0.10 ）方方程程中中的的变变量量-1.762.54810.3371.001.172.059.503.931.4454.3891.036

18、2.5381.0626.066x4x5步骤2bsewalddfsig.exp(b)下部上部95.0% ci 用于 exp(b)b ：偏回归系数 ，se：偏回归系数的标准误wald：用于检验总体偏回归系数与0有无显著差异， v=1时，w=（b/se）2 。exp(b)：相对危险度估计（rr值）（3）x4（治疗方式）对生存时间有影响，采用新疗法病人的死亡风险降至传（治疗方式）对生存时间有影响，采用新疗法病人的死亡风险降至传统疗法的统疗法的17.2%（rr的的95%ci为为0.0590.503）。）。（4）x5（淋巴结是否转移）对生存时间也有影响，有淋巴结转移病人的死亡（淋巴结是否转移）对生存时间也有影响，有淋巴结转移病人的死亡风险为无淋巴结转移病人的风险为无淋巴结转移病人的2.538倍（倍（ rr的的95%ci为为1.0626.066）。）。o4. 自/协变量的均值（选入2个变量的模型）协变量均值协变量均值.413.571x4x5均值o5. 生存曲线（选入2个变量的模型）company logop 练习练习1：某临床试验对20名第或第期黑色素瘤的患者进行随访研究，截至研究期结束，记录的生存资料见表1。应用spss软件建立数据文件，并计算100周生存率。表1 20名第或期黑色素瘤患者的治疗后生存时间（周）12.877.215.682.424.0+87.2+26.494