2

のパッケージを使用して、survivalこのRような I Kaplan-Meiser-object を生成しました

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

print(fit)I getを使用して

> 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)

しかし、非常に大きなデータ セットを使用すると、サバイバル オブジェクトを複数回再作成するのは好ましくありません。

4

1 に答える 1

2

使用できますgrep：

containsText <- "sex=1"
fit_sex1 <- fit[c(grep(containsText ,names(fit$strata)))]
print(fit_sex1)
于 2013-04-02T09:34:22.540 に答える