2 つのカテゴリ変数によって説明される MAM のグラフィカル表現を作成しようとしています。2つの間の相互作用は重要ではないため、各説明変数の効果を示すために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)
2 つのグラフをプロットするには、午前または午後には m4 を使用し、天気には m5 を使用する必要があると思います。
私はこれまでのところこれを持っていますが、行き詰まり、「$」を使用すると理解できないエラーが発生するため、役立つ答えがどこにも見つからないようです。
m5<-update(m5,REML=TRUE)
par(mfrow=c(1,1))
mam.graph1<-update(m5,~.-1)
ありがとう