Plot Staten Island neighborhood list compare Queen and Rook

library(stringr)
library(sf)
library(spdep)

# filter for si
si <- nyc_census_data_prepped[str_detect(nyc_census_data_prepped$NAM, "Richmond"),]

# store geometry of polygons
si_geom <- st_geometry(si)

# create neighbor list
si_nb_queen <- poly2nb(si)
si_nb_rook <- poly2nb(si, queen=FALSE)

# save plot as jpeg
jpeg("images/si_neighborhood_rook_queen.jpg", width = 1000, height = 800)

# set plot structure
par(mfrow = c(1,2))

# plot queen
plot(si_geom,
     main = "Queen",
     reset = FALSE,
     cex.main = 3)
plot(si_nb, si_centroids,
     add = TRUE,
     col = 2,
     lwd = 1.5)

# plot rook
plot(si_geom,
     main = "Rook",
     reset = FALSE,
     cex.main = 3)
plot(si_nb_rook, si_centroids,
     add = TRUE,
     col = 2,
     lwd = 1.5)

# close jpeg
dev.off()