生存分析SPSS单因素和多因素对生存率的可能分析学习教案.pptx

《生存分析SPSS单因素和多因素对生存率的可能分析学习教案.pptx》由会员分享，可在线阅读，更多相关《生存分析SPSS单因素和多因素对生存率的可能分析学习教案.pptx（31页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、会计学1生存分析生存分析(fnx)SPSS单因素和多因素对单因素和多因素对生存率的可能分析生存率的可能分析(fnx)第一页，共31页。1.何为生存分析(fnx)？生存分析(fnx)（survival analysis）是将事件的结果（终点事件）和出现结果经历的时间结合起来分析(fnx)的一种统计分析(fnx)方法。2.生存分析(fnx)的目的：描述生存过程：估计不同时间的总体生存率，计算中位生存期，绘制生存函数曲线。统计方法包括Kaplan-Meier（K-M）法、寿命表法。比较：比较不同处理组的生存率，如比较不同疗法治疗脑瘤的生存率，以了解哪种治疗方案较优。统计方法log-rank检验等。影

2、响因素分析(fnx)：研究某个或某些因素对生存率或生存时间的影响作用。如为改善脑瘤病人的预后，应了解影响病人预后的主要因素，包括病人的年龄、性别、病程、肿瘤分期、治疗方案等。统计方法cox比例风险回归模型等。预测：建立cox回归预测模型。生存生存(shngcn)分析的理论分析的理论复习复习第1页/共31页第二页，共31页。生存生存(shngcn)(shngcn)分析分析(Survival Analysis)(Survival Analysis)菜单菜单第2页/共31页第三页，共31页。寿命寿命(shumng)(shumng)表（表（Life TablesLife Tables）过）过程程Lif

3、e tables 过程用于（小样本和大样本资料）：估计某生存时间的生存率，以及中位生存时间。绘制(huzh)各种曲线：如生存函数、风险函数曲线等。对某一研究因素不同水平的生存时间分布的比较。控制另一个因素后对研究因素不同水平的生存时间分布的比较。对多组生存时间分布进行两两比较。（比较总体生存时间分布采用wilcoxon检验）第3页/共31页第四页，共31页。Company Logo实例实例(shl)分析分析例1：为了比较(bjio)不同手术方法治疗肾上腺肿瘤的疗效，某研究者随机将43例病人分成两组，甲组23例、乙组20例的生存时间（月）如下所示：其中(qzhng)有“+”者是删失数据，表示病人

4、仍生存或失访，括号内为死亡人数。（1）计算甲、乙两法术后10月的生存率和标准误。（2）估计两组的中位生存期。（3）绘制各组生存函数曲线。（4）比较两组的总体生存时间分布有无差别。第4页/共31页第五页，共31页。Company Logo一、建立数据文件（一、建立数据文件（data-01.sav）定义定义5个变量：个变量：生存时间变量：生存时间变量：t，值标签，值标签“生存时间（月）生存时间（月）”生存状态变量生存状态变量：status，取值，取值“1=死亡死亡(swng)，0=删失或存活删失或存活”频数变量：频数变量：freq，值标签，值标签“人数人数”分组变量：分组变量：group，取值，取

