生存分析報告地cox回歸模型案例-spss

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

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

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

4、術(shù)中放療)2分類0,無術(shù)中放擰Jt育術(shù)中放療T.f性別:!分類。二男.1:女I MU占位處2分更0,腰臟頭郃頭部以外tJ聯(lián)膽苗浸潤程度一庫片分類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對蘆度超標(biāo)準(zhǔn)caseno 數(shù)值(N) 80無無8 差右 /度量俗)time數(shù)值(N)81生存時間月)Is無36time數(shù)值(N)81生存時間月)Is無36度量網(wǎng)censsr數(shù)值(N)80刪關(guān)0.死亡 .無8W右*名義伊)無他函ESX數(shù)值(N) 80

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

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

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

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

9、5拒絕原假設(shè)，認(rèn)為有偏回歸系數(shù)不為零的因素，值得進(jìn) 一步分析。實用文檔方物中國變玷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一國1. 103這是多元回歸結(jié)果，第二列 B為偏回歸系數(shù)，最后三列為 OR值及其置信區(qū)問 由P值可以看出，在0.5的顯著水平下，只有trt有統(tǒng)計學(xué)差異，OR為2.265曲變流均闡均伯61.133sextrt.265bui.217.366絢.09&th(3J.205P,723stae.53Q這是協(xié)變量的平均值。實用文檔協(xié)變量均值處的生存函數(shù)生存時M C月）這是總體的生存函數(shù)，即累積生存率函數(shù)實用文檔O模式1