[R] model comparison with mixed effects glm

hadley wickham h.wickham at gmail.com
Tue Apr 4 19:20:25 CEST 2006

> ### 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)


More information about the R-help mailing list