การเปรียบเทียบโมเดล Weibull และ Log-Normal I
ในแบบฝึกหัดนี้ เราจะเพิ่มค่าความน่าจะเป็นการรอดชีวิตที่ถูกต้องลงใน data frame ซึ่ง data frame นี้จะถูกนำไปใช้พล็อตกราฟเส้นโค้งการรอดชีวิต surv_wide คือ data frame แบบ wide ที่มีข้อมูลการบำบัดด้วยฮอร์โมน รวมถึงเส้นโค้งการรอดชีวิตของโมเดล Weibull และ log-normal
แพ็กเกจ survival, survminer, และ reshape2 รวมถึงข้อมูล GBSG2 ได้ถูกโหลดไว้ให้แล้วในแบบฝึกหัดนี้
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Survival Analysis in R
คำแนะนำการฝึกหัด
- แปลง data frame ให้อยู่ในรูปแบบ long format
- ใช้
surv_long$surv_idเพื่อเพิ่มค่าความน่าจะเป็นการรอดชีวิตsurvที่ถูกต้องลงใน data frame
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Melt the data.frame into long format.
surv_long <- melt(surv_wide, ___ = c("horTh", "dist"), ___ = "surv_id", ___ = "time")
# Add column for the survival probabilities
surv_long$surv <- surv[as.numeric(surv_long$___)]
# Add columns upper, lower, std.err, and strata contianing NA values
surv_long[, c("upper", "lower", "std.err", "strata")] <- NA