Hello,

I run the following code.
1. REML fit was not tin the output
2. There was no residual in the output
3. I could not run anova(g)

> (g=lmer(total_count ~ c+(1|c:b:a), d2, REML=TRUE, family = "poisson"))
Generalized linear mixed model fit by the Laplace approximation
Formula: total_count ~ c + (1 | c:b:a)
   Data: d2
   AIC   BIC logLik deviance
 661.9 673.1 -326.9    653.9
Random effects:
 Groups Name        Variance Std.Dev.
 c:b:a  (Intercept) 1.8368   1.3553
Number of obs: 122, groups: c:b:a, 18
Fixed effects:
              Estimate Std. Error z value Pr(>|z|)
(Intercept)     1.6383     0.5605   2.923  0.00347 **
cLovers Point   0.2080     0.7924   0.262  0.79295
cPoint Pinos   -0.4282     0.7998  -0.535  0.59242
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Correlation of Fixed Effects:
            (Intr) cLvrsP
cLoversPont -0.707
cPointPinos -0.701  0.496
> anova(g)
Error in anova(g) : single argument anova for GLMMs not yet implemented
What might be wrong?
Thanks,
Y

	[[alternative HTML version deleted]]

