《醫(yī)學(xué)生存分析》

生存分析(SurvivalAnalysis)王靜安徽醫(yī)科大學(xué)流行病與衛(wèi)生統(tǒng)計(jì)學(xué)系5/17/20241安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作例1：25名惡性腫瘤病人隨機(jī)分成A（12人）、B（13人）兩組，分別采用手術(shù)及保守治療兩種方案進(jìn)行治療，之后隨訪觀察直到25名病人均死亡，記錄下他們的存活時(shí)間（單位：天）。見下表所示。5/17/20242安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/20243安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作例2：25名惡性腫瘤病人隨機(jī)分成A（12人）、B（13人）兩組，分別采用手術(shù)及保守治療兩種方案進(jìn)行治療，之后隨訪觀察10年，記錄他們的存活、死亡情況。見下表所示。5/17/20244安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/20245安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作對(duì)于例1，只比較兩組的生存時(shí)間長短即可，采用兩樣本秩和檢驗(yàn)。對(duì)于例2，若比較兩組的存活率，可采用兩樣本卡方檢驗(yàn)；但只考慮了生存效果評(píng)價(jià)的其中一個(gè)方面。而評(píng)價(jià)生存效果，不僅要看是否出現(xiàn)了某種結(jié)局（如：有效、治愈、死亡等），還要考慮出現(xiàn)這些結(jié)局所經(jīng)歷的時(shí)間長短。（舉例說明）5/17/20246安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作措施：采用生存分析方法。生存分析方法特點(diǎn)：（1）能將研究對(duì)象的隨訪結(jié)局和生存時(shí)間兩個(gè)因素同時(shí)結(jié)合起來考慮的一種統(tǒng)計(jì)方法；（2）能處理失訪等不完全數(shù)據(jù)，充分利用所獲得的信息，對(duì)生存時(shí)間的分布特征進(jìn)行描述、比較，對(duì)影響生存時(shí)間的主要因素進(jìn)行分析，達(dá)到全面評(píng)價(jià)和比較隨訪資料的目的。5/17/20247安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作生存分析方法的用途：用于隨訪研究（即：觀察結(jié)果并非在短期內(nèi)能夠確定，而需做長期隨訪觀察，如對(duì)一些慢性病或惡性腫瘤的預(yù)后及遠(yuǎn)期療效觀察等）。5/17/20248安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作隨訪研究的特點(diǎn)醫(yī)學(xué)隨訪研究：1、隊(duì)列研究——所有被觀察對(duì)象同時(shí)進(jìn)入研究；2、臨床隨訪研究（一般稱“臨床試驗(yàn)”）——被觀察對(duì)象逐個(gè)進(jìn)入研究。在大多數(shù)研究中，由于受經(jīng)費(fèi)和時(shí)間的限制，最終觀察時(shí)點(diǎn)是固定的，而不是無限制延長的，所以容易產(chǎn)生截尾值。5/17/20249安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作截尾值假設(shè)一組肺癌病人手術(shù)后，隨訪記錄他們的存活情況，死于肺癌是終點(diǎn)。當(dāng)觀察到規(guī)定的時(shí)點(diǎn)，除了死于肺癌，還有：

死于其他病、尚活著、因遷移等原因失去聯(lián)系（即失訪），后三者的觀察值都未達(dá)到終點(diǎn)，不能提供完全的信息，這種不完全數(shù)據(jù)稱為截尾值（censoring）。5/17/202410安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作

終點(diǎn)事件失訪

研究時(shí)間

觀察起點(diǎn)觀察終點(diǎn)

隊(duì)列研究中的截尾情況5/17/202411安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作

終點(diǎn)事件失訪

研究時(shí)間

觀察起點(diǎn)觀察終點(diǎn)

