if (require(maptools, quiet = TRUE) & require(spdep, quiet = TRUE)) {
data(mafragh)
maf.xy <- mafragh$xy
maf.flo <- mafragh$flo
maf.listw <- nb2listw(neig2nb(mafragh$neig))
s.label(maf.xy, neig = mafragh$neig, clab = 0.75)
maf.coa <- dudi.coa(maf.flo,scannf = FALSE)
multispati.randtest(maf.coa, maf.listw)
}
Run the code above in your browser using DataLab