if (requireNamespace("sf", quietly = TRUE)) {
if (requireNamespace("spdep", quietly = TRUE)) {
library(spdep)
nc.sids <- sf::st_read(system.file("shapes/sids.shp", package="spData")[1])
nc.sids <- as(nc.sids, "Spatial")
sp::proj4string(nc.sids) <- sp::CRS("+proj=longlat +ellps=clrk66")
row.names(nc.sids) <- as.character(nc.sids$FIPS)
rn <- row.names(nc.sids)
ncCC89_nb <- read.gal(system.file("weights/ncCC89.gal", package="spData")[1],
region.id=rn)
ncCR85_nb <- read.gal(system.file("weights/ncCR85.gal", package="spData")[1],
region.id=rn)
sp::plot(nc.sids, border="grey")
plot(ncCR85_nb, sp::coordinates(nc.sids), add=TRUE, col="blue")
sp::plot(nc.sids, border="grey")
plot(ncCC89_nb, sp::coordinates(nc.sids), add=TRUE, col="blue")
}
}
Run the code above in your browser using DataLab