得到一个“冲浪”对象，它将是相同大小的原始数据，如果有联系？(生存，R)

Question

我要用Kaplan估计器估计我的数据中每个人的生存概率。survfit(Surv(.))函数通过降序计算每个唯一时间的生存概率。

从一个survfit 对象中获取每个个体的生存概率的优雅方法是什么？

library(survival)
data = data.frame(cbind(id = 1:10, 
            time = c(2,3,4,5,2,3,8,9,10,11), 
            status = c(1,0,0,1,0,1,0,1,0,1)))
survfit(Surv(data$time, data$status)~1)$surv
 # 0.9000 0.7875 0.7875 0.6300 0.6300 0.4200 0.4200 0.0000
survfit(Surv(data$time, data$status)~1)$time
 #  2  3  4  5  8  9 10 11

提前感谢！

Answer 1

您可以将结果放入data.frame中，然后将其合并回

merge(data, 
   with(survfit(Surv(data$time, data$status)~1), data.frame(time, surv))
)

在这里，我只需使用with()轻松地从survfit结果中提取两列。