Momocs (version 1.1.6)

MANOVA_PW: Pairwise Multivariate analyses of variance

Description

A wrapper for pairwise MANOVAs on Coe objects. Calculates a MANOVA for every pairwise combination of the factor provided.

Usage

MANOVA_PW(x, ...)
# S3 method for PCA
MANOVA_PW(x, fac, verbose = FALSE, retain = 0.99, ...)

Arguments

x

a PCA object

...

more arguments to feed MANOVA

fac

a name (or its id) of a grouping factor in $fac or a factor or a formula.

verbose

to feed MANOVA

retain

the number of PC axis to retain (1:retain) or the proportion of variance to capture (0.99 par default).

Value

a list with the following components is returned (invisibly because $manovas may be very long, see examples):

• manovas a list containing all the raw manovas
• summary a matrix with all important statists
• stars.tab a table with 'significance star', discutable but useful: '***' if Pr(>F) < 0.001; '**' of < 0.01; '*' if < 0.05; '.' if < 0.10 and '-' if above.

See Also

MANOVA, manova. Other multivariate: CLUST, KMEANS, LDA, MANOVA, PCA

Examples

data(bot)
# we create a fake factor with 4 levels
bot$fac$fake <- factor(rep(letters[1:4], each=10))
bot.p <- PCA(efourier(bot, 8))
MANOVA_PW(bot.p, 'fake') # or MANOVA_PW(bot.p, 2)

# an example on open outlines
data(olea)
op <- PCA(npoly(olea))
MANOVA_PW(op, 'domes')
# to get the results
res <- MANOVA_PW(op, 'domes')
res$manovas
res$stars.tab
res$summary