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))
相同。
使用 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))
相同。