ade4 (version 1.7-22)

costatis.randtest: Monte-Carlo test on a Costatis analysis (in C).

Description

Performs a Monte-Carlo test on a Costatis analysis.

Usage

costatis.randtest(KTX, KTY, nrepet = 999, ...)

Value

a list of the class randtest

Arguments

KTX

an objet of class ktab

KTY

an objet of class ktab

nrepet

the number of permutations

...

further arguments passed to or from other methods

Author

Jean Thioulouse

References

Thioulouse J. (2011). Simultaneous analysis of a sequence of paired ecological tables: a comparison of several methods. Annals of Applied Statistics, 5, 2300-2325.

Examples

data(meau)
wit1 <- withinpca(meau$env, meau$design$season, scan = FALSE, scal = "total")
pcaspe <- dudi.pca(meau$spe, scale = FALSE, scan = FALSE, nf = 2)
wit2 <- wca(pcaspe, meau$design$season, scan = FALSE, nf = 2)
kta1 <- ktab.within(wit1, colnames = rep(c("S1","S2","S3","S4","S5","S6"), 4))
kta2 <- ktab.within(wit2, colnames = rep(c("S1","S2","S3","S4","S5","S6"), 4))
costatis1 <- costatis(kta1, kta2, scan = FALSE)
costatis.randtest(kta1, kta2)