臨床隨訪研究中的截尾情況5/17/202412安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作基本概念5/17/202413安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作1、生存時(shí)間（survivaltime）（狹義）從發(fā)病到死亡所經(jīng)歷的時(shí)間；（廣義）從某種起始事件到某種終點(diǎn)事件所經(jīng)歷的時(shí)間，又稱失效時(shí)間(failuretime)，表示為t。5/17/202414安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作“生存”和“死亡”擴(kuò)展為具有兩類互斥結(jié)果的事件，“生存”是某事件某狀態(tài)的持續(xù)，“死亡”是某終點(diǎn)事件的發(fā)生?！吧鏁r(shí)間”擴(kuò)展為某事件某狀態(tài)所持續(xù)的時(shí)間或過程。5/17/202415安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作2、起點(diǎn)事件和終點(diǎn)事件（endpointevent）終點(diǎn)事件也稱為觀察結(jié)果（outcome），是指研究者關(guān)心的研究對(duì)象的特定結(jié)局，當(dāng)被觀察對(duì)象出現(xiàn)終點(diǎn)事件通常記為1，出現(xiàn)截尾記為0。起點(diǎn)事件是反映研究對(duì)象生存過程的起始特征的事件。生存時(shí)間的起點(diǎn)和終點(diǎn)要有嚴(yán)格、統(tǒng)一的定義，以便于計(jì)算生存時(shí)間。5/17/202416安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作3、暴露因素及混雜因素除了要研究的危險(xiǎn)因素（暴露因素）外，還要控制混雜因素的影響。必須分清暴露因素和混雜因素。如：在比較A、B兩種治療方案治療某惡性疾病的生存情況的同時(shí)，病人的腎功能狀況也影響治療后的生存情況。這項(xiàng)研究中，暴露因素為A、B兩種治療方案，混雜因素是腎功能狀況。5/17/202417安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作隨訪資料的特點(diǎn)1、生存時(shí)間的分布呈偏態(tài)。一般為正偏態(tài)，因?yàn)樯鏁r(shí)間一般為單調(diào)遞減。102名黑色素瘤患者的生存時(shí)間資料，做直方圖。5/17/202418安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202419安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202420安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作2、數(shù)據(jù)中常含有截尾值。如41+（天），表示該病人至少活了41天；雖然不知道確切時(shí)間，但此類不完全數(shù)據(jù)包含一定信息，需要生存分析這特殊方法來利用這些信息。組別生存時(shí)間（天）A273541+……B9980+102……5/17/202421安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202422安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作生存分析的主要內(nèi)容及研究方法5/17/202423安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作

主要內(nèi)容研究方法1、生存過程的描述

乘積-極限法（Kaplan-Meier）

壽命表法（LifeTables）2、生存過程的比較

對(duì)數(shù)秩檢驗(yàn)（log-rank）3、影響因素的分析

Cox比例風(fēng)險(xiǎn)模型

