SPSS生存分析过程.docx

1、SPSS生存分析过程SPSS Survival(生存分析)菜单SPSS Survival菜单包括Life Tables过程、Kaplan-Meier过程、Cox Regression过程、Cox w/Time-Dep Cov过程。这里只介绍Life Tables过程和Kaplan-Meier过程。Life Tables过程Life Tables过程用于：1、估计某生存时间的生存率。2、绘制各种曲线如生存函数、风险函数曲线等。3、对某一研究因素不同水平的生存时间分布进行比较，控制另一因素后对研究因素不同水平的生存时间分布进行比较，包括从总体上比较和不同水平之间进行两两比较。一、建立数据文件 定义

2、两个列变量：时间变量： 取名“time”，label标上“survival time(week)”。生存状态变量：取名“status”，并赋值：0=“删失”，1=“死亡”。二、操作过程 从菜单选择1、Analyze=Survival =Life Tables2、Time框：选入time3、Display Time Intervals框：在by前面的框内填入生存时间上限，本例填入20（此区间必须包括生存时间的最大值）；在by后面的框内填入生存时间的组距，本例填入5，以保证结果列出“15-”的组段。4、Status框：选入status；击define events钮,在single value框右

3、边的空格中输入15、单击Option按钮，弹出对话框：Life Table(s) 输出寿命表，系统默认Plots: 选Survival（累积生存函数曲线）击Continue6、单击OK钮附：界面说明图1 寿命表主对话框【Time】框选入生存时间变量。【Display Time Intervals】框欲输出生存时间范围及组距。在by前面的框内填入生存时间上限，本例填入200（此区间必须包括生存时间的最大值）；在by后面的框内填入生存时间的组距，本例填入20，以保证结果列出“100-”的组段。【Status】框选入生存状态变量，并定义终结事件的标记值。选入变量“Status”后，【Define E

4、vent】钮被激活变黑，击该按钮，弹出定义终结事件标记值的对话框（图1）。对二分类变量，一般以死亡、复发、恶化等表示终结事件。本例以死亡为终结事件，其标记值为1，故在Single value 框内填入1。击Continue钮。若生存状态变量取值为一连续型变量，如反应变量为收缩压，则在下面的Range of values 框中输入140 through 400，此处上限输入400是我随便输入的一个上限，目的是为了定义高血压患者，实际上恐怕没有人的血压能达到400，这样才能保证包括所有的高血压病例，具体情况具体分析。图2 定义终结事件标记值的对话框【Factor】框定义第1层因素，即分组因素。【B

5、y Factor】框定义第2层因素，即分层因素。【Options】选项击Options按钮，弹出选项对话框。（图3）图3 寿命表选项对话框 Life Table(s)：输出寿命表，系统默认。Plot：统计图。Survival：累积生存函数曲线。Hazard：累积风险函数散点图。One minus survival：生存函数被1减后的曲线。Log survival：对数累积生存函数曲线。Density：密度函数散点图。Compare Levels of First Factor：对第1层因素不同水平的比较，即主对话框（图1）中的factor框中所选入的因素。 None：不做比较。系统默认。 Ov

6、erall：整体比较。 Pairwise：两两比较。 寿命表各个指标的意义Intrvl Start Time：生存时间的组段下限。Number Entrng this Intrvl：进入该组段的观察例数。Number Wdrawn During Intrval：该组段的删失例数。Number Exposed to Risk：暴露于危险因素的例数，即有效观察例数（校正观察人数）。Number of Termnl Events：终结事件的例数，即死亡例数。Propn Terminating：终结事件比例，即死亡比例。Propn Surviving：生存比例。Cumul Propn Surv at

7、End：至本组段上限的累积生存率。Probability Densty：概率密度。所有个体在时点t后单位时间内死亡概率的估计值。Hazard Rate：风险率。活过时点t的个体，在时点t后单位时间内死亡概率的估计值。SE of Cumul Surviving：累积生存率的标准误。SE of Probability Densty：概率密度的标准误。SE of Hazard Rate：风险率的标准误。 Kaplan-Meier过程采用乘积极限法（Product-limit estimates）来估计生存率，同时还可以对一个因素进行检验。适用于以个体为单位来收的未分组生存资料，是最基本的一种生存分析

8、方法。Kaplan-Meier法用于：1、估计某研究因素不同水平的中位生存时间。2、比较该研究因素不同水平的生存时间有无差异。3、控制一分层因素后对研究因素不同水平的生存时间比较（此时将按分层因素的不同水平对研究因素对生存时间的影响分别进行分析）。操作过程：1.Analyze=Survival =Kaplan-Meier2. Time框：选入 time3. Status框：选入status；击define events钮,在single value框右边的空格中输入14. Factor框：选入group；5. Compare factors列表框： Test Statistics:选择Log

9、rank、Breslow、Tarone-Ware Linear trend for factor levels:选Pooled over strata或Pairwise over strata6. Option列表框： Statistics: 选Survival table(s)、Mean and median、Survival Plots: 选Survival单击OK钮三、界面说明图1 Kaplan-Meier法主对话框【Time】框选入生存时间变量。【Status】框选入生存状态变量。【Factor】框选入分组变量。【Strata】框选入分层变量。【Lables Cases】框给个体标记。

10、【Compare Factor】选项击Compare Factor按钮，弹出选项对话框。（图2）图2 分组因素水平间比较对话框 Test Statistics:检验统计量。Log rank：检验生存分布是否相同，各时间点权重一样。Breslow：检验生存分布是否相同，以各时间点的观察例数为权重。Tarone-Ware：检验生存分布是否相同，以各时间点的观察例数的平方根为权重。Linear trend for factor levels：分组因素水平间的线性趋势检验。 Pooled over strata：水平间的整体比较。系统默认。 For each stratum：按分层变量，对每一层进行分

11、组因素各水平间的整体比较。 Pairwise over strata：分组因素各水平间的两两比较。 Pairwise for each stratum：按分层变量，对每一层进行分组因素各水平间的两两比较。【Save】 选项击Save按钮，弹出Save New Variables（保存新变量）对话框（图7）。图3 保存新变量对话框Survival：累积生存率估计。Standard error of survival：累积生存率估计的标准误。Hazard：累积风险函数估计。Cumulative events：终结事件的累积频数。在各水平内，按生存时间和生存状态排序。【Options】选项击Options按钮,弹出选项对话框（图8）。图4 K-M法选项对话框Statistics：统计量。Survival table(s)：生存分析表。Mean and median survival：平均生存时间和中位生存时间及其标准误和可信区间。Quartiles：生存时间的第25百分位数、中位生存时间、第75百分位数。Plot：统计图。Survival：累积生存函数曲线。One minus survival：生存函数被1减后的曲线。Hazard：累积风险函数散点图。Log survival：对数累积生存函数曲线。