我正在尝试创建由两个分类变量解释的 MAM 的图形表示。我意识到我需要制作 2 个箱线图来显示每个解释变量的效果，因为两者之间的交互作用微不足道。我如何绘制箱线图，因为我似乎无法提取手段等，而且我不确定它是否是因为我使用的是 lmer。
到目前为止，这是我的代码....
m1<-lmer(Bird.species~Weather*AM.or.PM+(1|Staff),family=poisson)
summary(m1)
m2<-update(m1,REML=FALSE)
m3<-update(m2,~.-Weather:AM.or.PM)
anova(m2,m3)
m4<-update(m3,~.-Weather)
anova(m3,m4)
m5<-update(m3,~.-AM.or.PM)
anova(m3,m5)
m6<-update(m5,~.-Weather)
anova(m5,m6)
m7<-update(m4,~.-AM.or.PM)
anova(m4,m6)
mam.lmer<-update(m3,REML=TRUE)
summary(mam.lmer)
我认为为了绘制这两个图表，我必须将 m4 用于上午或下午，将 m5 用于天气。
到目前为止我有这个但是后来卡住了，似乎无法在任何有帮助的地方找到答案，因为当我使用“$”时它会产生我不理解的错误。
m5<-update(m5,REML=TRUE)
par(mfrow=c(1,1))
mam.graph1<-update(m5,~.-1)
谢谢