卡普兰·迈耶与R

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

最后的问题是这些。当我使用plot（km1）时，我得到 , 但是我想

---

我怎样才能得到它？我怎样才能去掉那些NA

您可以使用

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

库（ggfortify）
图书馆（生存）
zeiten您可以使用
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') 

库（ggfortify）
图书馆（生存）
zeiten您想要什么来代替NAs？只想要有置信区间直到结束您想要什么来代替NAs？只想要有置信区间直到结束谢谢Marco！但我有两个错误：Warnmeldungen:1：删除了1行，其中包含缺少的值（geom_point）。2：删除了1行包含缺失值（geom_文本）。我的绘图与您的绘图不同：（@H.Down我的代码运行时没有错误或警告，我的绘图与您想要的绘图完全相同！您可以发布您尝试运行的代码吗？您更改了我的代码，我怀疑您有一些错误。@H.Down我怀疑您更改了我的代码，并且您现在正在使用
survMisc
的
自动绘图
！是吗t正确？完成。您能告诉我，我如何在绘图中插入图例吗？例如：红线是观察数据；橙色空格是95%的置信区间。谢谢！谢谢Marco！但我有两个错误：Warnmeldungen:1：删除了1行包含缺失值（几何点）.2：删除了1行包含缺失值（geom_文本）。我的绘图与您的不同：（@H.Down我的代码运行时没有错误或警告，我的绘图与您想要的绘图完全相同！您可以发布您尝试运行的代码吗？您更改了我的代码，我怀疑您有一些错误。@H.Down我怀疑您更改了我的代码，并且您现在正在使用
survMisc
的
自动绘图
！是吗是吗？完成了。你能告诉我，我如何在绘图中插入图例吗？比如：红线是观察数据；橙色空格是95%的置信区间。谢谢！