Вычислить вероятности после coxph

Question

Я пытаюсь выяснить, как рассчитать некоторые вероятности после запуска анализа выживания с помощью coxph от survival пакет. Я прочитал кучу постов, но не могу найти ни одного, который отвечал бы на этот конкретный вопрос.

У меня есть модель как

coxph(Surv(time, event) ~ x1 + x2)  

И мне интересно вычислить вероятность того, что человек, доживший до времени t умрет раньше времени t+1,

Это возможно? Я думаю, что мне просто нужно найти оценку функции опасности в момент времени t, но мне интересно, как это вычислить. Спасибо за любой совет.

Answer 1

Да, это можно сделать с помощью следующего кода:

cox.fit<- coxph(Surv(df$time,df$Status)~X1+X2+X3+X4))
summary(cox.fit)
predicted_cox.fit<-survfit(cox.fit)
Exp_Beta<- predict(predicted_cox.fit,type='lp')
basehaz_cox<-basehaz(cox.fit)

Predicted_probability_at_t1 <- as.data.frame((exp(-basehaz_cox(t1))^(exp(Exp_Beta))