Question

在survival中使用包R我生成了像这样的I Kaplan-Meiser对象

library(survival)
fit <- survfit(Surv(time, status) ~ ph.ecog + sex, data=lung)

使用print(fit)我

> print(fit)
Call: survfit(formula = Surv(time, status) ~ ph.ecog + sex, data = lung)

   1 observation deleted due to missingness 
                 records n.max n.start events median 0.95LCL 0.95UCL
ph.ecog=0, sex=1      36    36      36     28    353     303     558
ph.ecog=0, sex=2      27    27      27      9    705     350      NA
ph.ecog=1, sex=1      71    71      71     54    239     207     363
ph.ecog=1, sex=2      42    42      42     28    450     345     687
ph.ecog=2, sex=1      29    29      29     28    166     105     288
ph.ecog=2, sex=2      21    21      21     16    239     199     444
ph.ecog=3, sex=1       1     1       1      1    118      NA      NA

我现在如何才能访问这个整体生存对象的一部分？比如仅为sex=1打印结果？

当然我可以做点什么

fit_sex1 <- survfit(Surv(time, status) ~ ph.ecog + sex, data=lung[lung$sex == 1,])
print(fit_sex1)

但是使用非常大的数据集，多次重建生存对象是不利的。

Answer 1

您可以使用grep：

containsText <- "sex=1"
fit_sex1 <- fit[c(grep(containsText ,names(fit$strata)))]
print(fit_sex1)

只访问生存对象的一部分

1 个答案: