Exactly what rule has been broken when lme (package nlme) produces the
error message Incompatible formulas for groups in random and
correlation? Here is what I am trying to do:
lme(Y ~ trait-1, random = ~trait-1|sire,
corr = corSymm(form = ~trait|cow))
Trait, sire and cow are factors (cow
Not convinced that responses so far have addressed the problem. The
model is
y = mu + U + e
where e is a vector of independendent errors with variance ve, and U
is a vector of random effects with covariance matrix va*A, where A is a
known matrix (which we can assume is a correlation matrix). If
