การประมาณและการแสดงผลเส้นโค้งการรอดชีพ
มาดูข้อมูลการรอดชีพของผู้ป่วยมะเร็งเต้านมกัน
ในแบบฝึกหัดนี้ เราจะทำงานกับชุดข้อมูล GBSG2 อีกครั้ง
แพ็กเกจ survival และ survminer รวมถึงข้อมูล GBSG2 ได้ถูกโหลดไว้ให้แล้ว
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Survival Analysis in R
คำแนะนำการฝึกหัด
- ประมาณฟังก์ชันการรอดชีพสำหรับผู้ป่วยมะเร็งเต้านม
- แสดงผลฟังก์ชันการรอดชีพที่ประมาณได้โดยใช้ฟังก์ชัน
ggsurvplot() - เพิ่มตารางความเสี่ยง (risk table) ลงในกราฟเพื่อแสดงจำนวนผู้ป่วยที่อยู่ภายใต้การสังเกต โดยใช้อาร์กิวเมนต์
risk.table - เพิ่มเส้นแสดงเวลาการรอดชีพมัธยฐานลงในกราฟ โดยใช้อาร์กิวเมนต์
surv.median.line
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Kaplan-Meier estimate
km <- survfit(Surv(___, ___) ~ ___, data = ___)
# Plot of the Kaplan-Meier estimate
ggsurvplot(___)
# Add the risk table to plot
ggsurvplot(___, ___ = TRUE)
# Add a line showing the median survival time
ggsurvplot(___, ___ = TRUE, ___ = "hv")