[R] is there a way to use "aov" to do mixed linear models with both random and fixed effects?

Dieter Menne dieter.menne at menne-biomed.de
Mon Nov 10 08:56:46 CET 2008


Michael <comtech.usa <at> gmail.com> writes:

> There is no way to signal to "aov" the A and B are random effects and
> C is fixed effect; or A is random and B and C are fixed?

It is possible to do blocking with aov and the Error() term, but I would
recommend against it, even if this method is still taught in universities. Use
mixed models in package nlme, function lme, and you get a good documentation in
the book by Pinheiro/Bates.

Dieter



More information about the R-help mailing list