[R] MANOVA usage
Prof Brian Ripley
ripley at stats.ox.ac.uk
Thu Feb 22 11:47:44 CET 2007
On Thu, 22 Feb 2007, Aalim Weljie wrote:
> Hello,
>
> I had a couple questions about manova modeling in R.
>
> I have calculated a manova model, and generated a summary.manova output
> using both the Wilks test and Pillai test.
>
> The output is essentially the same, except that the Wilks lambda = 1 -
> Pillai. Is this normal? (The output from both is appended below.)
For a 1-df test, yes: all the statistics give the same test (as
?summary.manova does say).
> My other question is about the use of MANOVA. If I have one variable which
> has a higher F-stat value (~60) than the MANOVA fstat (~20), but the other
> 15 variables all have lower ANOVA F-stat values (~5-10), am I still okay in
> using the MANOVA? (If I remove the one highly significant variable, the
> MANOVA f-stat is still larger than the other univariate F-stats).
OK for what purpose? Using (M)ANOVA for variable selection has many
issues.
> Thanks very much for your help,
>
> Aalim
>
> Output from Manova tests:
>
>> summary.manova(fit)
> Df Pillai approx F num Df den Df Pr(>F)
> response 1 0.9725 19.8967 16 9 4.193e-05 ***
> Residuals 24
> ---
> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>> summary.manova(fit, test="Wilks")
> Df Wilks approx F num Df den Df Pr(>F)
> response 1 0.0275 19.8967 16 9 4.193e-05 ***
> Residuals 24
> ---
> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
>
--
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 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595
More information about the R-help
mailing list