卡普兰梅尔与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')