两组的 R 生存曲线看起来很奇怪

Question

我使用此代码估计了具有两组的 Kaplan-Meier：

test3 <- survfit(Surv(wave, con3) ~ sex, data = consult3)
ggsurvplot(test3, color = "#2E9FDF",
       risk.table = TRUE, risk.table.y.text.col = TRUE)

尽管它产生了一个 KM，但它看起来很奇怪，并且没有向我显示 sex==1 或 sex==2 的 2 条不同的生存曲线。有人知道我可能做错了什么吗？谢谢！！

Answer 1

最好包含一些示例数据，这样我们就可以找出问题所在。但是，让我们看看我们是否可以创建一些数据来复制您的问题：

library(survival)
library(survminer)

set.seed(69)

consult3 <- data.frame(sex = rep(1:2, each = 50),
                       con3 = c(rbinom(50, 1, 0.2), rbinom(50, 1, 0.4)),
                       wave = sample(3, 100, TRUE))

现在，使用您的代码，我们可以获得类似的结果：

test3 <- survfit(Surv(wave, con3) ~ sex, data = consult3)

ggsurvplot(test3, risk.table = TRUE, color = "#2E9FDF",
           risk.table.y.text.col = TRUE)

据我所知，问题在于您设置了单一的颜色美学。解决方案是删除它：

ggsurvplot(test3, risk.table = TRUE, risk.table.y.text.col = TRUE)

如果您想控制线条的颜色，请使用palette 而不是color：

ggsurvplot(test3, palette = c("red", "forestgreen"), alpha = 0.5,
           risk.table = TRUE, risk.table.y.text.col = TRUE)

^{由reprex package (v0.3.0) 于 2020 年 9 月 19 日创建}