R'de sansürlenmiş Cox orantılı tehlikeler modeli

Aralık sansürlenmiş hayatta kalma süreleri göz önüne alındığında, aralık sansürlü Cox PH modelini nasıl yapabilirim R? Bir rseek araması intcoxartık Rdepoda bulunmayan paketi açar . Neredeyse paketteki coxphfonksiyonun survivalsansürlenmiş sağkalım verilerini aralıksız işleyemediğine inanıyorum.

Ayrıca, verileri çarpıtmak ve sonra coxphişlevi kullanmak istemiyorum . Bu yöntem, aralık sansürünün belirsizliğini göz ardı ettiğiniz için katsayıların standart hatalarını hafife alır.

Yukarıda belirtildiği gibi, survreg işlevini kullanabilirsiniz. Yine de bir not: bu kesinlikle bir Cox PH modeli değil, daha çok yer ölçekli modeller. Varsayılan günlük dönüşümü kullanılarak, bu kıç modelidir. Üstel dağılım durumunda, oransal tehlikeler ve kıç modeli eşdeğerdir, bu nedenle dağıtım üstel olarak ayarlanırsa, bu, üstel bir taban çizgisine sahip oransal bir tehlike modelidir. Benzer şekilde, bir temel Weibull dağıtım kıç modeli kullanılırsa, parametre tahminleri, Weibull temel dağılımlı oransal tehlike modelinde kullanılanların doğrusal bir dönüşümüdür. Ancak genel olarak, survreg bir Cox PH modeline uymuyor.

Eğer intcox içinde uygulandığı gibi yarı parametrik bir model isteniyorsa, bir uyarı kelimesi: mevcut intcox sürümü ile ilgili birkaç sorun var (algoritma tipik olarak erken MLE'den önemli ölçüde sonlandırılıyor, sansürsüz gözlemlerle açıkça başarısız oluyor, standart hata yok otomatik olarak sunulur).

Kullanabileceğiniz yeni bir alternatif "icenReg" paketidir.

Önyargı kabulü: bu icenReg'in yazarıdır.

R'de sansürlenmiş aralıklarla analiz yapmak için bir Surv nesnesi oluşturmanız ve ardından survfit () kullanmanız gerekir. Birden fazla değişkeniniz varsa, intcox paketi sorunu çözer.