survival包里的方法quantile.survfit有反函数吗?

Is there an inverse function for the method quantile.survfit in the survival package?

使用 survfit 对象,我可以像这样提取分位数:

library('survival')
fit <- survfit(Surv(time, status) ~ 1, data = aml) 
quantile(fit)[[1]]
#25 50 75 
#12 27 43

有反函数吗?传递一个分位数向量,然后 return 概率?

你可以用stepfun构造估计的生存函数:

survfun <- stepfun(fit$time, c(1, fit$surv))

那么survfun是一个returns生存概率的函数。

另一种方法是在summary中使用times参数:

summary(fit, times = c(1, 2, 3))

这应该与 survfun(c(1,2,3)) 相同。