5/17/202424安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作生存分析對(duì)資料中應(yīng)變量的要求：5/17/202425安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作1、達(dá)到終點(diǎn)的例數(shù)所占的比例不能太少，即完全數(shù)據(jù)占大部分，截尾值所占的比例要＜10%；2、截尾原因無偏性；3、生存時(shí)間盡可能精確。5/17/202426安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作常用的研究指標(biāo)5/17/202427安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作1、生存函數(shù)（survivalfunction）又稱累積生存概率（cumulativesurvivalprobability）、生存率，表示個(gè)體生存時(shí)間T≥t的概率。5/17/202428安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作生存概率與生存率生存概率：在單位時(shí)段開始存活的個(gè)體到該時(shí)段結(jié)束時(shí)仍存活的可能性，用p表示；生存率：觀察個(gè)體活過t個(gè)單位時(shí)間的概率，用S(t)表示。如數(shù)據(jù)中無截尾值，則生存率計(jì)算公式為：S(t)=生存時(shí)間≥t的病人數(shù)/隨訪病人總數(shù)如有截尾值，則分時(shí)段計(jì)算不同單位時(shí)段的生存概率p1、p2、…pt，S(t)=p1×p2×…×pt，即累積生存概率。5/17/202429安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202430安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作可采用SPSS軟件中生存分析方法計(jì)算出生存率大小，并畫出生存函數(shù)圖（或生存率曲線）直接觀察。5/17/202431安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202432安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202433安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作2、半數(shù)生存時(shí)間(mediansurvivaltime)又稱中位生存時(shí)間，表示累積生存概率（生存率）為50%的時(shí)候所對(duì)應(yīng)的生存時(shí)間；為中位數(shù)指標(biāo)，因?yàn)樯鏁r(shí)間的分布是偏態(tài)的。在分析結(jié)果中給出中位生存時(shí)間；從生存曲線圖中也能估計(jì)出相同結(jié)果。5/17/202434安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202435安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作3、風(fēng)險(xiǎn)函數(shù)(hazardfunction)又稱危險(xiǎn)率函數(shù)等，一般用h(t)表示，表示生存時(shí)間已達(dá)t的個(gè)體在t時(shí)刻的瞬時(shí)風(fēng)險(xiǎn)率；在Cox比例風(fēng)險(xiǎn)模型中風(fēng)險(xiǎn)函數(shù)常用λ(t)表示。h(t)=死于區(qū)間（t，t+⊿t）的病人數(shù)在t時(shí)刻尚存的病人數(shù)×⊿t⊿t為極小時(shí)間段。5/17/202436安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202437安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作4、風(fēng)險(xiǎn)比（hazardratio）是指同一時(shí)點(diǎn)兩組的風(fēng)險(xiǎn)函數(shù)之比，即相對(duì)危險(xiǎn)度RR。風(fēng)險(xiǎn)比=第一組個(gè)體的h1(t)

第二組個(gè)體的h2(t)

（1）比例風(fēng)險(xiǎn)：風(fēng)險(xiǎn)比與時(shí)間無關(guān)，即任何時(shí)刻，兩組的風(fēng)險(xiǎn)比值是相等的；（2）非比例風(fēng)險(xiǎn)或時(shí)間依賴型：風(fēng)險(xiǎn)比與時(shí)間有關(guān)。5/17/202438安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作乘積—極限法（Kaplan-Meier）壽命表法（Lifetable）Analyze→Survival→Lifetable、Kaplan-Meier用途：1、做出生存表、估計(jì)生存率，并繪出兩組的生存曲線圖和風(fēng)險(xiǎn)函數(shù)圖；2、進(jìn)行兩組生存資料的比較。5/17/202439安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作

LifetableKaplan-Meier共同點(diǎn)：非參數(shù)分析方法、一般用于單因素分析適用條件大樣本資料大樣本、小樣本資料生存表人為劃分各時(shí)間段每個(gè)時(shí)間值為一段多組間的生存率比較方法（方法選擇途徑）

