แบบจำลองเวลาอันตรายที่ไม่ต่อเนื่อง (cloglog) ใน R

survivalแพคเกจในการRปรากฏขึ้นที่จะมุ่งเน้นในรูปแบบเวลาการอยู่รอดอย่างต่อเนื่อง ฉันสนใจที่จะประเมินรุ่นเวลาที่ไม่ต่อเนื่องของโมเดลความเป็นอันตรายตามสัดส่วนซึ่งเป็นรูปแบบบันทึกการใช้งานเสริม ฉันมีรูปแบบการเอาชีวิตรอดที่ค่อนข้างตรงไปตรงมาพร้อมการตัดสิทธิ์ที่เรียบง่าย

ฉันรู้ว่าวิธีหนึ่งในการประเมินโมเดลนี้คือการสร้างชุดข้อมูลที่มีแถวแยกต่างหากสำหรับแต่ละการสังเกตสำหรับแต่ละช่วงเวลาที่ไม่ใช่ "ตาย" จากนั้นสามารถใช้glmโมเดลที่มีcloglogลิงก์

วิธีการนี้ดูเหมือนไม่มีประสิทธิภาพของหน่วยความจำมาก ที่จริงแล้วมันน่าจะสร้างชุดข้อมูลที่ใหญ่เกินไปสำหรับหน่วยความจำในเครื่องของฉัน

วิธีที่สองคือการเขียนโค้ด MLE ด้วยตัวเอง นั่นจะง่ายพอ แต่ฉันหวังว่าจะมีแพ็คเกจที่มีโมเดลการเอาตัวรอดนี้ มันจะง่ายขึ้นสำหรับการทำงานร่วมกันและเพื่อหลีกเลี่ยงข้อผิดพลาดในการเขียนโปรแกรมเพื่อใช้แพคเกจ

ไม่มีใครรู้แพคเกจดังกล่าวหรือไม่

มีหลายแถวสำหรับการสังเกตแต่ละครั้งอาจดูเหมือนซ้ำซ้อน แต่เป็นไปได้ว่าไม่ใช่ หากมีตัวแปรแปรผันตามเวลาในโมเดลแต่ละเดือนการสังเกตการณ์จะต้องมีแถวของมันเอง ตัวอย่างหนึ่งของตัวแปรแปรผันตามเวลาคือเวลาที่ผ่านไป เนื่องจากตัวแปรนี้ควรจะรวมอยู่ในโมเดลอย่างแน่นอนจึงเหมาะสมที่จะมีแถวแยกต่างหากสำหรับแต่ละช่วงเวลาการสังเกต ดังนั้นวิธีแรกที่แนะนำน่าจะเป็นวิธีที่ดีที่สุด

โปรดทราบว่านี่เป็นสิ่งที่แตกต่างจากโมเดลอันตรายตามสัดส่วนเวลาต่อเนื่องที่มีการแจกแจงแบบ Weibull ที่นั่นแบบจำลองการเอาตัวรอดสามารถลดความซับซ้อนลงในบรรทัดเดียวสำหรับการสังเกตแต่ละครั้งหากเวลาที่ผ่านไปเป็นค่า covariate ที่แปรผันตามเวลาเพียงครั้งเดียว ( ตัวอย่างเช่นที่นี่ ) ผลลัพธ์ที่คล้ายกันนี้จัดทำขึ้นสำหรับโมเดลอันตรายตามสัดส่วนของ Cox