if (requireNamespace("spdep")) {
## generate adjacency matrix for districts of Bayern and Baden-Wuerttemberg
data("fluBYBW")
adjmat <- poly2adjmat(fluBYBW@map)
## same as already stored in the neighbourhood slot (in different order)
stopifnot(all.equal(adjmat,
neighbourhood(fluBYBW)[rownames(adjmat),colnames(adjmat)]))
## the neighbourhood graph can be plotted with spdep
plot(spdep::mat2listw(adjmat), coordinates(fluBYBW@map))
}
Run the code above in your browser using DataLab