[R-sig-ME] Non normal random effects
bates at stat.wisc.edu
Mon Nov 29 20:12:31 CET 2010
On Mon, Nov 29, 2010 at 9:02 AM, Gustavo Betini <betinig at uoguelph.ca> wrote:
> Lets say that I want to look at the interaction between Species and some other variable...Temperature. The right specification would be
> (Temperature|Species) + (Temperature|Station:Species)
> or simply
> (1|Species) + (Temperature|Station:Species)
That depends on whether you feel that the changes in temperature
coefficient would be associated with the Species or with the
Station:Species combination or both. I would try to start simply,
(1 + Temperature|Species) + (1|Station:Species)
It is difficult to estimate a large number of variance components
without a lot of data. Starting with the "sumo" model that has every
possible term and interaction then seeing if you can pare it back
makes more sense in fixed-effects models than in mixed-effects models.
Deletion strategies are easier to describe and evaluate when working
with fixed-effects only. Also, in the fixed-effects case having a
model with too many terms in it doesn't usually prevent estimation of
coefficients, it just make the estimators less precise.
In a mixed-effects model, adding interaction terms increases the
number of variance-component parameters quadratically, not linearly.
And it is a lot more difficult to determine that the
variance-covariance matrix of the random effects is singular than to
detect that a model matrix is rank deficient.
So I would recommend building mixed-effects models using forward
selection, not backward selection, and proceeding cautiously.
> Gustavo S. Betini
> Dept. of Integrative Biology
> University of Guelph, Canada
> On 10-11-27 12:04 PM, Douglas Bates wrote:
> Just in terms of the model specification, do you really mean
> (1|Species/Station)? That expands to
> (1|Species) + (1|Species:Station)
> and wouldn't reduce to (1|Station) in a model specification. I think
> you meant "species within station", which would be written as
> (1|Station/Species) although I prefer the more explicit form
> (1|Station) + (1|Station:Species)
> R-sig-mixed-models at r-project.org mailing list
More information about the R-sig-mixed-models