两组的 R 生存曲线看起来很奇怪
R Survival curve with two groups looks weird
我使用以下代码估计了两组的 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 条不同的生存曲线。有人知道我可能做错了什么吗?谢谢!!
最好包含一些样本数据,这样我们就能知道问题出在哪里。但是,让我们看看是否可以创建一些数据来重现您的问题:
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))
现在,使用您的代码,我们可以获得 similar-looking 结果:
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-09-19 创建
我使用以下代码估计了两组的 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 条不同的生存曲线。有人知道我可能做错了什么吗?谢谢!!
最好包含一些样本数据,这样我们就能知道问题出在哪里。但是,让我们看看是否可以创建一些数据来重现您的问题:
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))
现在,使用您的代码,我们可以获得 similar-looking 结果:
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)