procuste.randtest {ade4} | R Documentation |
performs a Monte-Carlo Test on the sum of the singular values of a procustean rotation.
procuste.randtest(df1, df2, nrepet = 999)
df1 |
a data frame |
df2 |
a data frame |
nrepet |
the number of permutations |
returns a list of class randtest
Jean Thioulouse Jean.Thioulouse@univ-lyon1.fr
Jackson, D.A. (1995) PROTEST: a PROcustean randomization TEST of community environment concordance. Ecosciences, 2, 297–303.
data(doubs) pca1 <- dudi.pca(doubs$env, scal = TRUE, scann = FALSE) pca2 <- dudi.pca(doubs$fish, scal = FALSE, scann = FALSE) protest1 <- procuste.randtest(pca1$tab, pca2$tab, 999) protest1 plot(protest1,main="PROTEST")