R cinsinden ayrık zamanlı tehlike modelleri (cloglog)

survivalPaket Rgöründüğünü sürekli zaman hayatta kalma modellerine odaklanmak. Tamamlayıcı log-log modeli olan oransal bir tehlike modelinin ayrık zamanlı bir versiyonunu tahmin etmekle ilgileniyorum. Basit sansür ile oldukça basit bir hayatta kalma modelim var.

Bu modeli tahmin etmenin bir yolunun, "ölü" olmadığı her dönem için her gözlem için ayrı bir satır içeren bir veri kümesi oluşturmak olduğunu biliyorum. Ardından, bağlantıya glmsahip bir model cloglogkullanılabilir.

Bu yaklaşım çok bellek yetersiz kalıyor; aslında, muhtemelen makinemdeki bellek için çok büyük bir veri kümesi üretecektir.

İkinci bir yaklaşım MLE'yi kendim kodlamak olacaktır. Bu yeterince basit olurdu, ancak bu hayatta kalma modelinin hazır olduğu bir paket olduğunu umuyorum. İşbirliği ve kodlama hatalarının bir paket kullanmasını önlemek daha kolay olurdu.

Böyle bir paket bilen var mı?

Her bir gözlem için birkaç sıraya sahip olmak gereksiz görünebilir, ancak büyük olasılıkla değildir. Modelde zamanla değişen değişkenler varsa, her gözlem ayı mutlaka kendi satırına ihtiyaç duyacaktır. Zamanla değişen bir değişkenin belirli bir örneği geçen süredir. Bu değişken neredeyse kesinlikle modele dahil edilmesi gerektiğinden, her gözlem dönemi için ayrı bir satır olması mantıklıdır. Bu nedenle, önerilen ilk yaklaşım muhtemelen en iyisidir.

Bunun, Weibull dağılımına sahip sürekli zaman orantılı bir tehlike modelinden farklı olduğunu unutmayın. Burada, geçen zaman sadece zamanla değişen tek değişken ise, hayatta kalma modeli her bir gözlem için tek bir satıra sadeleştirilebilir ( örneğin buraya bakınız ). Benzer bir sonuç Cox orantılı tehlike modeli için de geçerlidir.