ฉันต้องการพัฒนาแบบจำลองการทำนาย (Cox PH) สำหรับการเสียชีวิตแบบทุกสาเหตุในชุดข้อมูลของผู้เข้าร่วมซึ่งเกือบทุกคนเสียชีวิตเมื่อสิ้นสุดการติดตาม (เช่น 1 ปี)
แทนที่จะทำนายความเสี่ยงที่แน่นอนของการตายในเวลาหนึ่งฉันต้องการทำนายเวลาการอยู่รอด (เป็นเดือน) สำหรับแต่ละคน
เป็นไปได้หรือไม่ที่จะได้รับการทำนายเช่นนี้ใน R (จากเช่น coxph-object) และถ้าใช่ฉันจะทำเช่นนั้นได้อย่างไร?
ขอบคุณมากล่วงหน้า!
คำตอบ:
โมเดล Cox Proportional Hazards ไม่ได้จำลองความเป็นอันตรายพื้นฐานซึ่งเป็นสิ่งที่คุณต้องการในการทำนายเวลาการเอาชีวิตรอดแบบนั้น - นี่คือทั้งความแข็งแกร่งที่ยอดเยี่ยมของโมเดลและหนึ่งในข้อเสียที่สำคัญ
หากคุณมีความสนใจเป็นพิเศษในการได้รับการประมาณความน่าจะเป็นของการเอาชีวิตรอดที่จุดเวลาใดเวลาหนึ่งฉันจะชี้ให้คุณไปที่โมเดลการอยู่รอดแบบพารามิเตอร์ สิ่งเหล่านี้ถูกนำมาใช้ในsurvivalแพ็คเกจสำหรับ R และจะให้การแจกแจงเวลาการอยู่รอดแบบอิงพารามิเตอร์ซึ่งคุณสามารถเสียบเข้าไปในเวลาที่คุณสนใจและกลับมาเป็นโอกาสรอดชีวิต
@statBeginner ใช่มันจะ มันต้องมีสองขั้นตอน:
x <- survfit(cox.ph.model, newdata = dataset)
dataset$Results <- summary(x)$table[,"median"]
แต่ฉันไม่แน่ใจว่าเวลามัธยฐานในการเอาชีวิตรอดนั้นแม่นยำเพียงพอหรือไม่
แม้ว่าฉันจะเห็นด้วยกับประเด็นเหล่านี้ความอยู่รอดของค่ามัธยฐานจะเป็นประโยชน์ทางคลินิก
คุณอาจจะสนใจในการทำงานของเรา (และอื่น ๆ ) กำลังมองหาที่ใช้แบ่งเป็นพื้นฐานสำหรับการอยู่รอดช่วงเวลา - เราคิดว่าสิ่งเหล่านี้เป็นประโยชน์มากขึ้น