如何将多个glmer模型绘制到一个绘图中？

我有两个

glmer

模型，每个模型有两个协变量，我试图将它们绘制成一个图

MWE：

## generalized linear mixed model
library(lattice)

cbpp$response <- sample(c(0,1), replace=TRUE, size=nrow(cbpp))
gm1 <- glmer(response ~ size + incidence + (1 | herd),
              data = cbpp, family = binomial)

cbpp$obs <- 1:nrow(cbpp)
gm2 <- glmer(response ~ size + incidence + (1 | herd) +  (1|obs),
              family = binomial, data = cbpp)

无论如何，我不能用这个在一个图形中绘制三个模型，所以我尝试了三种不同的方法。首先，我看到了函数，它将多个模型作为输入。当我尝试将这两个模型作为参数传递时，调用

plot\u models（gm1，gm2）

我得到以下错误：

Error: $ operator not defined for this S4 class

Error in gList(list(ppt = list(data = list(x = c(-2, -1, 0, 1, 2, 3, 4,  : only 'grobs' allowed in "gList"

其次，我尝试使用

par

函数设置

mfrow

，然后再次调用

plot\u model

，但没有成功。我没有得到任何错误，但图一直显示为单独的数字

第三，我尝试使用

gridExtra

库。召唤

p1 <- plot_model(gm1, type = "pred")
p2 <- plot_model(gm2, type = "pred")
grid.arrange(p1, p2)

有人对此有见解吗

编辑

这似乎有效：

pp1 <- plot_model(gm1,type="pred")
pp2 <- plot_model(gm2,type="pred")
plot_grid(c(pp1,pp2))

---

pp1请编辑您的帖子，包括可复制的样本数据和代码，为我们提供一些有用的东西。在不了解导致错误的数据的情况下，很难（即不可能）帮助解决“我遇到错误XYZ”形式的问题。我建议看一下如何提供一个。这里只是一个注释，适用于我的语法是“plot_grid（list（pp1，pp2））”
pp1 <- plot_model(gm1,type="pred")
pp2 <- plot_model(gm2,type="pred")
plot_grid(c(pp1,pp2))