生存分析報(bào)告地cox回歸模型案例-spss

1、實(shí)用文檔一、生存分析基本概念1、事件(Event)指研究中規(guī)定的生存研究的終點(diǎn)，在研究開(kāi)始之前就已經(jīng)制定好。根據(jù)研究性 質(zhì)的不同，事件可以是患者的死亡、疾病的復(fù)發(fā)、儀器的故障，也可以是下崗 工人的再就業(yè)等等。2、生存時(shí)間(Survival time)指從某一起點(diǎn)到事件發(fā)生所經(jīng)過(guò)的時(shí)間。生存是一個(gè)廣義的概念，不僅僅指醫(yī) 學(xué)中的存活，也可以是機(jī)器出故障前的正常運(yùn)行時(shí)間，或者下崗工人再就業(yè)前 的待業(yè)時(shí)間等等。有的時(shí)候甚至不是通用意義上的時(shí)間，比如汽車在出故障前 的行駛里程，也可以作為生存時(shí)間來(lái)考慮。3、刪失(Sensoring )指由于所關(guān)心的事件沒(méi)有被觀測(cè)到或者無(wú)法觀測(cè)到，以至于生存時(shí)間無(wú)法記錄的

2、情況。常由兩種情況導(dǎo)致：(1)失訪；(2)在研究終止時(shí)，所關(guān)心的事件 還未發(fā)生。4、生存函數(shù)(Survival distribution function )又叫累積生存率，表達(dá)式為 S (t) =P(Tt),其中T為生存時(shí)間，該函數(shù)的意義 是生存時(shí)間大于時(shí)間點(diǎn)t的概率。t=0時(shí)S(t)=1 ,隨著t的增加S(t)遞減(嚴(yán) 格的說(shuō)是不增)，1-S(t)為累積分布函數(shù)，表示生存時(shí)間 T不超過(guò)t的概率。二、生存分析的方法1、生存分析的主要目的是估計(jì)生存函數(shù)，常用的方法有Kaplan-Meier法和壽命表法。對(duì)于分組數(shù)據(jù)，在不考慮其他混雜因素的情況下，可以用這兩種方法 對(duì)生存函數(shù)進(jìn)行組間比較。2、如

3、果考慮其他影響生存時(shí)間分布的因素，可以使用Cox回歸模型(也叫比例風(fēng)險(xiǎn)模型)，利用數(shù)學(xué)模型擬合生存分布與影響因子之間的關(guān)系，評(píng)價(jià)影響因 子對(duì)生存函數(shù)分布的影響程度。這里的前體是影響因素的作用不隨時(shí)間改變， 如果不滿足這個(gè)條件，則應(yīng)使用含有時(shí)間依存協(xié)變量的Cox回歸模型。卜面用一個(gè)例子來(lái)說(shuō)明SPSS Cox回歸模型的操作方法。例題要研究胰腺癌術(shù)中放療對(duì)患者生存時(shí)間的影響，收集了下面所示的數(shù)據(jù):實(shí)用文檔委1丈3腴苴1癌術(shù)中放擰效果研究數(shù)據(jù)的說(shuō)明變量名變埴說(shuō)明變麻類型分類轉(zhuǎn)的端附UHSTIH1患者編號(hào)hriw生一時(shí)同CM)連續(xù)ve川甘雕失2分關(guān)0：死亡J：洲失H卯手術(shù)時(shí)的年的連續(xù)trt處理組別(仃無(wú)

4、術(shù)中放療)2分類0,無(wú)術(shù)中放擰Jt育術(shù)中放療T.f性別:!分類。二男.1:女I MU占位處2分更0,腰臟頭郃頭部以外tJ聯(lián)膽苗浸潤(rùn)程度一庫(kù)片分類lltOsilil 3：t2.4r J13P行天腹膜轉(zhuǎn)移2分類0i3；J j 有嗎I力第LTNH分類2分類3dii 期、4； n .操作步驟：SPS皴量視圖名稱 I：類j1寬度F扇1標(biāo)噴 II 值 R ,列J對(duì)蘆度超標(biāo)準(zhǔn)caseno 數(shù)值(N) 80無(wú)無(wú)8 差右 /度量俗)time數(shù)值(N)81生存時(shí)間月)Is無(wú)36time數(shù)值(N)81生存時(shí)間月)Is無(wú)36度量網(wǎng)censsr數(shù)值(N)80刪關(guān)0.死亡 .無(wú)8W右*名義伊)無(wú)他函ESX數(shù)值(N) 80

5、手術(shù)時(shí)的年齡5 h 嗎 trtbuchpst數(shù)值(N)8Q有無(wú)術(shù)中放療Q無(wú)術(shù)中放黝值(N)80占位處但胰臟頭茴數(shù)值(N)80膜膽管it制程度乩CH0數(shù)值網(wǎng))80有無(wú)腹膜轉(zhuǎn)移也無(wú)數(shù)值(忖)80TN博分期p. nm.數(shù)值(忖)80性別無(wú)無(wú)無(wú)無(wú)要右型度量(S)專有 *名義網(wǎng)A名義州) 晶名叉串) d序號(hào)(。) 名義N) 46名尺叫右右右右右三二三.二號(hào)基菜單選擇:實(shí)用文檔點(diǎn)擊進(jìn)入Cox主對(duì)話框，如下，將time選入“時(shí)間”框，將代表刪失的censor 變量選入“狀態(tài)”框，其余分析變量選入“協(xié)變量”框。其余默認(rèn)就行。實(shí)用文檔點(diǎn)擊“狀態(tài)”框下方的“定義事件”，將事件發(fā)生的標(biāo)志設(shè)為值 0,即0代表事 件發(fā)

