Cox模型的基本假設(shè)

Cox比例風(fēng)險(xiǎn)模型的前提條件是：在某因素下，個(gè)體的風(fēng)險(xiǎn)比例在不同時(shí)間點(diǎn)都是定值。

Cox PH假設(shè)檢驗(yàn)

一般我們使用Schoenfeld殘差來(lái)檢驗(yàn)比例風(fēng)險(xiǎn)假設(shè)（還可以簡(jiǎn)單的通過(guò)生存曲線是否交叉來(lái)粗略判斷下），如果數(shù)據(jù)滿足Cox回歸模型，那么Schoenfeld殘差和生存時(shí)間的秩次之間無(wú)明顯的相關(guān)性，也就是說(shuō)Schoenfeld殘差與生存時(shí)間無(wú)明顯的一致變化趨勢(shì)，這里通過(guò)survival包的cox.zph()函數(shù)即可判斷。

library("survival")

library("survminer")

fit <- coxph(Surv(time, status) ~ age + sex, data =? data.set)

test.ph <- cox.zph(fit)

從檢驗(yàn)結(jié)果test.ph中協(xié)變量的P值可進(jìn)行判斷，這里的原假設(shè)為滿足Cox PH假設(shè)，因此，若P<0.05，則拒絕原假設(shè)，即不滿足Cox PH假設(shè)。

此外我們還可以通過(guò)Schoenfeld殘差與時(shí)間的關(guān)系圖來(lái)診斷模型是否滿足前提，理論上如果滿足的話，Schoenfeld殘差應(yīng)該與時(shí)間無(wú)明顯的趨勢(shì)，使用ggcoxzph函數(shù)，即ggcoxzph(test.ph)。

解決方法

若不滿足Cox PH假設(shè)，可以采用如下方法進(jìn)行解決：

Adding covariate*time interaction（增加協(xié)變量與時(shí)間的交互關(guān)系）

Stratification（分層）