我正在从我运行的许多GLMM中绘制一些模型预测,但是在分离因子水平以获得更合适的图时遇到了麻烦。这是我目前拥有的示例:
我正在使用sjPlot中的plot_model和ggplot2中的各种绘图修改器的组合来绘制我的模型,但是我用来解决此问题的每个函数都是空的。我得到的最接近的是使用facet_wrap将所有级别分离到单独的图中。另一个例子:
这是我当前正在运行的代码:
#topINBU is a glmer object
plot_model(topINBU, type = "int", title = "Indigo bunting",
legend.title = "Treatment", dot.size = 3, line.size = 1) +
labs(x = "", y = "Bird density (territories/ha)") +
theme(plot.title = element_text(hjust = 0.5, size = 15),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
panel.border = element_blank(),
axis.line = element_line(colour = "black"),
axis.title.x = element_text(margin = margin(t = 5, r = 0, b = 0, l = 0)),
axis.title.y = element_text(margin = margin(t = 0, r = 5, b = 0, l = 0)),
text = element_text(size = 15)) +
geom_line(size = 1, position = position_dodge(0.1))
我希望能够将它们分为两个图。其中一个显示了明确的趋势,不均衡的管理年龄和中间收获的趋势。另一个显示休假,没有收获,没有增长。任何帮助将不胜感激。