[R] Contrasts in manova
Hi everybody I am trying to find contrast in MANOVA. I used next code contrasts(ffage)-ctr contrasts(ffage) MANOVA.agec-manova(Y1~ffage,data=vol18.df) summary(MANOVA.agec, split =list (ffage=list(0-17 v over 18=0, 18-25 v over 26=1, 26-31 v over 32=2, 32-42 v over 43=3, 43-65 v 66+=4))) But the output was only the overall fit contrasts(ffage)-ctr contrasts(ffage) [,1] [,2] [,3] [,4] [,5] 050000 1 -14000 2 -1 -1300 3 -1 -1 -120 4 -1 -1 -1 -11 5 -1 -1 -1 -1 -1 MANOVA.agec-manova(Y1~ffage,data=vol18.df) summary(MANOVA.agec, split =list (ffage=list(0-17 v over 18=0, 18-25 v over 26=1, 26-31 v over 32=2, 32-42 v over 43=3, 43-65 v 66+=4))) Df Pillai approx F num Df den Df Pr(F) ffage 5 0.30607 2.1543 20520 0.002681 ** Residuals 130 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 I try to use linearHypothesis, but i couldnt find anything, maybe I didnt use in a corect way, *any advice is welcome* I always get the error Error in solve.default(wcrossprod(model.matrix(model), w = wts)) : Lapack routine dgesv: system is exactly singular Thanks in advance Paola -- View this message in context: http://r.789695.n4.nabble.com/Contrasts-in-manova-tp4648306.html Sent from the R help mailing list archive at Nabble.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.
[R] contrasts in MANOVA
Dear r-helpers, I have a query regarding use of contrasts in MANOVA. summary(manova(model)) gives me only result of test for overall difference. Would you be so kind and give me a hint how to get the same test statistics (e.g.Pillai's) and P values for the predefined contrasts? Best regards Ondrej Mikula -- Institute of Animal Physiology and Genetics Academy of Sciences of the Czech Republic Veveri 97, 60200 Brno, Czech Republic Institute of Vertebrate Biology Academy of Sciences of the Czech Republic Studenec 122, 67502 Konesin, Czech Republic [[alternative HTML version deleted]] __ 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.
Re: [R] contrasts in MANOVA
Dear Ondrej, You might use the linearHypothesis() function in the car package. Best, John John Fox Senator William McMaster Professor of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Ondrej Mikula Sent: October-18-11 7:34 AM To: r-help@r-project.org Subject: [R] contrasts in MANOVA Dear r-helpers, I have a query regarding use of contrasts in MANOVA. summary(manova(model)) gives me only result of test for overall difference. Would you be so kind and give me a hint how to get the same test statistics (e.g.Pillai's) and P values for the predefined contrasts? Best regards Ondrej Mikula -- Institute of Animal Physiology and Genetics Academy of Sciences of the Czech Republic Veveri 97, 60200 Brno, Czech Republic Institute of Vertebrate Biology Academy of Sciences of the Czech Republic Studenec 122, 67502 Konesin, Czech Republic [[alternative HTML version deleted]] __ 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. __ 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.
[R] Contrasts for MANOVA
Can anyone tell me how to run contrasts for MANOVA? Thanks! [[alternative HTML version deleted]] __ 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.
Re: [R] Contrasts for MANOVA
There are several ways to perform a MANOVA analysis in R. You'll have to be more specific about your approach. For example, which function are you using? -tgs On Fri, Oct 1, 2010 at 3:31 PM, Ali S alit...@yahoo.com wrote: Can anyone tell me how to run contrasts for MANOVA? Thanks! [[alternative HTML version deleted]] __ 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. [[alternative HTML version deleted]] __ 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.