我们可以在 R 森林图中用不同的颜色给不同的 rows/covariates/studies 上色吗？

Question

使用编程语言 R 中的 forestplot 包，我想制作一个每行颜色不同的森林图。每行，我指的是框和各自的置信区间。

以小插图 [https://cran.r-project.org/web/packages/forestplot/vignettes/forestplot.html] 为例，

library(forestplot)
# Cochrane data from the 'rmeta'-package
cochrane_from_rmeta <- 
  structure(list(
    mean  = c(NA, NA, 0.578, 0.165, 0.246, 0.700, 0.348, 0.139, 1.017, NA, 0.531), 
    lower = c(NA, NA, 0.372, 0.018, 0.072, 0.333, 0.083, 0.016, 0.365, NA, 0.386),
    upper = c(NA, NA, 0.898, 1.517, 0.833, 1.474, 1.455, 1.209, 2.831, NA, 0.731)),
    .Names = c("mean", "lower", "upper"), 
    row.names = c(NA, -11L), 
    class = "data.frame")

tabletext <- cbind(c("", "Study", "Auckland", "Block", "Doran", "Gamsu", "Morrison", "Papageorgiou", "Tauesch", NA, "Summary"),
  c("Deaths", "(steroid)", "36", "1", "4", "14", "3", "1", "8", NA, NA),
  c("Deaths", "(placebo)", "60", "5", "11", "20", "7", "7", "10", NA, NA),
  c("", "OR", "0.58", "0.16", "0.25", "0.70", "0.35", "0.14", "1.02", NA, "0.53"))

forestplot(tabletext, cochrane_from_rmeta, new_page = TRUE,
           is.summary = c(TRUE, TRUE, rep(FALSE, 8), TRUE),
           clip = c(0.1,2.5), xlog = TRUE, 
           col = fpColors(box="royalblue",line="darkblue", summary="royalblue"))

我希望每个书房都有自己的颜色（奥克兰可以是蓝色，布洛克可以是红色，多兰可以是绿色，等等）。我认为这可以通过更改 fpColors() 函数的参数来实现。

有什么办法吗？

Answer 1

看看fpShapesGp。有了这个，就可以用不同的颜色为行着色。

一个简单的代码示例：

styles <- fpShapesGp(
  lines = list(
    gpar(col = "black"),
    gpar(col = "blue"),
    gpar(col = "black"),
    gpar(col = "blue"),
    gpar(col = "black"),
    gpar(col = "blue"),
    gpar(col = "black"),
    gpar(col = "blue"),
    gpar(col = "black"),
    gpar(col = "blue"),
    gpar(col = "black")
  ),
  box = list(
    gpar(fill = "black"),
    gpar(fill = "blue"),
    gpar(fill = "black"),
    gpar(fill = "blue"),
    gpar(fill = "black"),
    gpar(fill = "blue"),
    gpar(fill = "black"),
    gpar(fill = "blue"),
    gpar(fill = "black"),
    gpar(fill = "blue"),
    gpar(fill = "black")
  ) 
)

forestplot(tabletext, cochrane_from_rmeta, new_page = TRUE,
           is.summary = c(TRUE, TRUE, rep(FALSE, 8), TRUE),
           clip = c(0.1,2.5), xlog = TRUE, 
           shapes_gp = styles)

结果你得到这个：

我们可以在 R 森林图中用不同的颜色给不同的 rows/covariates/studies 上色吗？

Can we color the different rows/covariates/studies in different colors in R forest plots?

r

survival-analysis

forestplot