คำนวณโมเดล Weibull
สำหรับโมเดล Weibull ที่มีตัวแปรร่วม เราสามารถคำนวณ survival curve ได้เช่นเดียวกับโมเดล Weibull ที่ไม่มีตัวแปรร่วม สิ่งที่ต้องทำเพิ่มเติมคือระบุค่าของตัวแปรร่วมสำหรับ survival curve ที่ต้องการในฟังก์ชัน predict() ซึ่งทำได้ผ่านอาร์กิวเมนต์ newdata
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Survival Analysis in R
คำแนะนำการฝึกหัด
- สร้างโมเดล Weibull สำหรับชุดข้อมูล
GBSG2โดยใช้ตัวแปรร่วมhorThเพื่อวิเคราะห์ผลของการรักษาด้วยฮอร์โมนต่อการรอดชีวิตของผู้ป่วย - คำนวณ survival curve สำหรับผู้ป่วยที่ได้รับการรักษาด้วยฮอร์โมน
- ตรวจสอบ survival curve ด้วย
str()
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Weibull model
wbmod <- survreg(Surv(___, ___) ~ ___, data = GBSG2)
coef(wbmod)
# Retrieve survival curve from model
surv <- seq(.99, .01, by = -.01)
t_yes <- predict(___, type = "quantile", p = ___,
newdata = data.frame(horTh = ___))
# Take a look at survival curve
str(t_yes)