บทสรุป: Cox Model
เพื่อปิดท้ายคอร์สนี้ เรามาย้อนกลับไปดูชุดข้อมูลมะเร็งปอดที่ใช้ไปในช่วง 2 บทที่ผ่านมากัน ชุดข้อมูลนี้บันทึกข้อมูลการรอดชีวิตของผู้ป่วยมะเร็งปอดระยะลุกลามจาก North Central Cancer Treatment Group เหตุการณ์ที่สนใจเก็บอยู่ในตัวแปร status โดยค่า 2 หมายความว่าผู้ป่วยรายนั้นไม่รอดชีวิต คะแนนสมรรถภาพ (ตัวแปร performance) วัดความสามารถในการทำกิจกรรมประจำวันของผู้ป่วย (แย่=0, ดี=100) ประเมินโดยแพทย์ เราต้องการทราบความสัมพันธ์ระหว่างคะแนนสมรรถภาพในระดับต่าง ๆ กับระยะเวลาการรอดชีวิต
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Survival Analysis in R
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Compute Cox model and survival curves
cxmod <- ___(___(___, ___) ~ ___, data = lung)
new_lung <- data.frame(performance = ___)
cxsf <- ___(cxmod, data = ___, newdata = ___, conf.type = "none")