我正在尝试创建我的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作为AM或PM，m5作为天气。
到目前为止我有这个但是后来卡住了似乎找不到任何有帮助的答案，因为当我使用“$”时它会产生我不会理解的错误。
m5<-update(m5,REML=TRUE)
par(mfrow=c(1,1))
mam.graph1<-update(m5,~.-1)
感谢