王靜老師最全的生存分析

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

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

終點事件失訪

研究時間

觀察起點觀察終點

隊列研究中的截尾情況2023/3/2910安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前10頁，總共108頁。

終點事件失訪

研究時間

觀察起點觀察終點

臨床隨訪研究中的截尾情況2023/3/2911安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前11頁，總共108頁?；靖拍?023/3/2912安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前12頁，總共108頁。1、生存時間（survivaltime）（狹義）從發(fā)病到死亡所經(jīng)歷的時間；（廣義）從某種起始事件到某種終點事件所經(jīng)歷的時間，又稱失效時間(failuretime)，表示為t。2023/3/2913安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前13頁，總共108頁。“生存”和“死亡”擴展為具有兩類互斥結(jié)果的事件，“生存”是某事件某狀態(tài)的持續(xù)，“死亡”是某終點事件的發(fā)生。

“生存時間”擴展為某事件某狀態(tài)所持續(xù)的時間或過程。2023/3/2914安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前14頁，總共108頁。2、起點事件和終點事件（endpointevent）終點事件也稱為觀察結(jié)果（outcome），是指研究者關(guān)心的研究對象的特定結(jié)局，當(dāng)被觀察對象出現(xiàn)終點事件通常記為1，出現(xiàn)截尾記為0。起點事件是反映研究對象生存過程的起始特征的事件。生存時間的起點和終點要有嚴(yán)格、統(tǒng)一的定義，以便于計算生存時間。2023/3/2915安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前15頁，總共108頁。3、暴露因素及混雜因素除了要研究的危險因素（暴露因素）外，還要控制混雜因素的影響。必須分清暴露因素和混雜因素。如：在比較A、B兩種治療方案治療某惡性疾病的生存情況的同時，病人的腎功能狀況也影響治療后的生存情況。這項研究中，暴露因素為A、B兩種治療方案，混雜因素是腎功能狀況。2023/3/2916安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前16頁，總共108頁。隨訪資料的特點1、生存時間的分布呈偏態(tài)。一般為正偏態(tài)，因為生存時間一般為單調(diào)遞減。102名黑色素瘤患者的生存時間資料，做直方圖。2023/3/2917安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前17頁，總共108頁。2023/3/2918安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前18頁，總共108頁。2023/3/2919安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前19頁，總共108頁。2、數(shù)據(jù)中常含有截尾值。如41+（天），表示該病人至少活了41天；雖然不知道確切時間，但此類不完全數(shù)據(jù)包含一定信息，需要生存分析這特殊方法來利用這些信息。組別生存時間（天）A273541+…

…B9980+102…

…2023/3/2920安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前20頁，總共108頁。2023/3/2921安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前21頁，總共108頁。生存分析的主要內(nèi)容及研究方法2023/3/2922安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前22頁，總共108頁。

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

乘積-極限法（Kaplan-Meier）壽命表法（LifeTables）2、生存過程的比較

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

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

2023/3/2923安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前23頁，總共108頁。生存分析對資料中應(yīng)變量的要求：2023/3/2924安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前24頁，總共108頁。1、達(dá)到終點的例數(shù)所占的比例不能太少，即完全數(shù)據(jù)占大部分，截尾值所占的比例要＜10%；2、截尾原因無偏性；3、生存時間盡可能精確。2023/3/2925安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前25頁，總共108頁。常用的研究指標(biāo)2023/3/2926安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前26頁，總共108頁。1、生存函數(shù)（survivalfunction）又稱累積生存概率（cumulativesurvivalprobability）、生存率，表示個體生存時間T≥t的概率。2023/3/2927安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前27頁，總共108頁。生存概率與生存率生存概率：在單位時段開始存活的個體到該時段結(jié)束時仍存活的可能性，用p表示；生存率：觀察個體活過t個單位時間的概率，用S(t)表示。如數(shù)據(jù)中無截尾值，則生存率計算公式為：S(t)=生存時間≥t的病人數(shù)/隨訪病人總數(shù)如有截尾值，則分時段計算不同單位時段的生存概率p1、p2、…pt，S(t)=p1×p2×…×pt，即累積生存概率。2023/3/2928安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前28頁，總共108頁。2023/3/2929安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前29頁，總共108頁?？刹捎肧PSS軟件中生存分析方法計算出生存率大小，并畫出生存函數(shù)圖（或生存率曲線）直接觀察。2023/3/2930安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前30頁，總共108頁。2023/3/2931安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前31頁，總共108頁。2023/3/2932安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前32頁，總共108頁。2、半數(shù)生存時間(mediansurvivaltime)又稱中位生存時間，表示累積生存概率（生存率）為50%的時候所對應(yīng)的生存時間；為中位數(shù)指標(biāo)，因為生存時間的分布是偏態(tài)的。在分析結(jié)果中給出中位生存時間；從生存曲線圖中也能估計出相同結(jié)果。2023/3/2933安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前33頁，總共108頁。2023/3/2934安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前34頁，總共108頁。3、風(fēng)險函數(shù)(hazardfunction)又稱危險率函數(shù)等，一般用h(t)表示，表示生存時間已達(dá)t的個體在t時刻的瞬時風(fēng)險率；在Cox比例風(fēng)險模型中風(fēng)險函數(shù)常用λ(t)表示。h(t)=死于區(qū)間（t，t+⊿t）的病人數(shù)在t時刻尚存的病人數(shù)×⊿t⊿t為極小時間段。2023/3/2935安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前35頁，總共108頁。2023/3/2936安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前36頁，總共108頁。4、風(fēng)險比（hazardratio）是指同一時點兩組的風(fēng)險函數(shù)之比，即相對危險度RR。風(fēng)險比=第一組個體的h1(t)第二組個體的h2(t)

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

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

