Skip to content
Prev 5797 / 20628 Next

Mixed effect logistic regression help

CJ Griffiths <Christine.Griffiths at ...> writes:
...

A little late in the response...

To my understanding, mixed models cannot accommodate the paired 1/0 structure of
your data.  Adding micro as a mixed effect is only accounting for the variance
in your 1/0 response among different levels of micro.  However, since you have
an equal number of 1s and 0s for each level, you see zero as the variance in
micro.  You need to use a conditional logistic model (e.g., clogit in package
survival) to do what you want to do to accommodate the paired case/control
nature of your samples.  However, clogit cannot accommodate the repeated
sampling of individuals (which lmer, nlme, etc. can, but not in this case).  The
only real application I've seen bridging these two 'worlds' is the following:

Duchesne, T., D. Fortin, and N. Courbin. 2010. Mixed conditional logistic
regression for habitat selection studies. Journal of Animal Ecology 79:548-555.

Unfortunately, it was implemented in Matlab, so not so useful to any of us in R
land.  Personally, I'd be very excited to see a mixed conditional logit package
or extension developed for R, but I'm neither statistician nor programmer.  

Giancarlo