Các mô hình nguy hiểm thời gian riêng biệt (guốc) trong R

Các survivalgói trong Rdường như tập trung vào các mô hình sống sót thời gian liên tục. Tôi quan tâm đến việc ước tính một phiên bản thời gian riêng biệt của một mô hình rủi ro tỷ lệ, mô hình log-log bổ sung. Tôi có một mô hình sinh tồn khá đơn giản, với sự kiểm duyệt đơn giản.

Tôi biết rằng một cách để ước tính mô hình này là tạo một tập dữ liệu có một hàng riêng cho mỗi lần quan sát cho từng giai đoạn mà nó không "chết". Sau đó, một glmmô hình với cloglogliên kết có thể được sử dụng.

Cách tiếp cận này có vẻ rất kém hiệu quả bộ nhớ; thật vậy, nó có thể sẽ tạo ra một tập dữ liệu quá lớn cho bộ nhớ trong máy của tôi.

Cách tiếp cận thứ hai sẽ là tự mình viết mã MLE. Điều đó đủ đơn giản, nhưng tôi hy vọng rằng có một gói có mô hình sinh tồn này đóng hộp. Nó sẽ chỉ dễ dàng hơn cho sự hợp tác và để tránh các lỗi mã hóa để sử dụng một gói.

Có ai biết một gói như vậy?

Có một số hàng cho mỗi quan sát có vẻ dư thừa, nhưng, có thể, không phải vậy. Nếu có bất kỳ đồng biến thời gian nào trong mô hình, thì mỗi tháng quan sát chắc chắn sẽ cần hàng riêng. Một ví dụ cụ thể về hiệp phương sai thời gian là thời gian trôi qua. Do biến này gần như chắc chắn được đưa vào mô hình, nên có một hàng riêng cho mỗi giai đoạn quan sát. Vì vậy, cách tiếp cận đầu tiên được đề xuất có khả năng là cách tốt nhất.

Lưu ý rằng điều này khác với mô hình mối nguy theo tỷ lệ thời gian liên tục với phân phối Weibull. Ở đó, mô hình sinh tồn có thể được đơn giản hóa thành một dòng duy nhất cho mỗi lần quan sát nếu thời gian trôi qua là hiệp phương sai duy nhất theo thời gian ( ví dụ, xem ở đây ). Một kết quả tương tự giữ cho mô hình nguy cơ tỷ lệ Cox.