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