[R] Simulate binary data for a logistic regression Monte Carlo
bolker at ufl.edu
Wed Apr 8 15:01:42 CEST 2009
Rolf Turner-3 wrote:
> On 8/04/2009, at 1:27 PM, Ben Bolker wrote:
>> I agree that that the individual-level random effect is probably
>> the issue.
>> I played with this some today but didn't manage to resolve it --
>> tried JAGS/R2jags and glmer from lme4 but didn't manage to
>> get an estimate of epsilon that matched the input value. I'm
>> a little worried about binary data with an underlying random
>> effect, I think there's an identifiability problem there ...
> Oh dear. This is getting hairier than I expected ... and a bit
> over my head. The toy model proposed *looks* simple enough at
> first blush. I would've thought that if a model which is that
> simple doesn't ``work'', then what hope is there for the more
> complicated models that one needs for analyzing real data?
I don't think this is impossible, but it's a bit subtle, and might
be worth taking up on r-sig-mixed-models instead ...
View this message in context: http://www.nabble.com/Simulate-binary-data-for-a-logistic-regression-Monte-Carlo-tp22928872p22949970.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help