What you should and can do depends on the expectations you have regarding the correlation structure of the data and is limited by your degrees of freedom. For example, what is wrong about:
reg<-lmer(response~femaleset+treatment+(1|group)) ? This assumes that there are constant group effects on a group's responses (i.e., it allows for differences in mean responses between groups) and that this constant group effect is independent of the effect of femaleset and treatment on a group. Best, Daniel -- View this message in context: http://r.789695.n4.nabble.com/lmer-with-2-random-effects-with-only-two-levels-tp3536791p3537184.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.