Momocs (version 1.4.1)

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, retain = 0.99, ...)

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

• stars.tab a table with 'significance stars', discutable but largely used: '' if Pr(>F) < 0.001; '' of < 0.01; '' if < 0.05; '.' if < 0.10 and '-' if above.

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.

retain

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

See Also

MANOVA, manova.

Other multivariate: CLUST(), KMEANS(), KMEDOIDS(), LDA(), MANOVA(), MDS(), MSHAPES(), NMDS(), PCA(), classification_metrics()

Examples

# 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
op <- PCA(npoly(olea))
MANOVA_PW(op, 'domes')
# to get the results
res <- MANOVA_PW(op, 'domes')
res$manovas
res$stars.tab
res$summary