เซ็นเซอร์ช่วงแบบ Cox ตามสัดส่วนความเป็นอันตรายใน R

เมื่อกำหนดช่วงเวลาการเอาตัวรอดให้ถูกต้องฉันจะทำการจำลอง Cox PH แบบช่วงเซ็นเซอร์ได้Rอย่างไร การค้นหา rseek จะเปิดแพ็คเกจintcoxซึ่งไม่มีอยู่ในที่Rเก็บอีกต่อไป ฉันเกือบจะเป็นบวกcoxphฟังก์ชั่นในsurvivalแพ็คเกจไม่สามารถจัดการกับข้อมูลการรอดชีวิตแบบเซ็นเซอร์

นอกจากนี้ฉันไม่ต้องการใส่ข้อมูลแล้วใช้coxphฟังก์ชัน วิธีนี้จะประเมินข้อผิดพลาดมาตรฐานต่ำกว่าค่าสัมประสิทธิ์เนื่องจากคุณไม่สนใจความไม่แน่นอนของการเซ็นเซอร์ช่วงเวลา

ตามที่ระบุไว้ข้างต้นคุณสามารถใช้ฟังก์ชัน survreg ข้อควรทราบ: นี่ไม่ได้เป็นรุ่น Cox PH เท่านั้น แต่เป็นรุ่นมาตราส่วนตำแหน่ง ใช้การแปลงค่าเริ่มต้นนี่คือรูปแบบของท้ายเรือ ในกรณีของการแจกแจงแบบเอ็กซ์โพเนนเชียลความแปรปรวนของสัดส่วนและ aft จะเท่ากันดังนั้นหากการแจกแจงถูกตั้งค่าเป็นเลขชี้กำลังนี่เป็นรูปแบบความเป็นอันตรายเชิงสัดส่วนที่มีพื้นฐานเป็นเลขชี้กำลัง ถ้าใช้โมเดลการแจกแจงพื้นฐาน Weibull ท้ายเรือการประมาณพารามิเตอร์เป็นเพียงการแปลงเชิงเส้นของสิ่งที่ใช้ในโมเดลอันตรายตามสัดส่วนด้วยการแจกแจงพื้นฐาน Weibull แต่โดยทั่วไปผู้รอดชีวิตไม่เหมาะกับรุ่น Cox PH

หากต้องการโมเดลกึ่งพารามิเตอร์ตามที่พบนำมาใช้ใน intcox คำเตือน: มีปัญหาหลายประการกับรุ่นปัจจุบันของ intcox (อัลกอริทึมโดยทั่วไปจะสิ้นสุดก่อนกำหนดอย่างมีนัยสำคัญห่างไกลจาก MLE ล้มเหลวทันทีด้วยการสังเกตวาบไม่มีข้อผิดพลาดมาตรฐาน นำเสนอโดยอัตโนมัติ)

ทางเลือกใหม่ที่คุณสามารถใช้ได้คือแพ็คเกจ "icenReg"

ค่าเข้าชมของอคติ: นี่คือผู้เขียนของ icenReg

ในการทำการวิเคราะห์แบบช่วงเซ็นเซอร์ใน R คุณต้องสร้างวัตถุ Surv และใช้ Survfit () หากคุณมีมากกว่าหนึ่งตัวแปรแพคเกจ Intcox จะแก้ปัญหาได้