data(tintoodiel)
estuary.pca <- dudi.pca(tintoodiel$tab, scan = FALSE, nf = 4)
if(requireNamespace("spdep", quietly = TRUE) & requireNamespace("adespatial", quietly = TRUE)) {
estuary.listw <- spdep::nb2listw(tintoodiel$nb)
estuary.pca.ms <- adespatial::multispati(estuary.pca, estuary.listw, scan = FALSE,
nfposi = 3, nfnega = 2)
summary(estuary.pca.ms)
par(mfrow = c(1, 2))
barplot(estuary.pca$eig)
barplot(estuary.pca.ms$eig)
par(mfrow = c(1, 1))
}
Run the code above in your browser using DataLab