在本篇文中我們將會簡單介紹存活分析中的Cox regression model(Cox proportional hazard model),用以分析會顯著影響死亡率的變數,以下詳細說明。
一、使用狀況
在過去的文章中,我們已經說明Life-table method的用途,以及Kaplan-Meier method與其檢定方式log-rank test跟 Wilcoxon test,這些方法都用於簡單的組別比較,屬於無母數的分析。
如果想要同時校正其他因素(干擾因子),我們可以使用有母數分析的預測方式,比較常見的是Cox proportional hazard model,此模型用於探討影響存活率的重要因子。
二、前提假設與理論說明
在探討Cox regression時,存活函數及風險函數這兩個部分要先釐清。
(一)存活函數 (Survival Function):
t 是時間點
S(0)= 1
S(∞) = 0
存活函數呈現的是在某特定時間點之下,個案可以活過特定時間點的機率是多少。
(二)風險函數 (Hazard Function):
風險函數h(t)值不隨時間變動。
t 是時間點,風險函數h(t)值與存活時間有關。
舉例來說,如果t=3個月,h(t)就是事件機率(死亡)除以3個月時的生存函數。因為t越大,存活函數S(t)越小(存活函數為遞減函數),假定死亡機率密度f(t)對任何時間點一樣,那麼風險函數h(t)值越大,預期存活時間短。
在Cox proportional hazard model之下
λ為一定值,由此可推導出Cox proportional hazard model為:
其Hazard Ratio計算方式如下:
假設今天只有一個變項x1,x1=1(treatment)跟x1=0(placebo),則風險比會如下:
三、SPSS 操作Example
(一)在SPSS中輸入欲分析之資料。
目前有14位分別使用不同藥物的兩組胰臟癌病人,我們追蹤了12個月,獲得他們到過世前存活時間的資料,我們利用這筆資料,分析胰臟癌的存活率。
Group表示治療組:第一組為新藥,第二組為舊藥;Sex表示性別:第一組為男性,第二組為女性。
(二)分析→存活分析→Cox回歸
(三)Cox回歸
時間:點選Survival time。
狀態:選event的變項,定義事件要定義什麼樣的值為event,所以填1。
共變量:選變項group跟sex。
另外在右上角的「類別」裡把「共變量」的將group跟sex放入「類別共變量」並在比對選「指標」,「參考類別」看你想要以哪個為Reference group,我選最後一個是把舊藥當成參考。
圖形裡面,勾選「存活分析」,「圖線繪製依據」選group
(四)結果:
結果如上表,我們可將他整理如下:
由上表可以看出,控制了性別之後,第一組新藥相較於第二組的舊藥的死亡風險比為0.107,且有達到顯著上的差異(p=0.046),但在性別上未達統計上顯著的差異(p=0.183)。由此可知,若控制住其他醫療或是個體差異,則新藥的效果顯著了降低了病人的死亡風險(實務分析中做這樣的推論則必須更為謹慎)。
(五)繪圖
點選「圖形」就會輸出Kaplan-Meier法的生存函數圖,由此圖也可知第一組存活的時間較長,死亡數目較少。