Skip to content
Prev 307 / 7420 Next

glmer , quasipoisson and standard errors of the coefficients

Albert Romero wrote:
There are several issues here (and you should think about asking this 
question on r-sig-mixed-models , where there
is more expertise).

  1. glmer with a quasipoisson link does not provide a likelihood 
(rather, a quasilikelihood),
so you shouldn't necessarily assume that you can do *any* 
likelihood-based inference with
the results from this analysis.  The most conservative approach is to 
use only the estimated
standard errors or Z statistics on the parameters (this is a Wald test) 
for inference.

  2. if you _do_ want to use the likelihood for inference, you need to 
convert it to a "quasi-likelihood"
by dividing it by the estimated scale parameter (which you can extract 
via lme4:::sigma(model) ).
This is done (for example) when you state test="F" in anova() on 
generalized linear models --
also note that in this case you should technically do an F test rather 
than a chi-squared test,
with denominator df equal to the residual df (although it's hard to 
figure out what these should
be ...)

  3. it is generally advised *NOT* to use the likelihood ratio test for 
testing fixed effects
(see Pinheiro and Bates 2000).

  So ... bottom line ... for now, I would use the Z tests (labeled as t 
tests), i.e.
means / standard errors ...

  good luck,
   Ben Bolker