Michael,
If you mean Manova in the car package (or manova in stats), just use
str() on the resulting object to see the slots. However, Manova()
doesn't return F values, since these depend on the type of test
statistic used for the multivariate test.
> library(car)
> ?Manova
> soils.mod <- lm(cbind(pH,N,Dens,P,Ca,Mg,K,Na,Conduc) ~ Block +
Contour*Depth,
+ data=Soils)
> man <-Manova(soils.mod)
> str(man)
List of 8
$ SSP :List of 4
..$ Block : num [1:9, 1:9] 1.2325 -0.0466 0.2716 -62.6558
-0.6192 ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
..$ Contour : num [1:9, 1:9] 0.2607 -0.0172 0.1002 23.9644
-1.5842 ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
..$ Depth : num [1:9, 1:9] 14.96 1.45 -4.30 1703.12
74.47 ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
..$ Contour:Depth: num [1:9, 1:9] 0.5159 0.0137 0.0457 -30.9419
1.5070 ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
.. .. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
$ SSPE : num [1:9, 1:9] 4.247 -0.049 -0.201 -118.777
9.368 ...
..- attr(*, "dimnames")=List of 2
.. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
.. ..$ : chr [1:9] "pH" "N" "Dens" "P" ...
$ df : Named num [1:4] 3 2 3 6
..- attr(*, "names")= chr [1:4] "Block" "Contour" "Depth" "Contour:Depth"
$ error.df: int 33
$ terms : chr [1:4] "Block" "Contour" "Depth" "Contour:Depth"
$ repeated: logi FALSE
$ type : chr "II"
$ test : chr "Pillai"
- attr(*, "class")= chr "Anova.mlm"
Michael Rennie wrote:
Hi there,
Does anyone know how to extract elements from the table returned by Manova()?
Using the univariate equivalent, Anova(), it's easy:
a.an<-Anova(lm(y~x1*x2))
a.an$F
This will return a vector of the F-values in order of the terms of the model.
However, a similar application using Manova():
m.an<-Manova(lm(Y~x1~x2))
m.an$F
Returns NULL. So does any attempt at calling the headers in the MANOVA table.
Any thoughts?
Mike
--
Michael Friendly Email: friendly AT yorku DOT ca
Professor, Psychology Dept.
York University Voice: 416 736-5115 x66249 Fax: 416 736-5814
4700 Keele Street http://www.math.yorku.ca/SCS/friendly.html
Toronto, ONT M3J 1P3 CANADA
______________________________________________
[email protected] 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.