6、生。在主對(duì)話框中點(diǎn)擊“分類”按鈕，進(jìn)入如下的對(duì)話框，將所有分類變量選入右 邊框中。實(shí)用文檔詼C”回歸：定義分類班變量4手術(shù)時(shí)的年般國(guó)羹分類協(xié)重量立）： 由的示符） 眄【指示新） 向繾示用 P舸符 M能示用 小英脂示符:4手術(shù)時(shí)的年般國(guó)羹.更仁對(duì)比世上（1叵畫(huà)耗考空另才斷鼎后一個(gè)心）第一個(gè)任）嫌取消j國(guó)記在主對(duì)話框中點(diǎn)擊“繪圖”按鈕，進(jìn)入如下的對(duì)話框，選擇繪圖的類型，這里 只選擇“生存函數(shù)”。由于我們關(guān)心的主要變量是trt （是否放療），所以將trt選入“單線”框中，繪制生存曲線。實(shí)用文檔在主對(duì)話框中點(diǎn)擊“選項(xiàng)”按鈕，進(jìn)入如下的對(duì)話框，設(shè)置如下，輸出RR的95% 置信區(qū)間?；氐街鹘缑?，點(diǎn)擊“確定

7、”輸出結(jié)果。結(jié)果輸出案例處理摘要R口分比分析中可用荊案例事件a州失11.2%合計(jì)S3100.0K刪除帆維例帶有賦欠值的案例00.0%帶有魚(yú)時(shí)間的里颯00. 0K尼中晌堆早事件之而則央的案例00. D%00.0%管計(jì)S3100.09ft，囚爻最：欠存時(shí)間1月）這是案例處理摘要，有一個(gè)刪失數(shù)據(jù)實(shí)用文檔分突變顯編班/冏打?qū)胔蜒率(l)cG)se/。=男4S12女350trtb。二丈術(shù)中放住5211=時(shí)術(shù)中放療G】0bn/。二腆幀頭部IS1仁弘那以外650chb1=CHO321cD2=CH160i0=CH217c014 二 CH326000pb&01i=W23astaeb 3=ELI期44i年工刊周3

8、9a分英變狼；百。（性別）b.示性參數(shù)國(guó)儲(chǔ)明 己歌記求了 SJ）奘玨,所以共系版不含與特示符M） 御的相同.L分類變量：（有光術(shù)中械瘴電,分類變量：bui （占仲，處）.分類嚏量：ch （朝朋笞浸潤(rùn)程度）分類變宛：（存無(wú)勝服陂感）h,分類空量；stag* （T刪分期）這是分類變量的編碼方式。桃烈祭敕的鐮到秋口-21時(shí)做械舞作1格體（和分）從上加歌開(kāi)蚣史&從上一耿”始肥改h hdf*書(shū)*方4Sif機(jī)方arSi55?IC 26093 士】6 1 169M4if； 1 ieg一曲口.劉購(gòu)以編口L方法三臨入這是對(duì)擬合模型的檢驗(yàn)，原假設(shè)是“所有影響因素的偏回歸系數(shù)均為0，這里可以看出P=0.0320.0

9、5拒絕原假設(shè)，認(rèn)為有偏回歸系數(shù)不為零的因素，值得進(jìn) 一步分析。實(shí)用文檔方物中國(guó)變玷EstlaiddfSi耳,e95.0* CI 用 P下部上部.016.CJ31.6631.172I.01&一由睚1.0d5e3E.03424E.01911.036d01.673trt.016.3266. 2ftB1.0122. 2fi5J. 16bni-.599.3502. 93。10B6.5HS.2771.090ch4一 6793.197ch (1)一.135.359.M21一 m.874.433i.765thQ).697.1632. 2601.1332.007S091.976uk (3).4541. 37212-ii1.575.73-3.366P.360.276L 66冷I.191.6%野.4041.205stage. 453.2612. 5961.636一國(guó)1. 103這是多元回歸結(jié)果，第二列 B為偏回歸系數(shù)，最后三列為 OR值及其置信區(qū)問(wèn) 由P值可以看出，在0.5的顯著水平下，只有trt有統(tǒng)計(jì)學(xué)差異，OR為2.265曲變流均闡均伯61.133sextrt.265bui.217.366絢.09&th(3J.205P,723stae.53Q這是協(xié)變量的平均值。實(shí)用文檔協(xié)變量均值處的生存函數(shù)生存時(shí)M C月）這是總體的生存函數(shù)，即累積生存率函數(shù)實(shí)用文檔O模式1