Takip sonunda (örneğin 1 yıl) hepsi ölmüş olan katılımcıların bir veri setinde tüm nedenlere bağlı mortalite için bir tahmin modeli (Cox PH) geliştirmek istiyorum.
Belli bir zaman noktasında ölme mutlak riskini tahmin etmek yerine, her bireyin hayatta kalma süresini (ay olarak) tahmin etmek istiyorum.
R'de (örneğin bir coxph-nesnesinden) bu tür tahminler elde etmek mümkün mü ve eğer evetse, bunu nasıl yapabilirim?
Şimdiden çok teşekkürler!
Yanıtlar:
Cox Orantılı Tehlikeler modeli, altta yatan tehlikeyi modellemez, bu da hayatta kalma süresini böyle tahmin etmek için ihtiyacınız olan şeydir - bu hem modelin büyük gücü hem de büyük dezavantajlarından biridir.
Belirli zaman noktalarında hayatta kalma olasılığına ilişkin tahminler elde etmekle özellikle ilgileniyorsanız, sizi parametrik hayatta kalma modellerine (yani hızlandırılmış başarısızlık süresi modelleri) yönlendiririm. Bunlar survivalR için pakete uygulanır ve size parametrik hayatta kalma süresi dağılımları verir, burada ilgilendiğiniz zamanı takıp hayatta kalma olasılığını geri alabilirsiniz.
@statBeginner Evet olacak. İki adım gerektirir:
x <- survfit(cox.ph.model, newdata = dataset)
dataset$Results <- summary(x)$table[,"median"]
fakat ortalama hayatta kalma süresinin yeterince doğru olup olmadığından emin değilim.
Bu noktaya katılmama rağmen, medyan sağkalım klinik olarak faydalıdır.
Medyanı hayatta kalma aralıkları için bir temel olarak kullanmaya çalışan çalışmalarımızla (ve diğerleriyle) ilgilenebilirsiniz - bunların daha yararlı olduğunu düşünüyoruz.