“Options”對(duì)話框“CompareFactors”對(duì)話框5/17/202440安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作見下例word文檔中的數(shù)據(jù)及對(duì)應(yīng)的SPSS數(shù)據(jù)文件KM1.sav，采用乘積極限法進(jìn)行分析。注意：1）SPSS數(shù)據(jù)文件格式；2）SPSS中乘積極限法的操作過程；3）乘積極限法的具體結(jié)果解釋。5/17/202441安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202442安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202443安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202444安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202445安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202446安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202447安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202448安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202449安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作注意事項(xiàng)：組間的生存率比較時(shí)，要求各組的生存曲線不能交叉，若出現(xiàn)交叉，提示可能存在混雜因素，應(yīng)采用多因素方法來校正混雜因素或分層做統(tǒng)計(jì)分析。5/17/202450安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作采用Lifetable法對(duì)某生存數(shù)據(jù)進(jìn)行分析5/17/202451安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202452安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202453安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作生存分析中的多因素分析方法5/17/202454安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202455安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202456安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作（生存分析方法中單因素與多因素分析的區(qū)別）若只比較男、女的生存情況優(yōu)劣，則可采用前面所述的單因素分析方法——乘積極限法、壽命表法。而在實(shí)際工作中，研究者關(guān)心的影響因素不止一個(gè)，如在此例中，除了性別，還有年齡等其他因素。因此，應(yīng)采用多因素分析方法。5/17/202457安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作（多元線性回歸、多元logistic回歸與生存分析中多因素回歸分析方法的區(qū)別）如果僅以生存時(shí)間為反應(yīng)變量做多元線性回歸，則生存時(shí)間并不一定服從近似正態(tài)分布，不滿足線性回歸的要求；若僅以是否達(dá)到終點(diǎn)為反應(yīng)變量做多元logistic回歸，則未能充分利用生存時(shí)間長短的信息。而且兩者均未考慮到資料中刪失數(shù)據(jù)提供的信息。5/17/202458安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作解決措施：采用生存分析中的多因素回歸。生存分析中的多因素回歸同時(shí)以生存時(shí)間、結(jié)局為反應(yīng)變量，有效利用不完全數(shù)據(jù)提供的信息，從而探索生存的影響因素。5/17/202459安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作生存分析中的多因素回歸包括：半?yún)?shù)方法：cox比例風(fēng)險(xiǎn)模型、非比例風(fēng)險(xiǎn)模型參數(shù)方法：指數(shù)分布模型、Weibull分布模型、Gompertz分布模型、對(duì)數(shù)正態(tài)分布模型、對(duì)數(shù)logistic分布模型，等。參數(shù)方法中均假設(shè)生存時(shí)間服從某分布，而半?yún)?shù)方法中不需指定生存時(shí)間的分布情況。5/17/202460安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作Cox比例風(fēng)險(xiǎn)模型5/17/202461安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作適用情況：用于分析帶有伴隨變量的生存時(shí)間資料，如：腫瘤和其它慢性病的預(yù)后分析，臨床療效評(píng)價(jià)和隊(duì)列研究的病因探索。優(yōu)點(diǎn)：適用條件寬，便于作多因素分析。5/17/202462安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作一、cox回歸模型簡介5/17/202463安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作1、模型結(jié)構(gòu)即Cox回歸，它的一般形式是：

λ(t)=[λ0(t)]?e(β1*x1+β2*x2+….+βm*xm

)

=[λ0(t)]?exp(β1*x1+β2*x2+….+βm*xm)5/17/202464安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作λ(t)表示t時(shí)刻暴露于各危險(xiǎn)因素（x1，x2，….xm）狀態(tài)下的風(fēng)險(xiǎn)函數(shù)。λ0(t)為基線風(fēng)險(xiǎn)函數(shù)，表示所有的危險(xiǎn)因素狀態(tài)都為0情況下的風(fēng)險(xiǎn)函數(shù)。λ(t)/λ0(t)=exp(β1*x1+β2*x2+….+βm*xm

)即：

