On Mar 19, 2011, at 8:44 AM, Juan Andres Hernandez wrote:
Dear R users, I use the excelent Anova function of the library car
because
the easy way to get sphericity correction. Unless I use the scan
function. I
have not been able to access the values of sum squares and
degrees of
freedom for each effect in the univariate summary table.
Example of the car library for Anova function:
library(car)
phase - factor(rep(c(pretest, posttest, followup), c(5, 5, 5)),
levels=c(pretest, posttest, followup))
hour - ordered(rep(1:5, 3))
idata - data.frame(phase, hour)
mod.ok - lm(cbind(pre.1, pre.2, pre.3, pre.4, pre.5,
post.1, post.2, post.3, post.4, post.5,
fup.1, fup.2, fup.3, fup.4, fup.5) ~
treatment*gender,
data=OBrienKaiser)
av.ok - Anova(mod.ok, idata=idata, idesign=~phase*hour)
summary(av.ok, multivariate=FALSE)
If you assign the above line to a new object you will get a
multivariate
summary and not the univariate one.
my.summary - summary(av.ok, multivariate=FALSE)
my.summary
Does anyone know how to assign the full univariate output to a new
object
or variable. The logical way does not produce the desired result.
Than'ks in advance.
summary,Anova.mlm appears to returns the object that was passed to it.
How to I know this (despite never having used the function before and
having tried unsuccessfully to access the answer with
help(summary.Anova.mlm))? I looked at the code. You are getting
console output that varies a bit with different parameters to summary,
but the object av.ok being passed as the value of summary and is not
changing.
all.equal(av.ok, my.summary)
[1] TRUE
If you want to capture the text output that appears on the console
display then you need to use either sink() or capture.output().
Juan Hernández Cabrera
Facultad de Psicología
Universidad de La Laguna
Spain
--
David Winsemius, MD
West Hartford, CT
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
