步骤 Sig. .000 2 a. 起始块编号 0，最初的对数似然函数:-2 倍对数似然值: 201.994 b. 起始块编号 1. 方法 = 向前逐步（似然比）
-2 倍对数 似然值 182.777
整体 (得分) 卡方 df 17.594 2
从上一块开始更改 卡方 df Sig. 19.217 2 .000
X1 62 40 50 33 57 48 28 54 35 47 49 43 48 44 60 40 32 44 48 72 42 63 55 39 44 42 74 61 45 38 62
X2 0 1 1 1 1 1 0 1 0 0 1 0 1 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0
实例分析
例1：为了比较不同手术方法治疗肾上腺肿瘤的疗效， 某研究者随机将43例病人分成两组，甲组23例、乙组20 例的生存时间（月）如下所示：
其中有“+”者是删失数据，表示病人仍生存或失访，括号内为死亡人数。
（1）计算甲、乙两法术后10月的生存率和标准误。 （2）估计两组的中位生存期。 （3）绘制各组生存函数曲线。 （4）比较两组的总体生存时间分布有无差别。
Y 0 1 1 0 0 0 1 1 1 1 0 0 1 1 0 1 1 1 0 1 1 0 1 1 0 0 1 1 0 1 1
一、建立数据文件（data-03.sav） 二、操作过程
主菜单：分析Analyze生存SurvivalCox Regression
Cox回归主对话框
二、操作过程
主对话框参数设置：

2-1.模型检验（全变量模型）
a 模 型系 数 的 综合 测 试 ,b
整体 (得分) 从上一步骤开始更改 -2 倍对数 似然值 卡方 df Sig. 卡方 df Sig. 180.052 20.174 6 .003 21.942 6 .001 a. 起始块编号 0，最初的对数似然函数:-2 倍对数似然值: 201.994 b. 起始块编号 1. 方法 = 输入
二、操作过程
2）
√
水平间的两两比较。
6. 单击Save按钮，弹出保存新变量Save new variables 对话框：
√ √
三、主要输出结果
1. 生存表： 略 2. 两组的中位生存期估计：
3. 绘制生存曲线：
4. 两组生存时间分布的比较：
Company Logo
Cox回归过程
Cox回归过程用于： 1. 多个因素对生存时间的影响作用分析和比较
Company Logo
一、建立数据文件（同前）
二、操作过程
主菜单：分析Analyze生存SurvivalKaplan-Meier
对话框参数设置： 1. 时间time框：选入 “t”。 2. 状态status框：选入“status”，击define events 钮，在single value 框填入“1”。 3. 因子factor框：选入“group”。 4. 单击选项option按钮，弹出对话框： 1）统计量： √生存分析表，系统默认。 √ 均值和中位生存时间，系统默认。 2）图： √生存函数 5. 单击比较因子Compare Factor按钮，弹出对话框： 1）检验统计量Test Statistics： 都用于检验时间分布是否相同。 √对数秩Log-rank：各时间点的权重一样。 Breslow：按各时间点的观察例数赋权。 Tarone-Ware：按各时间点观察例数的平方根赋权。
X6 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 1 0 0 0 0
t 120 40 26 120 120 120 3 120 7 18 120 120 15 4 120 16 24 19 120 24 2 120 12 5 120 120 7 40 108 24 16
Y 0 0 1 0 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 1 1 0 0 0 0 0 0
No
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63
生存分析SPSS过程
（SPSS of Survival Analysis）
邹莉玲 预防医学教研室
生存分析的理论复习
1. 何为生存分析？
生存分析（survival analysis）是将事件的结果（终点事件） 和出现结果经历的时间结合起来分析的一种统计分析方法。
2. 生存分析的目的：
(1)描述生存过程：估计不同时间的总体生存率，计算中位生存期， 绘制生存函数曲线。统计方法包括Kaplan-Meier（K-M）法、 寿命表法。 (2)比较：比较不同处理组的生存率，如比较不同疗法治疗脑瘤的 生存率，以了解哪种治疗方案较优。统计方法log-rank检验等。 (3)影响因素分析：研究某个或某些因素对生存率或生存时间的影 响作用。如为改善脑瘤病人的预后，应了解影响病人预后的主 要因素，包括病人的年龄、性别、病程、肿瘤分期、治疗方案 等。统计方法cox比例风险回归模型等。 (4)预测：建立cox回归预测模型。
2. 生存（或死亡）风险预测
实例分析

