Các giả định của mô hình
Bạn đã xem qua mô hình Cox PH ở bài tập mã hóa trước. Trong bài này, bạn sẽ kiểm tra xem mô hình của mình có phù hợp không. Mô hình của bạn vẫn đang được lưu trong đối tượng fitCPH.
Bài tập này là một phần của khóa học
Machine Learning cho Phân tích Marketing bằng R
Hướng dẫn bài tập
- Kiểm tra giả định hazard tỷ lệ của mô hình
fitCPHbằngcox.zph(). Lưu kết quả kiểm định vào một đối tượng têntestCPHvà in ra. - Có vẻ giả định bị vi phạm với một biến ở mức alpha 0,05. Đó là biến nào? Hãy vẽ hệ số beta theo thời gian cho biến này.
- Đánh giá mô hình bằng cross validation với hàm
validate()từ góirms.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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)