On May 3, 2013, at 14:59 , Ozgul Inceoglu wrote: > Dear All, I am trying to perform MANOVA. I have table with 504 > columns(species) and 36 rows) with two grouping (season and location) > > Zx <- Z[c(4:504)] > Zxm <- as.matrix(Z) > m<- manova(Zxm~Season*location, data=Z) > > when I do summary.aov, I get respond for each species but summary.manova > summary.manova(m) :" residuals have rank" 24<501. > > What can it be the reason for this error message?
Too many columns and too few rows. Multivariate tests require more degrees of freedom than response variables. > > Thank you, > > Ozgul > > Below you can see part of the table. > name Season location Acetobacter Aerococcus Alishewanella > Amaricoccus > xls-nord-01 J w 0 0,024078979 0 0 > bxls-sud-01 J w 0 0 0 0 > brux-nord-04 A w 0 0 0 0 > brux-sud-04 A w 0 0 0 0 > br-nord-07 Ju w 0 0 0 0 > br-sud-07 Ju w 0 0 0 0 > b-nord-10 O w 0 0 0 0 > bsud-10 O w 0,107836089 0 0,107836089 > 0,035945363 > Z1-01 J u 0 0 0 0,040567951 > Z3-01 J u 0 0 0 0 > Z5-01 J d 0,023116043 0 0 0 > Z7-01 J d 0,014130281 0 0 0 > Z9-01 J d 0 0 0 0 > Z10-01 J d 0 0 0 0 > Z12-01 J d 0 0 0 0 > Z1-04 A u 0 0 0 0 > Z3-04 A u 0 0 0 0 > Z5-04 A d 0 0 0 0 > Z7-04 A d 0 0 0 0 > Z9-04 A d 0 0,013839873 0 0 > Z10-04 A d 0 0 0 0 > Z12-04 A d 0 0 0 0 > Z1-07 Ju u 0 0 0 0 > Z3-07 Ju u 0 0 0 0 > Z5-07 Ju d 0 0 0 0 > Z7-07 Ju d 0 0 0 0 > Z9-07 Ju d 0 0 0 0 > Z10-07 Ju d 0 0 0 0 > Z12-07 Ju d 0 0,022301517 0 0 > Z1-10 O u 0 0 0 0 > Z3-10 O u 0 0 0 0 > Z5-10 O d 0 0 0 0 > Z7-10 O d 0 0 0,052924054 0 > Z9-10 O d 0 0 0,035050824 0 > Z10-10 O d 0 0 0 0,040783034 > Z12-10 O d 0 0 0 0 > > ______________________________________________ > 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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd....@cbs.dk Priv: pda...@gmail.com ______________________________________________ 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.