ln[λ(t)/λ0(t)]=β1*x1+β2*x2+….+βm*xm將cox回歸與線性回歸、logistic回歸做一比較。5/17/202465安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作λ0(t)的形式不限，使模型的適用范圍如同非參數(shù)法那么廣；模型中含有參數(shù)β，又使得其統(tǒng)計(jì)效率接近參數(shù)模型。（但參數(shù)β就不能用傳統(tǒng)的方法進(jìn)行估計(jì)和檢驗(yàn)）∴模型中含有參數(shù)β，但基線風(fēng)險(xiǎn)函數(shù)λ0(t)不要求服從特定分布形式，具有非參數(shù)的特點(diǎn),故Cox回歸又稱為半?yún)?shù)模型。5/17/202466安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作英國生物統(tǒng)計(jì)學(xué)家D.R.Cox于1972年建立了條件死亡概率和偏似然函數(shù)方法，解決了參數(shù)β的估計(jì)和檢驗(yàn)問題。故，該統(tǒng)計(jì)分析方法稱為Cox模型。5/17/202467安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作2、參數(shù)的統(tǒng)計(jì)學(xué)意義假設(shè)只有一個(gè)X危險(xiǎn)因素，0為“不暴露”，1為“暴露”，模型形式為：λ(t)=λ0(t)?exp(β*x)X=1時(shí)的風(fēng)險(xiǎn)函數(shù)為：λ1(t)=λ0(t)?exp(β*1)X=0時(shí)的風(fēng)險(xiǎn)函數(shù)為：λ2(t)=λ0(t)?exp(β*0)λ1(t)、λ2(t)分別表示暴露于危險(xiǎn)因素的兩種不同狀態(tài)下發(fā)病的風(fēng)險(xiǎn)函數(shù)。5/17/202468安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作暴露于危險(xiǎn)因素的兩種不同狀態(tài)下發(fā)病的風(fēng)險(xiǎn)比即為相對(duì)危險(xiǎn)度RR。RR=λ1(t)/λ2(t)=[λ0(t)?exp(β*1)]/[λ0(t)?exp(β*0)]=exp(β*1)/exp(β*0)=exp（β）RR表示暴露組與非暴露組的風(fēng)險(xiǎn)函數(shù)之比，即：暴露組發(fā)病的風(fēng)險(xiǎn)是非暴露組的RR倍。回歸系數(shù)β的流行病學(xué)含義是0、1協(xié)變量X的相對(duì)危險(xiǎn)度的自然對(duì)數(shù)。5/17/202469安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作如果X為連續(xù)性變量——年齡（歲），則λ1(t)、λ2(t)分別表示k歲、k+1歲的兩種不同狀態(tài)下發(fā)病的風(fēng)險(xiǎn)函數(shù)。RR=λ1(t)/λ2(t)=[λ0(t)?exp(β*（k+1）)]/[λ0(t)?exp(β*k)]=exp(β*（k+1）

/exp(β*k)=exp（β）β：X每增加一個(gè)單位時(shí)其相對(duì)危險(xiǎn)度的自然對(duì)數(shù)值。RR表示某危險(xiǎn)因素（或協(xié)變量）改變一個(gè)測量單位時(shí)風(fēng)險(xiǎn)函數(shù)改變多少倍。5/17/202470安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)β＞0，說明相應(yīng)協(xié)變量值的增加將增大所研究事件發(fā)生的可能性；當(dāng)β＜0，相應(yīng)協(xié)變量值的增加將減少所研究事件發(fā)生的可能性；當(dāng)β=0，相應(yīng)協(xié)變量與所研究事件的發(fā)生無關(guān)。5/17/202471安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作3、模型假設(shè)滿足比例風(fēng)險(xiǎn)（proportionalhazards）假定，簡稱PH假定。即：相對(duì)危險(xiǎn)度RR或風(fēng)險(xiǎn)比保持一個(gè)恒定的比例，與時(shí)間t無關(guān)。5/17/202472安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作Cox回歸無須對(duì)λ0(t)的形式作任何限制，因?yàn)椴挥绊懜魑ｋU(xiǎn)因素相對(duì)危險(xiǎn)度的估計(jì)，而相對(duì)危險(xiǎn)度正是多因素分析時(shí)最關(guān)注的問題?！郼ox模型巧妙地將非參數(shù)[λ0(t)]部分與參數(shù)（回歸系數(shù)β）的概念結(jié)合起來，這種靈活性使得它在生存分析的應(yīng)用中備受青睞。5/17/202473安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作比值比OR可以作為相對(duì)危險(xiǎn)度RR的估計(jì)值。5/17/202474安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作比例風(fēng)險(xiǎn)假設(shè)的檢驗(yàn)：可分組（若為0、1兩組）做二次對(duì)數(shù)生存率圖（LogminusLog，LML），進(jìn)行目測判斷，若兩條曲線有交叉或間距變動(dòng)大，則提示風(fēng)險(xiǎn)比隨時(shí)間而改變，不具備cox比例風(fēng)險(xiǎn)模型的適用性；反之，則資料適宜用cox比例風(fēng)險(xiǎn)模型。5/17/202475安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202476安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作4、協(xié)變量的篩選策略與其他回歸模型類似，如果協(xié)變量個(gè)數(shù)較多，則通常采用逐步法。5/17/202477安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作二、cox回歸模型的分析步驟及主要結(jié)果5/17/202478安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作SPSS中的程序?yàn)椋篈nalyze→Survival→CoxRegression將cox回歸的對(duì)話框與線性回歸、Kaplan-Meier、Lifetable做比較。見文件cox1.sav5/17/202479安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作示例：評(píng)價(jià)A、B兩治療方案對(duì)某病的治療效果。A組（group=0）12人，B組（group=1）13人；腎功能正常者為0，不正常者為1。應(yīng)變量是治療后的生存時(shí)間和觀察結(jié)果（達(dá)到終點(diǎn)還是截尾值）；group是研究因素x1，腎功能是混雜因素x2。5/17/202480安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作注意事項(xiàng)一般情況下，截尾值設(shè)為0，完全數(shù)據(jù)（即出現(xiàn)終點(diǎn)事件）為1；但在具體資料中，沒有限制，只要能有兩個(gè)分別表示完全數(shù)據(jù)和截尾的數(shù)值，結(jié)果不受影響。但是一旦設(shè)定后，分析時(shí)，“Status”對(duì)話框中的終點(diǎn)事件值一定要和資料中吻合，否則結(jié)果錯(cuò)誤。5/17/202481安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作SPSS分析結(jié)果如下：第一步、求Cox比例風(fēng)險(xiǎn)模型5/17/202482安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202483安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作λ(t)=λ0(t)·e(1.243group+4.105kidney)=λ0(t)·exp(1.243group+4.105kidney)group、kidney的P值分別為0.038、0.000，均＜0.05，說明不同治療方案及腎功能對(duì)病人的生存有影響。5/17/202484安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作group、kidney的OR值分別為3.466、60.670，均大于1，說明接受B治療方案的病人死亡的風(fēng)險(xiǎn)是接受A治療方案者的3.466倍、腎功能不正常的病人死亡的風(fēng)險(xiǎn)是腎功能正常者的60.670倍。5/17/202485安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作第二步、由模型可求出不同狀態(tài)下的相對(duì)危險(xiǎn)度RR值5/17/202486安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作1、腎功能正常者接受B治療方案比接受A治療方案在某時(shí)刻死亡的相對(duì)危險(xiǎn)度為：RR=[λ0(t)]?e(1.243*1+4.105*0)

