卡普兰梅尔与 R

Question

我是 R 的新手，所以我对置信区间有一些问题。这是我的代码：

    library(survival)
    zeiten <- c(7, 3, 24, 19, 8, 7, 11, 19)
    delta <- c(1, 1, 1, 1, 1, 1, 0, 0)
    survivalobjekt <- Surv(zeiten, delta)
    km1 <- survfit(survivalobjekt ~ 1)
    plot(x=km1, col=c("black", "red", "red"), xlab="x", ylab="S(x)", lty=1)

所以现在我们检查 km1 的摘要：

    Call: survfit(formula = survivalobjekt ~ 1)

    time n.risk n.event survival std.err lower 95% CI upper 95% CI
    3      8       1    0.875   0.117        0.673        1.000
    7      7       2    0.625   0.171        0.365        1.000
    8      5       1    0.500   0.177        0.250        1.000
   19      3       1    0.333   0.180        0.116        0.961
   24      1       1    0.000     NaN           NA           NA

问题最终是这些 NA。当我使用情节（km1）时，我得到 , 但我想得到 我怎么才能得到它？我怎样才能删除那些NA？

Answer 1

您可以使用 ggplot2 绘制 Kaplan-Meier。
根据您的需要绘制置信区间。

library(ggfortify)
library(survival)
zeiten <- c(7, 3, 24, 19, 8, 7, 11, 19)
delta <- c(1, 1, 1, 1, 1, 1, 0, 0)
km1 <- survfit(Surv(zeiten, delta) ~ 1)
autoplot(km1, censor.shape = '*', censor.size = 5,
        surv.colour = 'darkorange', conf.int.fill = 'orange',
        surv.size = 1, censor.colour = 'red')