`levels<-`(`*tmp*`, value = as.character(levels)) 中的錯誤：因子級別 [3] 重復

Question

我正在嘗試使用自動繪圖 plot 刻面生存曲線，但協變量和刻面的組合會在因子內重復水平

library(survival)
library(ggfortify)
fit <- survfit( Surv(time, status) ~ inst + sex,
                 data = lung )

autoplot(fit, facets = TRUE)


Error in `levels<-`(`*tmp*`, value = as.character(levels)) : 
  factor level [3] is duplicated

有沒有人用自動繪圖成功地繪制了多面生存曲線？ 我嘗試了 survminer，但 plot 看起來很可怕，協變量占據了 plot 區域的大部分。

Answer 1

我認為您應該再次查看ggsurvplot ，因為autoplot.survfit似乎不喜歡擁有多個自變量（無論您是否考慮）。

ggsurvplot function 返回一個 ggplot object，因此您無需滿足於默認選項。 您可以根據需要添加比例和樣式。 以您為例，我們可以這樣做：

library(survival)
library(ggfortify)
library(survminer)

fit <- survfit( Surv(time, status) ~ inst + sex,
                 data = lung )

p <- ggsurvplot(fit, facet.by = "inst", conf.int = TRUE) + 
  theme(strip.background = element_blank(),
        axis.line.x = element_line())

p$facet <- facet_wrap(.~inst, ncol = 3, nrow = 6, scales = "free")

p