[R-sig-ME] nested mixed effects logistic regression binomial glm) results differ by function.

[Lize van der Merwe]

Please advise:

I have a dichotomous outcome on 2500 individuals. From 18 geographical
areas, and many households nested within areas. I need to assess the
association between various predictors and my outcome, adjusting for the
correlation within households, as well as within areas. The following R
functions provide dramatically different results.

glmer(CC~predictor+1|area/household,family=binomial)

and

glmmPQL(CC~predictor, random=~1|area/household),family=binomial)

Why? Which is correct?

Thanks in advance.  (I posted this on another site too.)

Lize

 


	[[alternative HTML version deleted]]