เปรียบเทียบโมเดล Weibull กับค่าประมาณ Kaplan-Meier
มาพล็อต survival curve ที่ได้จากโมเดล Weibull สำหรับข้อมูล GBSG2 กัน!
ในแบบฝึกหัดนี้ แพ็กเกจ survival และ survminer รวมถึงข้อมูล GBSG2 ถูกโหลดไว้ให้แล้ว
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Survival Analysis in R
คำแนะนำการฝึกหัด
- คำนวณโมเดล Weibull สำหรับข้อมูล
GBSG2 - คำนวณ survival curve จากโมเดล
- พล็อต survival curve ที่ได้จากทั้งสองค่าประมาณ
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Weibull model
wb <- ___(___(___, ___) ~ ___, ___)
# Retrieve survival curve from model
surv <- seq(.99, .01, by = -.01)
# Get time for each probability
t <- predict(___, type = ___, p = ___, newdata = data.frame(___))
# Create data frame with the information needed for ggsurvplot_df
surv_wb <- data.frame(time = ___, surv = ___,
upper = NA, lower = NA, std.err = NA)
# Plot
ggsurvplot_df(fit = ___, surv.geom = ___)