我在R中的生存包中使用了数据,aml并计算了生存函数 “survfit”。由于survfit的结果没有显示平均值, 我使用以下代码打印平均值:
> print(survfit(Surv(aml1$time,aml1$status)~1),show.rmean=T)
(我使用的数据是aml1 <-aml[aml$x="Maintained",]
)
上面的代码在我的朋友电脑上工作,但不是我的。
所以,我考虑下载一些额外的包来使用print(...show.rmean=T)
。
但是,“打印”是基本的,所以我不需要任何东西来运行打印。
那为什么我不能得到平均值?
答案 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