模型假設
你已經在上一題的程式練習中看過 Cox PH 模型。這一題要檢查你的模型是否合適。你的模型仍然儲存在物件 fitCPH 中。
本練習屬於課程
R 的行銷分析機器學習
練習說明
- 使用
cox.zph()檢查模型fitCPH的比例風險假設。將檢定結果存到testCPH，並印出它。 - 在顯著水準 0.05 下，似乎有一個變數違反了假設。是哪一個？請為這個變數繪製係數 beta 隨時間變化的圖。
- 使用
rms套件的validate()函式進行交叉驗證，來驗證模型。
動手互動練習
試著完成這個範例程式碼，體驗一下這個練習。
# Check proportional hazard assumption and print result
___ <- ___(fitCPH)
___(testCPH)
# Plot time-dependent beta
plot(___, var = "___")
# Load rms package
library(rms)
# Validate model
___(fitCPH, method = "___",
B = 10, dxy = TRUE, pr = FALSE)