19、生存分析SPSS.课件.ppt

1、生存分析SPSS过程（SPSS of Survival Analysis） 1. 1. 何为生存分析？何为生存分析？生存分析生存分析（survival analysis）是将是将事件的结果事件的结果（终点（终点事件）和出现结果经历的事件）和出现结果经历的时间时间结合起来分析的一种统计分析方法结合起来分析的一种统计分析方法。2. 2. 生存分析的目的：生存分析的目的：(1)(1)描述生存过程描述生存过程：估计不同时间的总体生存率，计算中位生存期，估计不同时间的总体生存率，计算中位生存期，绘制生存函数曲线。统计方法包括绘制生存函数曲线。统计方法包括Kaplan-MeierKaplan-Meier（

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

3、、性别、病程、肿瘤分期、治疗方案等。统计方法统计方法coxcox比例风险回归模型等。比例风险回归模型等。(4)(4)预测：预测：建立建立coxcox回归预测模型。回归预测模型。生存分析的理论复习生存分析的理论复习3、死亡概率、生存概率 死亡概率(mortality probability)：记为q，是指死于某时段内的可能性大小。 生存概率(survival probability)：记为p，表示某单位时段开始时存活的个体到该时段结束时仍存活的机会大小。 生存率（survival rate):S(ti)是指病人经历ti个单位时间之后仍存活的概率。截尾例数年初人口数校正人口数某年年初观察例数某年内

4、死亡数2/1q某年年初人口数某年活满一年人数qp14、生存率须分时段校正如果有截尾数据分母必观察总例数时刻仍存活的例数iiittTPtS)()(5、生存率估计的概率乘法原理： 假定病人在各个时段生存的事件独立，生存概率为P，则应用概率乘法得生存率估计的公式为：iiippptTPtS21)()(生存分析生存分析( (Survival Analysis) )菜单菜单寿命表寿命表（Life Tables）过程过程Life tables 过程用于（小样本和大样本资料）：1. 估计某生存时间的生存率，以及中位生存时间。2. 绘制各种曲线：如生存函数、风险函数曲线等。3. 对某一研究因素不同水平的生存时间

5、分布的比较。4. 控制另一个因素后对研究因素不同水平的生存时间分布的比较。5. 对多组生存时间分布进行两两比较。（比较总体生存时间分布采用wilcoxon检验）序号i确诊后年数ti期内死亡数di期内删失数ci期初病例数ni期初有效例数ni死亡概率qi生存概率pi生存率S(ti) 生存率标准误SE(S(ti) (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)109003743740.24060.75940.75940.0221217602842840.26760.73240.55620.025732 5102082080.24520.75480.41980.025543 25121

6、571510.16560.83440.35030.024854 205120117.50.17020.82980.29370.023965 799590.50.07730.92270.26820.023576 497974.50.05370.94630.25380.023387 136664.50.01550.98450.24990.023398 356259.50.05040.94960.23730.023210910255451.50.03880.96120.22810.0232表21-3 寿命表法估计生存率计算表例19-1 374名食管癌患者随访资料，取时间均为期1年，整理结果见下表，计

7、算生存率？Company Logo一、建立数据文件（一、建立数据文件（data-01.sav）定义定义3个变量：个变量：生存时间变量：t，值标签“生存时间（年）”生存状态变量 ：status，取值“1=死亡，0=删失或存活”频数变量：freq，值标签“人数”二、操作过程二、操作过程Dataweight“f”主菜单：分析主菜单：分析Analyze生存生存Survival寿命表寿命表Life tables对话框参数设置：对话框参数设置：1.时间time框：选入 “t”。2.显示时间间隔Display time intervals框：步长by前面填入最大生存时间的上限（必须包括生存时间最大值），步长

8、by后面填入生存时间的组距。本例上限填“10年”，组距填“1”。3.状态status框：选入“status”，击define events 钮，在single value框填入“1”4.因子factor框：“无”。5.单击选项option按钮，弹出对话框： 1）寿命表，系统默认。 2）图： 生存函数 三、主要输出结果三、主要输出结果生存率生存概率Company Logo绘制生存曲线：绘制生存曲线：Kaplan-Meier 过程过程Kaplan-Meier过程用于(尤其小样本资料）：1. 估计各生存时间的生存率以及中位生存时间。2. 绘制各种曲线：如生存函数、风险函数曲线等。3. 比较某研究因素

9、不同水平的生存时间有无差异。4. 控制某个分层因素后对研究因素不同水平的生存时间分布进行比较。5. 对多组生存时间分布进行两两比较。（各总体分布比较采用Log-rank等非参数方法）2、Kaplan-meier法例21-2 14例膀胱肿瘤小于3.0cm患者和16例膀胱肿瘤大于或等于3.0cm患者的生存时间（月）如下，试估计两组生存率。 0.10 ） 4. 自/协变量的均值（选入4个变量的模型） 5. 生存曲线（选入2个变量的模型）协变量均值处的生存函数Company Logop 练习练习1：某临床试验对20名第或第期黑色素瘤的患者进行随访研究，截至研究期结束，记录的生存资料见表1。应用SPSS

10、软件建立数据文件，并计算100周生存率。表1 20名第或期黑色素瘤患者的治疗后生存时间（周）12.815.624.0+26.429.230.8+39.242.058.4+72.0+77.282.487.2+94.4+97.2+106.0+114.8+117.2+140.0+168.0+练习题练习题 练习练习2：某研究者随访收集了某地男性心绞痛患者2418例，其生存时间数据整理后见表2。（1）建立SPSS数据文件。（2）估计该地男性心绞痛患者的2年生存率和标准误，以及中位生存期。 练习3： 某研究者分别用免疫疗法、药物+免疫疗法治疗黑色素瘤患者，经随访得到各患者的生存时间（月）见数据文件ex-03.sav。（1）试用Kaplan-Meier过程计算其生存率和标准误。（2）对两组的生存率进行log-rank检验。（3）绘制生存曲线，计算中位生存期。