[λ0(t)]?e(1.243*0+4.105*0)=3.4665/17/202487安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作2、腎功能不正常者接受B治療方案比接受A治療方案在某時(shí)刻死亡的相對(duì)危險(xiǎn)度為：RR=[λ0(t)]?e(1.243*1+4.105*1)

[λ0(t)]?e(1.243*0+4.105*1)=3.4665/17/202488安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作3、腎功能不正常者接受B治療方案，比腎功能正常者接受A治療方案在某時(shí)刻死亡的相對(duì)危險(xiǎn)度為：RR=[λ0(t)]?e(1.243*1+4.105*1)

[λ0(t)]?e(1.243*0+4.105*0)=210.300可見在Cox模型的建立中沒有估計(jì)基線風(fēng)險(xiǎn)λ0(t)；在危險(xiǎn)因素的RR分析中也沒有涉及基線風(fēng)險(xiǎn)。所估計(jì)的風(fēng)險(xiǎn)比與時(shí)間無關(guān)。5/17/202489安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作第三步、畫出風(fēng)險(xiǎn)函數(shù)曲線5/17/202490安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202491安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作第四步、畫出二次對(duì)數(shù)生存率圖（LogminusLog，LML）從LML圖中，目測判斷PH假設(shè)是否成立。5/17/202492安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作5/17/202493安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作目的:探討影響胃癌患者長期生存的預(yù)后因素。方法:采用Cox比例