Modele dyskretnych zagrożeń czasowych (cloglog) w R

Wydaje się, że survivalpakiet Rkoncentruje się na modelach ciągłego przetrwania w czasie. Jestem zainteresowany oszacowaniem dyskretnej wersji proporcjonalnego modelu hazardu w czasie, komplementarnego modelu log-log. Mam dość prosty model przetrwania z prostą cenzurą.

Wiem, że jednym ze sposobów oszacowania tego modelu jest utworzenie zestawu danych, który ma osobny wiersz dla każdej obserwacji dla każdego okresu, w którym nie jest „martwy”. Następnie można użyć glmmodelu z cloglogłączem.

To podejście wydaje się bardzo nieefektywne pod względem pamięci; w rzeczy samej, prawdopodobnie stworzyłby zestaw danych, który jest zbyt duży dla pamięci na moim komputerze.

Drugim podejściem byłoby samodzielne kodowanie MLE. To byłoby dość proste, ale mam nadzieję, że istnieje pakiet z tym modelem przetrwania. Byłoby po prostu łatwiej współpracować i uniknąć błędów w kodowaniu, aby użyć pakietu.

Czy ktoś wie o takiej paczce?

Posiadanie kilku wierszy dla każdej obserwacji może wydawać się zbędne, ale prawdopodobnie tak nie jest. Jeśli w modelu występują zmienne towarzyszące w czasie, wówczas każdy miesiąc obserwacji z pewnością będzie potrzebował własnego wiersza. Jednym z konkretnych przykładów zmiennych towarzyszących w czasie jest upływający czas. Ponieważ zmienna ta prawie na pewno powinna zostać uwzględniona w modelu, sensowne jest posiadanie osobnego wiersza dla każdego okresu obserwacji. Zatem pierwsze sugerowane podejście jest prawdopodobnie najlepsze.

Należy zauważyć, że różni się to od modelu zagrożeń proporcjonalnego do ciągłego czasu z rozkładem Weibulla. Tam model przeżycia można uprościć do jednej linii dla każdej obserwacji, jeżeli upływający czas jest jedyną zmienną zmienną w czasie (patrz tutaj na przykład). Podobny wynik dotyczy proporcjonalnego modelu hazardu Coxa.