2010-01-24 61 views
1

我使用了數據,aml在R中的生存包中,並通過 「survfit」計算了生存函數。由於survfit的結果不顯示平均值, 我用下面的代碼打印的意思是:我無法打印survfit的平均值(...)

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

(在我使用的數據是,aml1 <-aml[aml$x="Maintained",]

上面的代碼在我的工作朋友電腦,但不是我的。 所以,我想下載一些額外的包使用print(...show.rmean=T)。 但是,「打印」是基本的,所以我不需要任何東西來運行打印。 那麼爲什麼我不能得到平均值?

+0

更新您的包存活...它爲我工作! – 2012-10-30 16:21:44

回答

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