Я запустил модель GLMM, используя функцию glmer
из пакета lme4
. Я хотел бы знать значения р фиксированных эффектов. Поэтому я попробовал функцию mixed
из пакета afex
. Но результаты очень странные для меня, потому что кадры данных в результатах mixed
отличаются от кадров в anova
.
mixed(glmm.c,data=dt2, method="LRT")
...
Df full model: 16
Effect df Chisq p.value
1 year 1 0.98 .32
2 waters 3 1.10 .78
3 round 1 5.57 * .02
4 year:waters 3 1.10 .78
5 waters:round 3 1.27 .74
6 year:round 1 5.56 * .02
anova(glmm.c)
Analysis of Variance Table
Df Sum Sq Mean Sq F value
year 3 46.856 15.6187 15.6187
waters 3 12.478 4.1595 4.1595
round 5 4.613 0.9226 0.9226
year:waters 9 24.097 2.6775 2.6775
waters:round 15 27.808 1.8538 1.8538
year:round 15 27.659 1.8439 1.8439
Может ли быть кто-то, кто мог бы сказать мне причину?