在单一地块混合模型中绘制多个随机效应

Question

library("lme4")
data(sleepstudy)             

fit1 <- lmer(Reaction ~ Days + (1|Subject), sleepstudy) 

想象随机效果，

library(sjPlot)
plot_model(fit1,type = "re",facet.grid = FALSE) 

在我的原始数据中，我有三个随机组。但是，如果我想绘制随机效果，它们都有三个不同的图。如何将它们放在1 X 3面板或3 X 1面板中的所有单个图中。

Answer 1

您可以使用gridExtra::grid.arrange()

fit1 <- lmer(Reaction ~ (1|Days) + (1|Subject), sleepstudy) 

library(sjPlot)
p <- plot_model(fit1, type = "re", facet.grid=FALSE) 

library(gridExtra)
grid.arrange(p[[1]], p[[2]])

产品：

您还可以考虑使用lattice::qqmath()。

library(lattice)
p2 <- qqmath(ranef(fit1, condVar=TRUE))
grid.arrange(p2[[1]], p2[[2]])

产品：

注意： 要指定列，请使用ncol选项。比较例如grid.arrange(p2[[1]], p2[[2]], ncol=2)与grid.arrange(p2[[1]], p2[[2]], ncol=1)。