5、值“1=甲组，甲组，2=乙组乙组”生存时间序号变量生存时间序号变量(可无）：可无）：i第5页/共31页第六页，共31页。二、操作过程二、操作过程主菜单：分析主菜单：分析Analyze生存生存Survival寿命表寿命表Life tables对话框参数设置：对话框参数设置：时间时间time框：选入框：选入“t”。显示时间间隔显示时间间隔Display time intervals框：步长框：步长by前面填入最大生存时间的上限（必须包括生存前面填入最大生存时间的上限（必须包括生存时间最大值），步长时间最大值），步长by后面填入生存时间的组距。后面填入生存时间的组距。本例上限填本例上限填“60”，组

6、距填，组距填“1”。状态状态status框：选入框：选入“status”，击，击define events 钮，在钮，在single value框填入框填入“1”因子因子factor框：选入框：选入“group”，定义最小值，定义最小值“1”，最大值，最大值“2”。单击选项单击选项option按钮，弹出对话框：按钮，弹出对话框：1）寿命表，系统默认。寿命表，系统默认。2）图：）图：生存函数生存函数(hnsh)3）比较第一个因子的水平：）比较第一个因子的水平：整体整体比较比较第6页/共31页第七页，共31页。三、主要输出三、主要输出(shch)结果结果10月生存率的估计：月生存率的估计：甲法甲法

7、 48%，标准误，标准误 0.1 乙法乙法 30%，标准误，标准误 0.1 两组的中位生存期估计：两组的中位生存期估计：第7页/共31页第八页，共31页。Company Logo3.绘制生存绘制生存(shngcn)曲线：曲线：第8页/共31页第九页，共31页。Company Logo4.两组生存时间分布两组生存时间分布(fnb)的比较：的比较：第9页/共31页第十页，共31页。Kaplan-Meier 过程过程(guchng)Kaplan-Meier过程用于(尤其小样本资料）：估计各生存时间的生存率以及中位生存时间。绘制各种曲线：如生存函数、风险函数曲线等。比较某研究因素不同水平的生存时间有无

8、差异。控制(kngzh)某个分层因素后对研究因素不同水平的生存时间分布进行比较。对多组生存时间分布进行两两比较。（各总体分布比较采用Log-rank等非参数方法）第10页/共31页第十一页，共31页。Company Logo实例实例(shl)分析分析例2：（数据同例1）为了比较不同手术方法治疗肾上腺肿瘤的疗效，某研究者随机(su j)将43例病人分成两组，甲组23例、乙组20例的生存时间（月）如下所示：其中有“+”者是删失数据，表示病人仍生存或失访，括号内为死亡人数。（1）计算甲、乙两法各生存时间的生存率和标准误。（2）估计(gj)两组的中位生存期。（3）绘制各组生存函数曲线。（4）比较两组的

9、总体生存时间分布有无差别。第11页/共31页第十二页，共31页。一、建立数据文件（同前）一、建立数据文件（同前）二、操作过程二、操作过程主菜单主菜单(ci dn)：分析：分析Analyze生存生存SurvivalKaplan-Meier对话框参数设置：对话框参数设置：时间时间time框：选入框：选入“t”。状态状态status框：选入框：选入“status”，击，击define events 钮，在钮，在single value框填入框填入“1”。因子因子factor框：选入框：选入“group”。单击选项单击选项option按钮，弹出对话框：按钮，弹出对话框：1）统计量：）统计量：生存分析表

10、，系统默认。生存分析表，系统默认。均值和中位生存时间，系统默均值和中位生存时间，系统默认。认。2）图：）图：生存函数生存函数5.单击比较因子单击比较因子Compare Factor按钮，弹出对话按钮，弹出对话框：框：1）检验统计量）检验统计量Test Statistics：都用于检验都用于检验时间分布是否相同。时间分布是否相同。对数秩对数秩Log-rank：各时间点的权重一：各时间点的权重一样。样。Breslow：按各时间点的观察例数赋：按各时间点的观察例数赋权。权。Tarone-Ware：按各时间点观察例数：按各时间点观察例数的平方根赋权。的平方根赋权。第12页/共31页第十三页，共31页。

11、二、操作过程二、操作过程2）水平水平(shupng)间的两两比较。间的两两比较。6.单击Save按钮，弹出保存(bocn)新变量Save new variables 对话框：第13页/共31页第十四页，共31页。三、主要输出结果三、主要输出结果(ji gu)生存表：生存表：略略 两组的中位生存期估计：两组的中位生存期估计：第14页/共31页第十五页，共31页。3.绘制生存绘制生存(shngcn)曲线：曲线：第15页/共31页第十六页，共31页。Company Logo4.两组生存时间两组生存时间(shjin)分布的比较：分布的比较：第16页/共31页第十七页，共31页。Cox回归回归(hugu

12、)过程过程Cox回归过程用于：1.多个(du)因素对生存时间的影响作用分析和比较 2.生存（或死亡）风险预测第17页/共31页第十八页，共31页。实例实例(shl)分析分析n n例例3 3：为探讨某恶性肿瘤的预后，某研究者收集了：为探讨某恶性肿瘤的预后，某研究者收集了6363例患者的例患者的生存时间、生存结局及影响因素。影响因素包括病人年龄、生存时间、生存结局及影响因素。影响因素包括病人年龄、性别、组织学类型、治疗方式、淋巴结转移、肿瘤浸润程度，性别、组织学类型、治疗方式、淋巴结转移、肿瘤浸润程度，生存时间以月计算。变量的赋值和所收集的资料分别生存时间以月计算。变量的赋值和所收集的资料分别(f

13、nbi)(fnbi)见表见表17-817-8和表和表17-917-9。试用。试用CoxCox回归模型进行分析。回归模型进行分析。第18页/共31页第十九页，共31页。NoX1X2X3X4X5X6tYNoX1X2X3X4X5X6tY15400110520336200010120025701000510344011101 40135800011351355010010 26144311110103036331100012005480100071375711100120064001000600384810010120074401000580392800010 31836000112914054101

14、10120193911101700413501011 71104201001670424700010 18111420100066043491011012001242101108704443010001200135111100850454811000 151145501001820464400010 411549111017604760111001200165211101740484000010 161174811100630493201001 2411854101111010504400011 1911938010001000514810010120020401110166152720101

15、0 241213800010930534200010 212219000102415463101101200236710110930555501100 121243700110900563900010 512543100101515744000101200264900010315842111001200275011111870597400011 712853111001200606101010 401293211100120061451011010803046010011200623801000 2413143101101200636200010 1613244101101200表表17-9

16、63名某恶性肿瘤名某恶性肿瘤(xng zhng li)患者的生存时间（月）及影响因素患者的生存时间（月）及影响因素第19页/共31页第二十页，共31页。一、建立数据文件（一、建立数据文件（data-03.sav）二、操作过程二、操作过程主菜单主菜单(ci dn)：分析：分析Analyze生存生存SurvivalCox RegressionCox回归回归(hugu)主对话框主对话框第20页/共31页第二十一页，共31页。二、操作过程二、操作过程主对话框参数设置：主对话框参数设置：主对话框参数设置：主对话框参数设置：时间时间时间时间timetime框：选入框：选入框：选入框：选入“t”“t”。状态

17、状态状态状态statusstatus框：选入框：选入框：选入框：选入“y”“y”，击，击，击，击define events define events 钮，在钮，在钮，在钮，在single valuesingle value框框框框填入填入填入填入“1”“1”协变量协变量协变量协变量CovariatesCovariates框：选入框：选入框：选入框：选入“x1x6”“x1x6”。方法方法方法方法MethodMethod框：选择框：选择框：选择框：选择(xu(xu nz)nz)自变量进入自变量进入自变量进入自变量进入CoxCox模型的方法，模型的方法，模型的方法，模型的方法，SPSSSPSS提供

18、以下提供以下提供以下提供以下7 7种方法。种方法。种方法。种方法。第21页/共31页第二十二页，共31页。n n1.分析(fnx)例数描述三、主要输出三、主要输出(shch)结果结果第22页/共31页第二十三页，共31页。n n2-1.模型检验(jinyn)（全变量模型）结果结果(ji gu)提示：提示：（1）对模型总体检验有显著意义（）对模型总体检验有显著意义（P=0.003），即至少），即至少有有1个自变量的总体回归系数不为个自变量的总体回归系数不为0。第23页/共31页第二十四页，共31页。n n2-2.模型(mxng)检验（逐步回归模型(mxng)，Method=向前法LR，自变量进入

19、P0.05，剔除P0.10）（2）采用）采用(ciyng)逐步回归法进行逐步回归法进行Cox模型分析的结果提示：模型拟合自变量进入和模型分析的结果提示：模型拟合自变量进入和剔除的检验水准分别为剔除的检验水准分别为0.05和和0.1时，筛选后的最佳模型包含两个协变量，时，筛选后的最佳模型包含两个协变量，即即X4（治疗方式）和（治疗方式）和X5（淋巴结是否转移），该拟合模型总体检验提示（淋巴结是否转移），该拟合模型总体检验提示具有统计学意义（整体卡方具有统计学意义（整体卡方=17.594，P0.10）B：偏回归系数，SE：偏回归系数的标准误Wald：用于检验(jinyn)总体偏回归系数与0有无显著

20、差异，v=1时，W=（B/SE）2。Exp(B)：相对危险度估计（RR值）（3）X4（治疗方式）对生存时间有影响（治疗方式）对生存时间有影响(yngxing)，采用新疗法病人的死亡风险降至传统疗法，采用新疗法病人的死亡风险降至传统疗法的的17.2%（RR的的95%CI为为0.0590.503）。）。（4）X5（淋巴结是否转移）对生存时间也有影响（淋巴结是否转移）对生存时间也有影响(yngxing)，有淋巴结转移病人的死亡风险为，有淋巴结转移病人的死亡风险为无淋巴结转移病人的无淋巴结转移病人的2.538倍（倍（RR的的95%CI为为1.0626.066）。）。第25页/共31页第二十六页，共31

21、页。n n4.自/协变量的均值(jn zh)（选入2个变量的模型）第26页/共31页第二十七页，共31页。o5.生存曲线（选入2个变量(binling)的模型）第27页/共31页第二十八页，共31页。Company Logop 练习练习1：某临床试验对：某临床试验对20名第名第或第或第期黑色素瘤的患者进行期黑色素瘤的患者进行随访研究，截至研究期结束，记录的生存随访研究，截至研究期结束，记录的生存(shngcn)资料见表资料见表1。应用。应用SPSS软件建立数据文件，并计算软件建立数据文件，并计算100周生存周生存(shngcn)率。率。p表表1 20名第名第或或期黑色素瘤患者的治疗后生存期黑色

22、素瘤患者的治疗后生存(shngcn)时间（周）时间（周）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+练习题练习题第28页/共31页第二十九页，共31页。n n练习练习练习练习2 2：某研究者随访收集：某研究者随访收集：某研究者随访收集：某研究者随访收集了某地男性心绞痛患者了某地男性心绞痛患者了某地男性心绞痛患者了某地男性心绞痛患者24182418例，其生存时间数据例，其生存时间数据例，其生存时间数据例，其生存时间数据整理后见表整理后见表整理后见表

23、整理后见表2 2。n n（1 1）建立）建立）建立）建立SPSSSPSS数据文件。数据文件。数据文件。数据文件。n n（2 2）估计）估计）估计）估计(gj)(gj)该地男性该地男性该地男性该地男性心绞痛患者的心绞痛患者的心绞痛患者的心绞痛患者的2 2年生存率和年生存率和年生存率和年生存率和标准误，以及中位生存期。标准误，以及中位生存期。标准误，以及中位生存期。标准误，以及中位生存期。第29页/共31页第三十页，共31页。n n练习3：某研究者分别用免疫疗法、药物+免疫疗法治疗黑色素瘤患者，经随访得到(d do)各患者的生存时间（月）见数据文件ex-03.sav。n n（1）试用Kaplan-Meier过程计算其生存率和标准误。n n（2）对两组的生存率进行log-rank检验。n n（3）绘制生存曲线，计算中位生存期。第30页/共31页第三十一页，共31页。