我无法打印幸存者的平均值(...)

时间:2010-01-24 05:57:34

标签: r

我在R中的生存包中使用了数据,aml并计算了生存函数 “survfit”。由于survfit的结果没有显示平均值, 我使用以下代码打印平均值:

> print(survfit(Surv(aml1$time,aml1$status)~1),show.rmean=T)

(我使用的数据是aml1 <-aml[aml$x="Maintained",]

上面的代码在我的朋友电脑上工作,但不是我的。 所以,我考虑下载一些额外的包来使用print(...show.rmean=T)。 但是,“打印”是基本的,所以我不需要任何东西来运行打印。 那为什么我不能得到平均值?

1 个答案:

答案 0 :(得分:3)

print是一个通用函数,其选项可能特定于它正在打印的东西。在这种情况下,选项由print.survfit提供。虽然我的survival版本实际上没有show.rmean。也许你想要print.rmean?有关详细信息,请参阅print.survfit的帮助。

> print(survfit(Surv(aml1$time,aml1$status)~1), print.rmean=T)
Call: survfit(formula = Surv(aml1$time, aml1$status) ~ 1)

   records      n.max    n.start     events     *rmean *se(rmean)     median 
     23.00      23.00      23.00      18.00      36.36       9.85      27.00 
   0.95LCL    0.95UCL 
     18.00      45.00 
    * restricted mean with upper limit =  161