Multivariate spatial autocorrelation test
This function performs a multivariate autocorrelation test.
multispati.rtest(dudi, listw, nrepet = 99, ...)
dudi |
an object of class |
listw |
an object of class |
nrepet |
the number of permutations |
... |
further arguments passed to or from other methods |
We note X the data frame with the variables, Q the column weight matrix and D the row weight matrix associated to the duality diagram dudi. We note L the neighbouring weights matrix associated to listw. This function performs a Monte-Carlo Test on the multivariate spatial autocorrelation index :
r = t(X)DLXQ / t(X)DXQ
Returns an object of class randtest (randomization tests).
Daniel Chessel
Sébastien Ollier sebastien.ollier@u-psud.fr
Smouse, P. E. and Peakall, R. (1999) Spatial autocorrelation analysis of individual multiallele and multilocus genetic structure. Heredity, 82, 561–573.
if (requireNamespace("spdep", quietly = TRUE)) {
data(mafragh)
maf.listw <- spdep::nb2listw(neig2nb(mafragh$neig))
maf.pca <- dudi.pca(mafragh$env, scannf = FALSE)
multispati.rtest(maf.pca, maf.listw)
maf.pca.ms <- multispati(maf.pca, maf.listw, scannf = FALSE)
plot(maf.pca.ms)
}Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.