“Options”對話框“CompareFactors”對話框2023/3/2939安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前39頁，總共108頁。見下例word文檔中的數(shù)據(jù)及對應(yīng)的SPSS數(shù)據(jù)文件KM1.sav，采用乘積極限法進行分析。注意：1）SPSS數(shù)據(jù)文件格式；2）SPSS中乘積極限法的操作過程；3）乘積極限法的具體結(jié)果解釋。2023/3/2940安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前40頁，總共108頁。2023/3/2941安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前41頁，總共108頁。2023/3/2942安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前42頁，總共108頁。2023/3/2943安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前43頁，總共108頁。2023/3/2944安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前44頁，總共108頁。2023/3/2945安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前45頁，總共108頁。2023/3/2946安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前46頁，總共108頁。2023/3/2947安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前47頁，總共108頁。2023/3/2948安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前48頁，總共108頁。注意事項：組間的生存率比較時，要求各組的生存曲線不能交叉，若出現(xiàn)交叉，提示可能存在混雜因素，應(yīng)采用多因素方法來校正混雜因素或分層做統(tǒng)計分析。2023/3/2949安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前49頁，總共108頁。采用Lifetable法對某生存數(shù)據(jù)進行分析2023/3/2950安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前50頁，總共108頁。2023/3/2951安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前51頁，總共108頁。2023/3/2952安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前52頁，總共108頁。生存分析中的多因素分析方法2023/3/2953安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前53頁，總共108頁。2023/3/2954安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前54頁，總共108頁。2023/3/2955安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前55頁，總共108頁。（生存分析方法中單因素與多因素分析的區(qū)別）若只比較男、女的生存情況優(yōu)劣，則可采用前面所述的單因素分析方法——乘積極限法、壽命表法。而在實際工作中，研究者關(guān)心的影響因素不止一個，如在此例中，除了性別，還有年齡等其他因素。因此，應(yīng)采用多因素分析方法。2023/3/2956安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前56頁，總共108頁。（多元線性回歸、多元logistic回歸與生存分析中多因素回歸分析方法的區(qū)別）如果僅以生存時間為反應(yīng)變量做多元線性回歸，則生存時間并不一定服從近似正態(tài)分布，不滿足線性回歸的要求；若僅以是否達(dá)到終點為反應(yīng)變量做多元logistic回歸，則未能充分利用生存時間長短的信息。而且兩者均未考慮到資料中刪失數(shù)據(jù)提供的信息。2023/3/2957安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前57頁，總共108頁。解決措施：采用生存分析中的多因素回歸。生存分析中的多因素回歸同時以生存時間、結(jié)局為反應(yīng)變量，有效利用不完全數(shù)據(jù)提供的信息，從而探索生存的影響因素。2023/3/2958安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前58頁，總共108頁。生存分析中的多因素回歸包括：半?yún)?shù)方法：cox比例風(fēng)險模型、非比例風(fēng)險模型參數(shù)方法：指數(shù)分布模型、Weibull分布模型、Gompertz分布模型、對數(shù)正態(tài)分布模型、對數(shù)logistic分布模型，等。參數(shù)方法中均假設(shè)生存時間服從某分布，而半?yún)?shù)方法中不需指定生存時間的分布情況。2023/3/2959安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前59頁，總共108頁。Cox比例風(fēng)險模型2023/3/2960安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前60頁，總共108頁。適用情況：用于分析帶有伴隨變量的生存時間資料，如：腫瘤和其它慢性病的預(yù)后分析，臨床療效評價和隊列研究的病因探索。優(yōu)點：適用條件寬，便于作多因素分析。2023/3/2961安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前61頁，總共108頁。一、cox回歸模型簡介2023/3/2962安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前62頁，總共108頁。1、模型結(jié)構(gòu)即Cox回歸，它的一般形式是：

