[R-meta] Two large forest plots in one image

Sun, Aug 25, 2024 1:40 PM

Dear Guido,

Apologies for the late reply.

I have been trying to work it out with your code, but unfortunately,
it was not coming out nicely.

However, I managed to use the {magik} package to work around it with
visually appealing images.

Kindest regards,

Emanuel



Below is the code:

Save each forest plot using the Plots -> Export -> Save as Image.

Save each forest plot in png format using a constant width of 950 and
an appropriate height so they fit well.



library(magick)



# Load the images

img1 <- image_read("image_1.png")

img2 <- image_read("image_2.png")

img3 <- image_read("image_3.png")



# Annotate each image with a letter

img1_annotated <- image_annotate(img1, "A", size = 40, gravity =
"northwest", color = "black")

img2_annotated <- image_annotate(img2, "B", size = 40, gravity =
"northwest", color = "black")

img3_annotated <- image_annotate(img3, "C", size = 40, gravity =
"northwest", color = "black")



# Display the annotated images (optional)

img1_annotated

img2_annotated

img3_annotated



# Combine images into a grid

combined_image <- image_append(c(img1_annotated,
img2_annotated,img3_annotated), stack = TRUE)



# add extra white space around the combined plot

padded_legend <- image_border(combined_image, color = "white",
geometry = "50x20")

padded_legend



# Save the combined image

image_write(padded_legend, path = "combined_plots_123.png")


On Mon, Jul 29, 2024 at 5:49?PM Dr. Guido Schwarzer

<guido.schwarzer at uniklinik-freiburg.de> wrote:

Emanuel,

I typically generate a PDF-file using pushViewport() from the grid package. The code below works for two forest plots either side by side or one below the other.

Best,
Guido


library("meta")
library("grid")

E1 <- structure(list(study = c("Barratt", "Carli", "Gupta", "Taqi",
"Chen", "H?bner", "Zutshi", "Park", "Klotz", "Mohamad", "Chen",
"Tsui"), year = c(2000, 2000, 2000, 2000, 2000, 2000, 2000, 2000,
2000, 2000, 2000, 2000), ev.EZ = c(1, 3, 1, 1, 0, 2, 0, 1, 1,
0, 0, 0), n.EZ = c(55, 41, 60, 9, 26, 65, 31, 25, 21, 28, 25,
57), ev.PPA = c(1, 4, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0), n.PPA = c(55,
42, 60, 11, 27, 57, 28, 22, 21, 28, 25, 54), TRX = c(1, 2, 1,
1, 1, 1, 1, 2, 2, 2, 2, 2), TIME = c("Update", "Update", "Update",
"Old", "Old", "Old", "Old", "Old", "Old", "Old", "Old", "Old"
)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA,
-12L))


E2 <- structure(list(study = c("Steinberg", "Parker", "Madej", "Schumann",
"Carli", "Tsui", "Xu", "Park", "Cho", "Xu", "Yosunkaya", "Mann",
"Jayr", "Aceto"), year = c(2000, 2000, 2000, 2000, 2000, 2000,
2000, 2000, 2000, 2000, 2000, 2000, 2000, 2000), ev.EZ = c(8,
3, 6, 16, 10, 10, 4, 11, 7, 5, 16, 30, 16, 2), n.EZ = c(51, 55,
41, 60, 20, 31, 20, 51, 20, 36, 21, 57, 38, 20), ev.PPA = c(9,
11, 8, 27, 11, 10, 3, 3, 9, 11, 17, 29, 14, 0), n.PPA = c(51,
55, 42, 60, 20, 33, 21, 48, 10, 32, 21, 54, 46, 19), TRX = c(2,
1, 1, 2, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2), TIME = c("Update", "Update",
"Update", "Update", "Old", "Old", "Old", "Old", "Old", "Old",
"Old", "Old", "Old", "Old")), class = c("tbl_df", "tbl", "data.frame"
), row.names = c(NA, -14L))


M3 <-  metabin(ev.EZ, n.EZ, ev.PPA, n.PPA,
               data = E1,
               studlab = paste0(study, " (", year, ")"),
               sm = "RR",
               method = "MH",
               subgroup = TIME)

M4 <- metabin(ev.EZ, n.EZ, ev.PPA, n.PPA,
              data = E2,
              studlab = paste0(study, " (", year, ")"),
              sm = "RR",
              method = "MH",
              subgroup = TRX)


height <- 14
width <- 9
#
ncols <- 1
nrows <- 2
#
firstrow <- 0.5
firstcol <- 0.5

pdf("two-forestplots.pdf", width = width, height = height)

pushViewport(
  viewport(
    layout = grid.layout(
    nrow = nrows, ncol = ncols,
    widths = unit(2 * c(firstcol, 1 - firstcol) * width / ncols,
                  "inches"),
    heights = unit(2 * c(firstrow, 1 - firstrow) * height / nrows,
                  "inches"))
  )
  )

pushViewport(viewport(layout.pos.col = 1, layout.pos.row = 1))

forest(M3, new = FALSE,
       common = FALSE,
       test.overall.random=TRUE,
       prediction=TRUE,
       print.Q=TRUE,
       print.Q.pval=TRUE,
       just = "right",
       layout = "RevMan5",
       pooled.events=TRUE,
       sortvar = year)

grid.text("Figure A", 0.025, 0.96, gp = grid::gpar(cex = 1.5), just = "left")

popViewport()

pushViewport(viewport(layout.pos.col = 1, layout.pos.row = 2))

forest(M4, new = FALSE,
       common = FALSE,
       test.overall.random=TRUE,
       prediction=TRUE,
       print.Q=TRUE,
       print.Q.pval=TRUE,
       just = "right",
       layout = "RevMan5",
       pooled.events=TRUE,
       sortvar = year)

grid.text("Figure B", 0.025, 0.98, gp = grid::gpar(cex = 1.5), just = "left")

popViewport()

dev.off()