1

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

兆瓦：

## 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)

我正在尝试针对每个模型的每个协变量绘制预测值。我找到了sjPlot库和plot_model函数，它们可以在使用type = "pred". 在每个模型上单独调用此函数效果很好，并且为每个模型生成两个单独的数字，如下所示：

在此处输入图像描述

但是我不熟悉 R，我很难在同一个图中绘制 4 个图。

该plot_model函数有一个grid参数，该参数仅适用于具有泊松分布的模型。对于gm1and gm2，我在调用时收到以下错误plot_model(gm1, type = "pred", grid = TRUE)：

Error in if (attr(x, "logistic", exact = TRUE) == "1" && attr(x, "is.trial",  : missing value where TRUE/FALSE needed

无论如何，我无法使用它在一个图中绘制三个模型，所以我尝试了三种不同的方法。首先，我看到了plot_models函数，它以多个模型作为输入。当我尝试将这两个模型作为参数传递时，调用plot_models(gm1, gm2)我得到以下错误：

Error: $ operator not defined for this S4 class

其次，我尝试使用par设置mfrow然后plot_model再次调用的功能没有成功。我没有收到任何错误，但这些图一直显示为单独的数字。

第三，我尝试使用该gridExtra库。打电话

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

导致以下错误：

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

有人对此有见解吗？

编辑

4

1 回答 1

3

这似乎有效：

pp1 <- plot_model(gm1,type="pred")
pp2 <- plot_model(gm2,type="pred")
plot_grid(c(pp1,pp2))
于 2019-05-04T01:15:07.080 回答