[R] Random Effects for One-Way Anova

[Peter Dalgaard]

Prof Brian Ripley <ripley at stats.ox.ac.uk> writes:

> > My questions is:
> > which formula should I use if I want to use a random effects model?
> > I think I should use "lme", but I don't know how.
> 
> aov() is the function for anova modelling. so you want
> 
> aov(A ~ B, data = DATA)         # fixed effects
> aov(A ~ Error(B), data = DATA)  # random effects
> 
> but a 1-way ANOVA with fixed or random effects is the same analysis, and 
> only the interpretation differs.

Not quite. Only if the design is balanced and you're not looking to
compute things like the variance of the estimated overall mean. 

In the unbalanced case, I don't believe aov() is doing the right thing
(it might not be very wrong, though, if the imbalance is slight) and
lme/lmer is closer on the mark.

-- 
   O__  ---- Peter Dalgaard             Øster Farimagsgade 5, Entr.B
  c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
 (*) \(*) -- University of Copenhagen   Denmark          Ph:  (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45) 35327907