[R] model comparison with mixed effects glm

[hadley wickham]

> ### To get around that, I computed 2*log(likelihood ratio) manually:
>
> lglk0 <- logLik(fit0)
> lglk.ID1. <- logLik(Fit.ID1.)
> pchisq(as.numeric(chisq.ID.), 1, lower=FALSE)
>  > [1] 0.008545848

(I think you're missing a line in there)

But isn't this rather perilous unless you are confident that the two
models are using exactly the same formulation of the likelihood?  (ie.
that they are truly nested)

Hadley