λ(t)=[λ0(t)]?

e(β1*x1+β2*x2+….+βm*xm

)

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

)即：ln[λ(t)/λ0(t)]=β1*x1+β2*x2+….+βm*xm將cox回歸與線性回歸、logistic回歸做一比較。2023/3/2964安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前64頁，總共108頁。λ0(t)的形式不限，使模型的適用范圍如同非參數(shù)法那么廣；模型中含有參數(shù)β，又使得其統(tǒng)計效率接近參數(shù)模型。（但參數(shù)β就不能用傳統(tǒng)的方法進行估計和檢驗）∴模型中含有參數(shù)β，但基線風(fēng)險函數(shù)λ0(t)不要求服從特定分布形式，具有非參數(shù)的特點,故Cox回歸又稱為半?yún)?shù)模型。2023/3/2965安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前65頁，總共108頁。英國生物統(tǒng)計學(xué)家D.R.Cox于1972年建立了條件死亡概率和偏似然函數(shù)方法，解決了參數(shù)β的估計和檢驗問題。故，該統(tǒng)計分析方法稱為Cox模型。2023/3/2966安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前66頁，總共108頁。2、參數(shù)的統(tǒng)計學(xué)意義假設(shè)只有一個X危險因素，0為“不暴露”，1為“暴露”，模型形式為：λ(t)=λ0(t)?exp(β*x)X=1時的風(fēng)險函數(shù)為：λ1(t)=λ0(t)?exp(β*1)X=0時的風(fēng)險函數(shù)為：λ2(t)=λ0(t)?exp(β*0)λ1(t)、λ2(t)分別表示暴露于危險因素的兩種不同狀態(tài)下發(fā)病的風(fēng)險函數(shù)。2023/3/2967安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前67頁，總共108頁。暴露于危險因素的兩種不同狀態(tài)下發(fā)病的風(fēng)險比即為相對危險度RR。RR=λ1(t)/λ2(t)=[λ0(t)?exp(β*1)]/[λ0(t)?exp(β*0)]=exp(β*1)/exp(β*0)=exp（β）RR表示暴露組與非暴露組的風(fēng)險函數(shù)之比，即：暴露組發(fā)病的風(fēng)險是非暴露組的RR倍?；貧w系數(shù)β的流行病學(xué)含義是0、1協(xié)變量X的相對危險度的自然對數(shù)。2023/3/2968安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前68頁，總共108頁。如果X為連續(xù)性變量——年齡（歲），則λ1(t)、λ2(t)分別表示k歲、k+1歲的兩種不同狀態(tài)下發(fā)病的風(fēng)險函數(shù)。RR=λ1(t)/λ2(t)=[λ0(t)?exp(β*（k+1）)]/[λ0(t)?exp(β*k)]=exp(β*（k+1）/exp(β*k)=exp（β）β：X每增加一個單位時其相對危險度的自然對數(shù)值。RR表示某危險因素（或協(xié)變量）改變一個測量單位時風(fēng)險函數(shù)改變多少倍。2023/3/2969安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前69頁，總共108頁。當(dāng)β＞0，說明相應(yīng)協(xié)變量值的增加將增大所研究事件發(fā)生的可能性；當(dāng)β＜0，相應(yīng)協(xié)變量值的增加將減少所研究事件發(fā)生的可能性；當(dāng)β=0，相應(yīng)協(xié)變量與所研究事件的發(fā)生無關(guān)。2023/3/2970安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前70頁，總共108頁。3、模型假設(shè)滿足比例風(fēng)險（proportionalhazards）假定，簡稱PH假定。即：相對危險度RR或風(fēng)險比保持一個恒定的比例，與時間t無關(guān)。2023/3/2971安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前71頁，總共108頁。Cox回歸無須對λ0(t)的形式作任何限制，因為不影響各危險因素相對危險度的估計，而相對危險度正是多因素分析時最關(guān)注的問題?！郼ox模型巧妙地將非參數(shù)[λ0(t)]部分與參數(shù)（回歸系數(shù)β）的概念結(jié)合起來，這種靈活性使得它在生存分析的應(yīng)用中備受青睞。2023/3/2972安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前72頁，總共108頁。比值比OR可以作為相對危險度RR的估計值。2023/3/2973安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前73頁，總共108頁。比例風(fēng)險假設(shè)的檢驗：可分組（若為0、1兩組）做二次對數(shù)生存率圖（LogminusLog，LML），進行目測判斷，若兩條曲線有交叉或間距變動大，則提示風(fēng)險比隨時間而改變，不具備cox比例風(fēng)險模型的適用性；反之，則資料適宜用cox比例風(fēng)險模型。2023/3/2974安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前74頁，總共108頁。2023/3/2975安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前75頁，總共108頁。4、協(xié)變量的篩選策略與其他回歸模型類似，如果協(xié)變量個數(shù)較多，則通常采用逐步法。2023/3/2976安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前76頁，總共108頁。二、cox回歸模型的分析步驟及主要結(jié)果2023/3/2977安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前77頁，總共108頁。SPSS中的程序為：Analyze→Survival→CoxRegression將cox回歸的對話框與線性回歸、Kaplan-Meier、Lifetable做比較。見文件cox1.sav2023/3/2978安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前78頁，總共108頁。示例：評價A、B兩治療方案對某病的治療效果。A組（group=0）12人，B組（group=1）13人；腎功能正常者為0，不正常者為1。應(yīng)變量是治療后的生存時間和觀察結(jié)果（達(dá)到終點還是截尾值）；group是研究因素x1，腎功能是混雜因素x2。2023/3/2979安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前79頁，總共108頁。注意事項一般情況下，截尾值設(shè)為0，完全數(shù)據(jù)（即出現(xiàn)終點事件）為1；但在具體資料中，沒有限制，只要能有兩個分別表示完全數(shù)據(jù)和截尾的數(shù)值，結(jié)果不受影響。但是一旦設(shè)定后，分析時，“Status”對話框中的終點事件值一定要和資料中吻合，否則結(jié)果錯誤。2023/3/2980安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前80頁，總共108頁。SPSS分析結(jié)果如下：第一步、求Cox比例風(fēng)險模型2023/3/2981安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前81頁，總共108頁。2023/3/2982安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前82頁，總共108頁。λ(t)=λ0(t)·e(1.243group+4.105kidney)=λ0(t)·exp(1.243group+4.105kidney)group、kidney的P值分別為0.038、0.000，均＜0.05，說明不同治療方案及腎功能對病人的生存有影響。2023/3/2983安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前83頁，總共108頁。group、kidney的OR值分別為3.466、60.670，均大于1，說明接受B治療方案的病人死亡的風(fēng)險是接受A治療方案者的3.466倍、腎功能不正常的病人死亡的風(fēng)險是腎功能正常者的60.670倍。2023/3/2984安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前84頁，總共108頁。第二步、由模型可求出不同狀態(tài)下的相對危險度RR值2023/3/2985安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前85頁，總共108頁。1、腎功能正常者接受B治療方案比接受A治療方案在某時刻死亡的相對危險度為：RR=[λ0(t)]?

e(1.243*1+4.105*0)

[λ0(t)]?e(1.243*0+4.105*0)=3.4662023/3/2986安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前86頁，總共108頁。2、腎功能不正常者接受B治療方案比接受A治療方案在某時刻死亡的相對危險度為：RR=[λ0(t)]?e(1.243*1+4.105*1)

[λ0(t)]?e(1.243*0+4.105*1)=3.4662023/3/2987安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前87頁，總共108頁。3、腎功能不正常者接受B治療方案，比腎功能正常者接受A治療方案在某時刻死亡的相對危險度為：RR=[λ0(t)]?e(1.243*1+4.105*1)

[λ0(t)]?e(1.243*0+4.105*0)=210.300可見在Cox模型的建立中沒有估計基線風(fēng)險λ0(t)；在危險因素的RR分析中也沒有涉及基線風(fēng)險。所估計的風(fēng)險比與時間無關(guān)。2023/3/2988安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前88頁，總共108頁。第三步、畫出風(fēng)險函數(shù)曲線2023/3/2989安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前89頁，總共108頁。2023/3/2990安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前90頁，總共108頁。第四步、畫出二次對數(shù)生存率圖（LogminusLog，LML）從LML圖中，目測判斷PH假設(shè)是否成立。2023/3/2991安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前91頁，總共108頁。2023/3/2992安徽醫(yī)科大學(xué)流統(tǒng)系王靜制作當(dāng)前92頁，總共108頁。目的:探討影響胃癌患者長期生存的預(yù)后因素。方法:采用Cox比例風(fēng)險模型對可能影響胃癌患者預(yù)后的11項臨床病理參數(shù)進行單因素和多因素回歸分析。結(jié)果:Cox單因素回歸分析表明胃癌的分化程度、浸潤深度、淋巴結(jié)轉(zhuǎn)移數(shù)、脈管內(nèi)是否有癌栓、臨床分期、腫瘤直徑、淋巴細(xì)胞浸潤量與胃癌的預(yù)后有關(guān)(P<0.05),而性別、