p = c("sf", "sp", "spdep")
if(any(!unlist(lapply(p, requireNamespace, quietly=TRUE)))) {
m = which(!unlist(lapply(p, requireNamespace, quietly=TRUE)))
message("Can't run examples, please install ", paste(p[m], collapse = " "))
} else {
invisible(lapply(p, require, character.only=TRUE))
data(meuse, package = "sp")
meuse <- st_as_sf(meuse, coords = c("x", "y"), crs = 28992, agr = "constant")
# Autocorrelation (univariate)
morans.plot(meuse$zinc, coords = st_coordinates(meuse)[,1:2])
# Cross-correlation of: x influencing y and y influencing x
opar <- par(no.readonly=TRUE)
par(mfrow=c(1,2))
morans.plot(x=meuse$zinc, y=meuse$copper, coords = st_coordinates(meuse)[,1:2],
scale.morans = TRUE)
morans.plot(x=meuse$zinc, y=meuse$copper, coords = st_coordinates(meuse)[,1:2],
scale.morans = TRUE, type.ac="yx")
par(opar)
}
Run the code above in your browser using DataLab