資源描述:
《同濟(jì)醫(yī)學(xué)院SPSS課件-SPSS Cox回歸》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、生存分析第10章非參數(shù)檢驗(yàn)SPSS應(yīng)用華中科技大學(xué)公共衛(wèi)生學(xué)院流行病學(xué)與衛(wèi)生統(tǒng)計(jì)學(xué)系蔣紅衛(wèi)jhwccc@21cn.com內(nèi)容基本概念基本步驟基本操作基本結(jié)果解釋生存分析(survivalanalysis)是將事件的結(jié)果(終點(diǎn)事件)和出現(xiàn)這一結(jié)果所經(jīng)歷的時間結(jié)合起來分析的一種統(tǒng)計(jì)分析方法。生存分析不同于其它多因素分析的主要區(qū)別點(diǎn)就是生存分析考慮了每個觀測出現(xiàn)某一結(jié)局的時間長短?;靖拍钌娣治龅哪康拿枋錾孢^程,估計(jì)不同時間的總體生存率,計(jì)算中位生存期,繪制生存函數(shù)曲線。Kaplan-Meier(K-M)法壽命表法比較:比較不同處理
2、組的生存率,如比較不同療法治療腦瘤的生存率,以了解哪種治療方案較優(yōu)。log-rank檢驗(yàn)影響因素分析:研究某個或某些因素對生存率或生存時間的影響作用。如,了解影響腦瘤病人預(yù)后的主要因素,包括年齡、性別、腫瘤分期、治療方案等。cox比例風(fēng)險(xiǎn)回歸模型預(yù)測cox回歸預(yù)測模型生存分析菜單壽命表法壽命表法用于(小樣本和大樣本資料)估計(jì)某生存時間的生存率,以及中位生存時間繪制各種曲線:如生存函數(shù)、風(fēng)險(xiǎn)函數(shù)曲線等對某一研究因素不同水平的生存時間分布的比較控制另一個因素后對研究因素不同水平的生存時間分布的比較對多組生存時間分布進(jìn)行兩兩比較(比較總
3、體生存時間分布采用wilcoxon檢驗(yàn))例1:為了比較不同手術(shù)方法治療腎上腺腫瘤的療效,43例病人隨機(jī)分為兩組其中有“+”者是刪失數(shù)據(jù),表示病人仍生存或失訪,括號內(nèi)為死亡人數(shù)。例1文件格式定義5個變量生存時間變量:Time,值標(biāo)簽“生存時間(月)”生存狀態(tài)變量:status,取值“0=死亡,1=刪失或存活”頻數(shù)變量(可無):f,值標(biāo)簽“人數(shù)”分組變量:group,取值“0=甲組,1=乙組”生存時間序號變量(可無):IDSPSS基本操作SPSS基本操作生存時間最大生存時間組距生存狀態(tài)定義終點(diǎn)事件分組因素分組定義生存狀態(tài)定義生存狀態(tài)定
4、義SPSS基本操作給定選項(xiàng)SPSS基本操作生成壽命表生存曲線作圖等各組比較的方式不作比較整體性比較兩兩比較Kaplan-Meier法Kaplan-Meier過程用于(尤其小樣本資料)估計(jì)各生存時間的生存率以及中位生存時間。繪制各種曲線:如生存函數(shù)、風(fēng)險(xiǎn)函數(shù)曲線等。比較某研究因素不同水平的生存時間有無差異。控制某個分層因素后對研究因素不同水平的生存時間分布進(jìn)行比較。對多組生存時間分布進(jìn)行兩兩比較(各總體分布比較采用Log-rank等非參數(shù)方法)同例1:為了比較不同手術(shù)方法治療腎上腺腫瘤的療效,43例病人隨機(jī)分為兩組其中有“+”者是刪
5、失數(shù)據(jù),表示病人仍生存或失訪,括號內(nèi)為死亡人數(shù)。例2SPSS基本操作SPSS基本操作生存時間生存狀態(tài)定義生存狀態(tài)分組因素SPSS基本操作SPSS基本操作比較方法SPSS基本操作比較方法比較方式幾點(diǎn)說明檢驗(yàn)統(tǒng)計(jì)量TestStatistics:都用于檢驗(yàn)時間分布是否相同對數(shù)秩Log-rank:各時間點(diǎn)的權(quán)重一樣。Breslow:按各時間點(diǎn)的觀察例數(shù)賦權(quán)Tarone-Ware:按各時間點(diǎn)觀察例數(shù)的平方根賦權(quán)Lineartrendforfactorlevels:線性趨勢檢驗(yàn)Pooledoverstrata:水平間整體性比較Foreachs
6、tratum:每層進(jìn)行各分組水平間的整體性比較Pairwiseoverstrata:分組各水平間的兩兩比較Pairwiseforeachstratum:每層上各水平間的兩兩比較SPSS基本操作選項(xiàng)SPSS基本操作統(tǒng)計(jì)指標(biāo)作圖Cox回歸在醫(yī)學(xué)中,一方面要看治療結(jié)局的好壞,另一方面還要看生存時間的長短。生存時間的長短不僅與治療措施有關(guān),還可能與病人的體質(zhì),年齡,病情的輕重等多種因素有關(guān)。如何找出其中哪些因素與生存時間有關(guān)、哪些無關(guān)呢?由于失訪、試驗(yàn)終止等原因造成某些時間的不完全,不能用多元線性回歸分析。1972年英國統(tǒng)計(jì)學(xué)家CoxDR
7、.提出一種比例危險(xiǎn)模型方法,能處理多個因素對生存時間影響的問題。其目的多個因素對生存時間的影響作用分析和比較生存(或死亡)風(fēng)險(xiǎn)預(yù)測Cox回歸Cox回歸屬半?yún)?shù)模型,模型基本形式:h(t,x)為具有協(xié)變量x的個體在時刻t的風(fēng)險(xiǎn)函數(shù)(風(fēng)險(xiǎn)率,瞬時死亡率),h0(t)稱基準(zhǔn)風(fēng)險(xiǎn)率,即所有協(xié)變量均為0時,個體時刻t的風(fēng)險(xiǎn)函數(shù)上式的右側(cè)分兩部分:h0(t)與時間有關(guān)的任意函數(shù),其分布和形狀無明確假定,是非參數(shù)部分.另一部分是參數(shù)部分,其參數(shù)可以通過樣本估計(jì).因?yàn)闊o須估計(jì)h0(t),故Cox回歸的模型擬和不是直接用生存時間作為因變量,而是以風(fēng)
8、險(xiǎn)函數(shù)與基礎(chǔ)風(fēng)險(xiǎn)函數(shù)的比值為因變量。Cox比例風(fēng)險(xiǎn)模型Cox回歸-多因素分析數(shù)據(jù)格式設(shè)含有p個變量x1,x2,…,xp及時間T和結(jié)局C的n個觀察對象.其數(shù)據(jù)結(jié)構(gòu)為:IDx1x2…xpTimeStatus1x11x12…x1py102x21x22…x