Re: [R] lmer: random factor nested in a fixed factor

2008-10-07 Thread Hank Stevens
hi folks, I believe that the construction y ~ A + (A|B) (where A is a categorical variable for a fixed effect and B is a categorical variable for a random effect) will generate random effects for all of the fixed effect coefficients. That is, in addition fitted a systematic fixed effect coef

Re: [R] lmer: random factor nested in a fixed factor

2008-10-06 Thread Christian Ritz
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 whi

[R] lmer: random factor nested in a fixed factor

2008-10-06 Thread Marie-Agnes Coutellec
Hi all, I try to build a model using lmer, with 2 crossed fixed factors (poptype and matingtype) and 2 random factors (pop and family) which I want to nest within poptype. Which of these formulae should I use ? 1. lmer(y~poptype*matingtype + (1|poptype/pop/fam)) 2. lmer(y~poptype*matingtype +