9

我创建了这个模型：

model <- survfit(Surv(time,status)~c$sex)
model

输出是：

Call: survfit(formula = Surv(time, status) ~ c$sex)

             records n.max n.start events median 0.95LCL 0.95UCL
c$sex=female      15    15      15      8    720     517      NA
c$sex=male        28    28      28     23    234     145     712

所以，我想提取男性的中位数和女性的中位数，但不知道该怎么做。这是我的尝试：

>model$median
NULL

>summary(model)$table[, "median"]
c$sex=female c$sex=male 
       720.0            234.5

我只想要每个值（“720”和“234.5”），有人可以帮助我吗？

提前致谢

4

2 回答 2

12

你已经得到它了。您在屏幕上看到的只是names长度为 2 的向量的属性。

fit <- survfit(Surv(time, status) ~ x, data = aml)
summary(fit)$table
#                records n.max n.start events median 0.95LCL 0.95UCL
#x=Maintained         11    11      11      7     31      18      NA
#x=Nonmaintained      12    12      12     11     23       8      NA

#  Access first value like any other vector
summary(fit)$table[,'median'][1]
#x=Maintained 
#          31

要打印没有名称，请使用unname()...

unname(summary(fit)$table[,'median'])
# [1] 31 23

但是你不需要unname()他们把它们当作一个数值来使用，那只是一种审美欲望......

sum( summary(fit)$table[,'median'] )
[1] 54

为了进一步证明（！）它是一个向量，请使用str()...

str(summary(fit)$table[,'median'])
# Named num [1:2] 31 23
# - attr(*, "names")= chr [1:2] "x=Maintained" "x=Nonmaintained"
于 2013-10-23T17:11:42.417 回答
5

这也有效：

> library(survMisc)
> fit <- survfit(Surv(time, status) ~ x, data = aml)
> median(fit)
                median
x=Maintained        31
x=Nonmaintained     23

并且没有名称（即删除 a 的结构data.frame）：

> unname(unlist(median(fit)))
[1] 31 23

如果您还想要置信区间（默认为“log”），那就太好了：

> median(fit, CI=TRUE)
                median lower upper
x=Maintained        31    13    NA
x=Nonmaintained     23     5    43
于 2013-10-25T03:18:41.853 回答