Kaplan Meier與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

您可以繪製Kaplan-Meier usin ggplot2。
置信區間根據您的需求繪製。

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