Dear Agnes, I think your model specification should look like this:
YourModel1 <- lmerlmer(y ~ poptype*matingtype + (1|poptype:pop) + (1|poptype:fam), data = ...) The "1" in front of "|" refers to models that are random intercepts models as opposed to general random coefficients models in which case "1" would need to be replaced by a variable that is quantitative. So, the "(poptype|/pop/fam)" construction is definitely not relevant to your problem, unless "poptype" is a quantitative variable... The combined factor "poptype:pop" corresponds to the factor pop, but taking the nesting structure into account. Similar for the construction "poptype:fam". >From the summary output of the lmer() fit you should be able to check whether >or not the correct number of groups are used for these random factors. Kind regards Christian ______________________________________________ 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.
