st_coordinates()
sf
inla()
Non-spatial criteria – continued:
nc[!(nc$FIPS %in% c("37125", "37051")), ] |> ggplot() + geom_sf(aes(fill = SID79)) + theme_bw()