ข้อสมมติของโมเดล
ในแบบฝึกหัดการเขียนโค้ดก่อนหน้า คุณได้ทำความรู้จักกับโมเดล Cox PH ไปแล้ว ในแบบฝึกหัดนี้ เราจะมาตรวจสอบว่าโมเดลที่สร้างขึ้นมีความเหมาะสมเพียงใด โดยโมเดลยังคงถูกเก็บไว้ในออบเจกต์ fitCPH
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Machine Learning for Marketing Analytics ใน R
คำแนะนำการฝึกหัด
- ตรวจสอบข้อสมมติ proportional hazard ของโมเดล
fitCPHโดยใช้cox.zph()เก็บผลการทดสอบไว้ในออบเจกต์ชื่อtestCPHแล้วแสดงผล - ดูเหมือนว่าข้อสมมติจะถูกละเมิดสำหรับตัวแปรหนึ่งที่ระดับนัยสำคัญ 0.05 ตัวแปรใดที่ว่า? ให้พล็อตค่าสัมประสิทธิ์ beta ตามเวลาสำหรับตัวแปรนั้น
- ตรวจสอบความถูกต้องของโมเดลโดยใช้ cross validation ผ่านฟังก์ชัน
validate()จากแพ็กเกจrms
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# 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)