生存時(shí)間資料分析方法課件

第一節(jié)生存時(shí)間資料的特點(diǎn)第一節(jié)生存時(shí)間資料的特點(diǎn)1一、基本概念終點(diǎn)事件與起始事件終點(diǎn)事件又稱死亡事件、失效事件，是指研究者所關(guān)心的特定結(jié)局；起始事件是反映研究對象生存過程的起始特征的事件。一、基本概念終點(diǎn)事件與起始事件2生存時(shí)間(survivaltime)又稱失效時(shí)間，為終點(diǎn)事件與起始事件之間的時(shí)間間隔，常用符號(hào)t表示。刪失值在隨訪研究中，由于某種原因未能明確地觀察到隨訪對象發(fā)生終點(diǎn)事件的時(shí)間，無法得知隨訪對象的確切生存時(shí)間，這種現(xiàn)象稱為刪失，亦稱截尾或終檢，其有關(guān)生存時(shí)間的數(shù)據(jù)稱為刪失值、截尾數(shù)據(jù)(censoreddata)或終檢數(shù)據(jù)。生存時(shí)間(survivaltime)3二、資料特點(diǎn)效應(yīng)變量有2個(gè)生存時(shí)間結(jié)局存在截尾數(shù)據(jù)分布類型復(fù)雜二、資料特點(diǎn)效應(yīng)變量有2個(gè)4三、統(tǒng)計(jì)描述三、統(tǒng)計(jì)描述5

1.死亡率、死亡概率、生存概率

死亡率(mortalityrate,deathrate)表示在某單位時(shí)間里的死亡頻率或強(qiáng)度，記為m，年死亡率的計(jì)算公式為：m=1000‰某年內(nèi)死亡數(shù)某年平均人口數(shù)

1.死亡率、死亡概率、生存概率

