Grid.arrange объекты лесного фонда

Question

У меня есть пара лесных участков, сделанных с использованием пакета R forestplot. Я хотел бы grid.arrange их - но, похоже, это не работает легко.

Пример:

library(forestplot) 
# A basic example, create some fake data 
row_names <- list(list("test = 1", expression(test >= 2))) 
test_data <- data.frame(coef=c(1.59, 1.24), 
        low=c(1.4, 0.78), 
        high=c(1.8, 1.55)) 

forestplot(row_names, 
     test_data$coef, 
     test_data$low, 
     test_data$high, 
     zero = 1, 
     cex = 2, 
     lineheight = "auto", 
     xlab = "Lab axis txt") 

Хорошо это нарисует сюжет. Теперь предположим, что я хочу, чтобы захватить, что в объект сюжете бок о бок с другим сюжетом:

fp1 <- forestplot(row_names, 
     test_data$coef, 
     test_data$low, 
     test_data$high, 
     zero = 1, 
     cex = 2, 
     lineheight = "auto", 
     xlab = "Lab axis txt") 

Следующая выдает ошибку:

> grid.arrange(fp1, fp1) 
Hit <Return> to see next plot: 
Error in gList(list(path = "GRID.VP.7537", name = "arrange.1-1-1-1", n = 2L, : 
only 'grobs' allowed in "gList" 

Так Клири FP1 не является Grob - но как Я достигаю этого другими способами?

Answer 1

См. Второй пример на страницах справки ?forestplot. который показывает, как это сделать.

---

forestplot, кажется, не вернуть сюжет: смотреть на str(fp1).

Пара вариантов: grid, чтобы создать пространство для проецирования (v1), или захватить график, а затем объединить (v2).

---

v1 с использованием сетки

library(grid) 
library(forestplot) 

# Create 2 rows by one columns viewport 
grid.newpage() 
pushViewport(viewport(layout=grid.layout(2, 1))) 

# Plot in viewport position 1x1 
pushViewport(viewport(layout.pos.row=1, layout.pos.col=1)) 
forestplot(row_names, 
     test_data$coef, 
     test_data$low, 
     test_data$high, 
     zero = 1, 
     cex = 2, 
     lineheight = "auto", 
     xlab = "Lab axis txt") 
upViewport(1) 

# Plot in viewport position 2x1 
pushViewport(viewport(layout.pos.row=2, layout.pos.col=1)) 
forestplot(row_names, 
     test_data$coef, 
     test_data$low, 
     test_data$high, 
     zero = 1, 
     cex = 2, 
     lineheight = "auto", 
     xlab = "Lab axis txt", 
     new_page = FALSE) 
upViewport(1) 

---

v2, захватить сюжет

fp1 <- grid.grabExpr(print(forestplot(row_names, 
     test_data$coef, 
     test_data$low, 
     test_data$high, 
     zero = 1, 
     cex = 2, 
     lineheight = "auto", 
     xlab = "Lab axis txt"))) 

gridExtra::grid.arrange(fp1, fp1)