例3：为探讨某恶性肿瘤的预后，某研究者收集了63 例患者的生存时间、生存结局及影响因素。影响因素 包括病人年龄、性别、组织学类型、治疗方式、淋巴 结转移、肿瘤浸润程度，生存时间以月计算。变量的 赋值和所收集的资料分别见表17-8和表17-9。试用 Cox回归模型进行分析。
Company Logo
二、操作过程
主菜单：分析Analyze生存Survival寿命表Life tables
对话框参数设置： 1. 时间time框：选入 “t”。 2. 显示时间间隔Display time intervals框：步长by前面填入最大生存时 间的上限（必须包括生存时间最大值），步长by后面填入生存时 间的组距。本例上限填“60”，组距填“1”。 3. 状态status框：选入“status”，击define events 钮，在single value 框填入“1” 4. 因子factor框：选入“group”，定义最小值“1”，最大值“2”。 5. 单击选项option按钮，弹出对话框： 1）√寿命表，系统默认。 2）图： √生存函数 3）比较第一个因子的水平： √整体比较
X3 0 1 0 1 1 0 0 0 1 0 0 1 1 0 1 0 1 0 0 1 0 0 1 0 0 1 0 1 0 1 0
X4 0 1 0 0 1 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 1 1 0 0 1 0 0 1 0 0
X5 1 0 1 0 0 1 1 1 1 1 1 0 0 1 0 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1
三、主要输出结果
பைடு நூலகம்
1.分析例数描述
案 例 处理 摘 要 N 分析 中可 用的 案例 删除 的案 例 事件a 删失 合计 带有缺失值的案例 带有负时间的案例 层中的最早事件之 前删失的案例 合计 0 合计 a. 因变量: t 63 .0% 100.0% 26 37 63 0 0 0 百分比 41.3% 58.7% 100.0% .0% .0% .0%
X2 0 0 0 1 0 0 0 0 1 0 0 1 1 0 1 1 1 1 0 1 0 0 1 0 1 0 1 1 1 0 1 1
X3 0 1 0 1 1 1 1 0 1 1 1 0 1 1 1 1 1 0 1 1 0 0 0 0 0 0 1 1 1 1 0 0
X4 1 0 0 1 0 0 0 0 1 0 0 1 1 0 1 1 1 1 0 1 0 0 1 1 0 0 1 1 1 0 1 1
生存分析(Survival Analysis)菜单
寿命表（Life Tables）过程
Life tables 过程用于（小样本和大样本资料）：
1. 估计某生存时间的生存率，以及中位生存时间。
2. 绘制各种曲线：如生存函数、风险函数曲线等。 3. 对某一研究因素不同水平的生存时间分布的比较。 4. 控制另一个因素后对研究因素不同水平的生存时间分 布的比较。 5. 对多组生存时间分布进行两两比较。 （比较总体生存时间分布采用wilcoxon检验）
Company Logo
一、建立数据文件（data-01.sav）
定义5个变量：
生存时间变量：t，值标签“生存时间（月）”
生存状态变量 ：status，取值“1=死亡，0=删失或存活” 频数变量：freq，值标签“人数” 分组变量：group，取值“1=甲组，2=乙组” 生存时间序号变量(可无）：i
（2）采用逐步回归法进行Cox模型分析的结果提示：模型拟合自变量进入和 剔除的检验水准分别为0.05和0.1时，筛选后的最佳模型包含两个协变量， 即X4（治疗方式）和X5（淋巴结是否转移），该拟合模型总体检验提示
具有统计学意义（整体卡方=17.594，P<0.001）。

3. 参数估计（逐步回归模型，Method=向前法 LR，自变量进入P≤0.05，剔除P>0.10 ）
X5 1 0 1 1 0 0 0 1 0 0 0 1 0 0 0 0 0 1 0 0 1 1 1 1 1 1 1 0 0 0 1 1
X6 0 0 1 0 0 0 0 1 1 1 0 0 0 1 1 1 0 1 0 1 0 0 0 0 0 0 1 0 0 1 0 0
t 52 51 35 103 7 60 58 29 70 67 66 87 85 82 76 74 63 101 100 66 93 24 93 90 15 3 87 120 120 120 120 120
实例分析
例2：（数据同例1）为了比较不同手术方法治疗肾上腺 肿瘤的疗效，某研究者随机将43例病人分成两组，甲组 23例、乙组20例的生存时间（月）如下所示：