[R] Repeated aov residuals

[ripley@stats.ox.ac.uk]

On Tue, 2 Apr 2002, Yvonnick NOEL wrote:

> Are there any access functions to the various residual variables that should
> result from a repeated measures ANOVA ? MyAOVObject$residuals does not exist,
> and simply printing MyAOVObject gives a very long print of all fields in the
> result list, many of which I can't see what they are exactly :
> $error.qr$qraux, for instance.

Do this by stratum: see e.g. the MASS script which contains

oats.aov <- aov(Y ~ Nf*V + Error(B/V), data = oats, qr = T)
plot(fitted(oats.aov[[4]]), studres(oats.aov[[4]]))
abline(h=0, lty=2)
oats.pr <- proj(oats.aov)
qqnorm(oats.pr[[4]][,"Residuals"], ylab="Stratum 4 residuals")
qqline(oats.pr[[4]][,"Residuals"])

You may need the book to understand the context, but a multi-stratum aov
object is a list of aov objects, one for each stratum.

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._