死亡率(mortality6死亡概率(mortalityprobability)是指死于某時(shí)段內(nèi)的可能性大小，記為q。年死亡概率的估計(jì)值為；

q=某年內(nèi)死亡數(shù)某年年初人口數(shù)死亡概率(mortalityprobability)某年內(nèi)7生存概率(survivalprobability)與死亡概率相對立，記為p。年生存概率表示年初人口往后活滿一年的機(jī)會(huì)大小，其估計(jì)值為；

p=1-q=某年活滿1年人數(shù)某年年初人口數(shù)生存概率(survivalprobability)某年活滿82.生存率、生存曲線、半數(shù)生存期生存率(survivalrate)是指觀察對象（患者）能存活到某一時(shí)點(diǎn)tk的概率，記為S（tk

）。其實(shí)質(zhì)是累積生存概率。常用的是一年生存率、五年生存率和十年生存率等。

n年生存率=100%活滿n年的例數(shù)觀察例數(shù)2.生存率、生存曲線、半數(shù)生存期生存率(survivalr9生存率估計(jì)的概率乘法原理；

S（tk

）=P（ttk

）=P1P2…Pk生存率估計(jì)的概率乘法原理；10生存曲線(survivalcurve)是指將各個(gè)時(shí)點(diǎn)tk的生存率S（tk

）連接在一起的曲線圖(縱軸：生存率，橫軸：生存時(shí)間)。階梯型折線型生存曲線(survivalcurve)11半數(shù)生存期(mediansurvivaltime)又稱中數(shù)生存期、中位生存時(shí)間，表示剛好(有且只有)50％的個(gè)體可活這么長時(shí)間。它反映生存期的平均水平。半數(shù)生存期(mediansurvivaltime)12第二節(jié)小樣本生存率的Kaplan-Meier估計(jì)直接用概率乘法原理估計(jì)生存率，又稱乘積極限法。第二節(jié)小樣本生存率的Kaplan-Meier估計(jì)直接用概率13SPSS演示SPSS演示14小樣本生存率的Kaplan-Meier估計(jì),例12-2

ViewVariable:ViewData:小樣本生存率的Kaplan-Meier估計(jì),例12-2Vi15AnalyzeSurvivalKaplan-Meier…Time框:tStatus框:cDefineEvent…SingleValue:1ContinueOptions…SurvivaltablesMeanandmediansurvivalSurvivalContinueOKAnalyzeSurvivalKaplan-Mei16生存時(shí)間資料分析方法課件17生存時(shí)間資料分析方法課件18第三節(jié)大樣本生存率的壽命表法估計(jì)第三節(jié)大樣本生存率的壽命表法估計(jì)19SPSS演示SPSS演示20大樣本生存率的壽命表法估計(jì),例12-3

ViewVariable:ViewData:大樣本生存率的壽命表法估計(jì),例12-3ViewData:21DataWeighCases…WeighCasesby:FrequencyVariable框：fOKAnalyzeSurvivalLifeTables…Time框:tDisplayTimeIntervals:10by1Status框:cDefineEvent…SingleValue:1ContinueOptions…LifetablesSurvivalContinueOKDataWeighCases…Options…22生存時(shí)間資料分析方法課件23生存時(shí)間資料分析方法課件24第四節(jié)生存曲線比較的假設(shè)檢驗(yàn)第四節(jié)生存曲線比較的假設(shè)檢驗(yàn)25Logrank檢驗(yàn)的基本思想Logrank檢驗(yàn)的基本思想26SPSS演示SPSS演示27生存曲線比較的Logrank檢驗(yàn),例12-4

ViewVariable:ViewData:生存曲線比較的Logrank檢驗(yàn),例12-4ViewD28AnalyzeSurvivalKaplan-Meier…Time框:tStatus框:cDefineEvent…SingleValue:1ContinueFactor框:gCompareFactor…LogrankContinueOptions…SurvivaltablesMeanandmediansurvivalSurvivalContinueOKAnalyzeSurvivalKaplan-Mei29生存時(shí)間資料分析方法課件30生存時(shí)間資料分析方法課件31第五節(jié)Cox比例風(fēng)險(xiǎn)模型第五節(jié)Cox比例風(fēng)險(xiǎn)模型32Cox回歸是生存分析的重要方法之一，其優(yōu)點(diǎn)是適用范圍很廣以及便于做多因素分析。Cox回歸假定病人的風(fēng)險(xiǎn)函數(shù)如下式所示：h(t)=h0(t)exp(b1X1+b2X2++bpXp)生存函數(shù)如下式所示：S(t)=exp[-h0(t)exp(b1X1+b2X2++bpXp)dt]Cox回歸是生存分析的重要方法之一，其優(yōu)點(diǎn)是適用范圍很廣以及33SPSS演示例1：為了研究惡性膽道梗阻住院患者生存時(shí)間的影響因素，從三所醫(yī)院獲得1995-2000年95例住院患者的數(shù)據(jù)，試用Cox回歸模型做生存分析。SPSS演示例1：為了研究惡性膽道梗阻住院患者生存時(shí)間的影響34

ViewVariable:ViewVariable:35ViewData:ViewData:36AnalyzeSurvivalCoxRegression…Time:timeStatus:status(1)Covariates:groupgenderagecostdurationcomplicaobstructsizestagehbstbaltastalprgtwbcMethod:Enter/ForwardLRPlots:SurvivalOptions…CIforexp(B):95%ContinueOKAnalyzeSurvivalCoxRegres37生存時(shí)間資料分析方法課件38生存時(shí)間資料分析方法課件39生存時(shí)間資料分析方法課件40生存時(shí)間資料分析方法課件41例2：某醫(yī)師在研究白血病患者的生存率時(shí)收集了33名患者的資料，按Ag陰陽性把33個(gè)數(shù)據(jù)文件分為兩組，并分別記錄了兩組患者的生存月數(shù)，是否死亡及白細(xì)胞數(shù)，試做Cox回歸分析。例2：某醫(yī)師在研究白血病患者的生存率時(shí)收集了33名患者的資料42

ViewVariable:ViewVariable:43ViewData:ViewData:44AnalyzeSurvivalCoxRegression…Time:生存時(shí)間Status:是否死亡(1)Covariates:Ag陰陽性白細(xì)胞數(shù)Method:EnterCategorical:CategoricalCovariates:Ag陰陽性

Plots:SurvivalHazardOptions…CIforexp(B):95%ContinueOKAnalyzeSurvivalCoxRegres45生存時(shí)間資料分析方法課件46生存時(shí)間資料分析方法課件47生存時(shí)間資料分析方法課件48生存時(shí)間資料分析方法課件49第一節(jié)生存時(shí)間資料的特點(diǎn)第一節(jié)生存時(shí)間資料的特點(diǎn)50一、基本概念終點(diǎn)事件與起始事件終點(diǎn)事件又稱死亡事件、失效事件，是指研究者所關(guān)心的特定結(jié)局；起始事件是反映研究對象生存過程的起始特征的事件。一、基本概念終點(diǎn)事件與起始事件51生存時(shí)間(survivaltime)又稱失效時(shí)間，為終點(diǎn)事件與起始事件之間的時(shí)間間隔，常用符號(hào)t表示。刪失值在隨訪研究中，由于某種原因未能明確地觀察到隨訪對象發(fā)生終點(diǎn)事件的時(shí)間，無法得知隨訪對象的確切生存時(shí)間，這種現(xiàn)象稱為刪失，亦稱截尾或終檢，其有關(guān)生存時(shí)間的數(shù)據(jù)稱為刪失值、截尾數(shù)據(jù)(censoreddata)或終檢數(shù)據(jù)。生存時(shí)間(survivaltime)52二、資料特點(diǎn)效應(yīng)變量有2個(gè)生存時(shí)間結(jié)局存在截尾數(shù)據(jù)分布類型復(fù)雜二、資料特點(diǎn)效應(yīng)變量有2個(gè)53三、統(tǒng)計(jì)描述三、統(tǒng)計(jì)描述54

1.死亡率、死亡概率、生存概率

死亡率(mortalityrate,deathrate)表示在某單位時(shí)間里的死亡頻率或強(qiáng)度，記為m，年死亡率的計(jì)算公式為：m=1000‰某年內(nèi)死亡數(shù)某年平均人口數(shù)

1.死亡率、死亡概率、生存概率

死亡率(mortality55死亡概率(mortalityprobability)是指死于某時(shí)段內(nèi)的可能性大小，記為q。年死亡概率的估計(jì)值為；

q=某年內(nèi)死亡數(shù)某年年初人口數(shù)死亡概率(mortalityprobability)某年內(nèi)56生存概率(survivalprobability)與死亡概率相對立，記為p。年生存概率表示年初人口往后活滿一年的機(jī)會(huì)大小，其估計(jì)值為；

p=1-q=某年活滿1年人數(shù)某年年初人口數(shù)生存概率(survivalprobability)某年活滿572.生存率、生存曲線、半數(shù)生存期生存率(survivalrate)是指觀察對象（患者）能存活到某一時(shí)點(diǎn)tk的概率，記為S（tk

）。其實(shí)質(zhì)是累積生存概率。常用的是一年生存率、五年生存率和十年生存率等。

n年生存率=100%活滿n年的例數(shù)觀察例數(shù)2.生存率、生存曲線、半數(shù)生存期生存率(survivalr58生存率估計(jì)的概率乘法原理；

S（tk

）=P（ttk

）=P1P2…Pk生存率估計(jì)的概率乘法原理；59生存曲線(survivalcurve)是指將各個(gè)時(shí)點(diǎn)tk的生存率S（tk

）連接在一起的曲線圖(縱軸：生存率，橫軸：生存時(shí)間)。階梯型折線型生存曲線(survivalcurve)60半數(shù)生存期(mediansurvivaltime)又稱中數(shù)生存期、中位生存時(shí)間，表示剛好(有且只有)50％的個(gè)體可活這么長時(shí)間。它反映生存期的平均水平。半數(shù)生存期(mediansurvivaltime)61第二節(jié)小樣本生存率的Kaplan-Meier估計(jì)直接用概率乘法原理估計(jì)生存率，又稱乘積極限法。第二節(jié)小樣本生存率的Kaplan-Meier估計(jì)直接用概率62SPSS演示SPSS演示63小樣本生存率的Kaplan-Meier估計(jì),例12-2

ViewVariable:ViewData:小樣本生存率的Kaplan-Meier估計(jì),例12-2Vi64AnalyzeSurvivalKaplan-Meier…Time框:tStatus框:cDefineEvent…SingleValue:1ContinueOptions…SurvivaltablesMeanandmediansurvivalSurvivalContinueOKAnalyzeSurvivalKaplan-Mei65生存時(shí)間資料分析方法課件66生存時(shí)間資料分析方法課件67第三節(jié)大樣本生存率的壽命表法估計(jì)第三節(jié)大樣本生存率的壽命表法估計(jì)68SPSS演示SPSS演示69大樣本生存率的壽命表法估計(jì),例12-3

ViewVariable:ViewData:大樣本生存率的壽命表法估計(jì),例12-3ViewData:70DataWeighCases…WeighCasesby:FrequencyVariable框：fOKAnalyzeSurvivalLifeTables…Time框:tDisplayTimeIntervals:10by1Status框:cDefineEvent…SingleValue:1ContinueOptions…LifetablesSurvivalContinueOKDataWeighCases…Options…71生存時(shí)間資料分析方法課件72生存時(shí)間資料分析方法課件73第四節(jié)生存曲線比較的假設(shè)檢驗(yàn)第四節(jié)生存曲線比較的假設(shè)檢驗(yàn)74Logrank檢驗(yàn)的基本思想Logrank檢驗(yàn)的基本思想75SPSS演示SPSS演示76生存曲線比較的Logrank檢驗(yàn),例12-4

ViewVariable:ViewData:生存曲線比較的Logrank檢驗(yàn),例12-4ViewD77AnalyzeSurvivalKaplan-Meier…Time框:tStatus框:cDefineEvent…SingleValue:1ContinueFactor框:gCompareFactor…LogrankContinueOptions…SurvivaltablesMeanandmediansurvivalSurvivalContinueOKAnalyzeSurvivalKaplan-Mei78生存時(shí)間資料分析方法課件79生存時(shí)間資料分析方法課件80第五節(jié)Cox比例風(fēng)險(xiǎn)模型第五節(jié)Cox比例風(fēng)險(xiǎn)模型81Cox回歸是生存分析的重要方法之一，其優(yōu)點(diǎn)是適用范圍很廣以及便于做多因素分析。Cox回歸假定病人的風(fēng)險(xiǎn)函數(shù)如下式所示：h(t)=h0(t)exp(b1X1+b2X2++bpXp)生存函數(shù)如下式所示：S(t)=exp[-h0(t)exp(b1X1+b2X2++bpXp)dt]Cox回歸是生存分析的重要方法之一，其優(yōu)點(diǎn)是適用范圍很廣以及82SPSS演示例1：為了研究惡性膽道梗阻住院患者生存時(shí)間的影響因素，從三所醫(yī)院獲得1995-2000年95例住院患者的數(shù)據(jù)，試用Cox回歸模型做生存分析。SPSS演示例1：為了研究惡性膽道梗阻住院患者生存時(shí)間的影響83

ViewVariable:ViewVariable:84ViewData:ViewData:85AnalyzeSurvivalCoxRegression…Time:time