procuste.rtest {ade4} R Documentation

Monte-Carlo Test on the sum of the singular values of a procustean rotation (in R).

Description

performs a Monte-Carlo Test on the sum of the singular values of a procustean rotation.

Usage

procuste.rtest(df1, df2, nrepet = 99)

Arguments

 df1 a data frame df2 a data frame nrepet the number of permutations

Value

returns a list of class rtest

Author(s)

Daniel Chessel
Anne B Dufour dufour@biomserv.univ-lyon1.fr

References

Jackson, D.A. (1995) PROTEST: a PROcustean randomization TEST of community environment concordance. Ecosciences, 2, 297–303.

Examples

data(doubs)
pca1 <- dudi.pca(doubs\$mil, scal = TRUE, scann = FALSE)
pca2 <- dudi.pca(doubs\$poi, scal = FALSE, scann = FALSE)
proc1 <- procuste(pca1\$tab, pca2\$tab)
protest1 <- procuste.rtest(pca1\$tab, pca2\$tab, 999)
protest1
plot(protest1)

Worked out examples

> ### Name: procuste.rtest
> ### Title: Monte-Carlo Test on the sum of the singular values of a
> ###   procustean rotation (in R).
> ### Aliases: procuste.rtest
> ### Keywords: multivariate nonparametric
>
> ### ** Examples
>
> data(doubs)
> pca1 <- dudi.pca(doubs\$mil, scal = TRUE, scann = FALSE)
> pca2 <- dudi.pca(doubs\$poi, scal = FALSE, scann = FALSE)
> proc1 <- procuste(pca1\$tab, pca2\$tab)
> protest1 <- procuste.rtest(pca1\$tab, pca2\$tab, 999)
> protest1
Monte-Carlo test
Observation: 0.6562
Call: procuste.rtest(df1 = pca1\$tab, df2 = pca2\$tab, nrepet = 999)
Based on 999 replicates
Simulated p-value: 0.001
> plot(protest1)
>
>
>
>