การประมาณค่ามัธยฐานการรอดชีพจากโมเดล Weibull
ตอนนี้เราจะประมาณค่าการรอดชีพของผู้ป่วยมะเร็งเต้านมในชุดข้อมูล GBSG2 โดยใช้โมเดล Weibull (ฟังก์ชัน survreg())
จำไว้ว่าโมเดล Weibull ประมาณฟังก์ชันการรอดชีพแบบเส้นโค้งต่อเนื่อง แทนที่จะเป็นฟังก์ชันแบบขั้นบันไดอย่างวิธี Kaplan-Meier
ฟังก์ชัน predict() ที่ใช้ type = "quantile" ช่วยให้คำนวณควอนไทล์ของฟังก์ชันการแจกแจงได้
เราจะใช้สิ่งนี้เพื่อคำนวณค่ามัธยฐานการรอดชีพ
แพ็กเกจ survival และชุดข้อมูล GBSG2 ถูกโหลดไว้ให้แล้วในแบบฝึกหัดนี้
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Survival Analysis in R
คำแนะนำการฝึกหัด
- สร้างโมเดล Weibull สำหรับผู้ป่วยมะเร็งเต้านม
- คำนวณค่ามัธยฐานการรอดชีพจากโมเดลนี้โดยใช้ฟังก์ชัน
predict()กับtype = "quantile"
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Weibull model
wb <- ___(Surv(time, cens) ~ 1, data = GBSG2)
# Compute the median survival from the model
predict(wb, type = "quantile", p = ___, newdata = data.frame(1))