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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._