Используя объект survfit, я могу извлечь квантили следующим образом:
library('survival')
fit <- survfit(Surv(time, status) ~ 1, data = aml)
quantile(fit)[[1]]
#25 50 75
#12 27 43
Есть ли функция для выполнения обратного? Передать вектор квантилей и вернуть вероятности?
Вы можете построить предполагаемую функцию выживания с помощью stepfun:
survfun <- stepfun(fit$time, c(1, fit$surv))
Тогда survfun — это функция, которая возвращает вероятности выживания.
Другой способ — использовать аргумент times в summary:
summary(fit, times = c(1, 2, 3))
Это должно дать то же самое, что и survfun(c(1,2,3)).
