Hi, I am trying to run a logistic regression to look at the risk of malaria infection in individuals. I want to account for intra household correlation and so want to include a household level random effect. I have been using the lmer command in lme4 package but am getting some strange results that are completely different to those I get using STATA.
Can I just check that this is the correct code lmer(IsPos ~ Dist + (1 | HouseID), family=binomial) Where IsPos is a binary vector of positives and negatives, Dist is the variable of interest and HouseID is the household ID number. For those STATA users, the equivalent I'm using is xtlogit IsPos Dist, i(HouseID) Any help would be much appreciated..! Thanks! Hugh -- View this message in context: http://r.789695.n4.nabble.com/random-effects-using-lmer-tp4566030p4566030.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
