multispati.rtest {ade4} | R Documentation |
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 |
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("maptools", quiet = TRUE) & requireNamespace("spdep", quiet = TRUE)) { data(mafragh) maf.listw <- nb2listw(neig2nb(mafragh$neig)) maf.pca <- dudi.pca(mafragh$mil, scannf = FALSE) multispati.rtest(maf.pca, maf.listw) maf.pca.ms <- multispati(maf.pca, maf.listw, scannf=FALSE) plot(maf.pca.ms) }