# NOT RUN {
library(ggplot2)
library(sf)
# Convert all 2008 plot stems to sf object
census_2008_bw_sf <- census_2008_bw %>%
st_as_sf(coords = c("gx", "gy"))
# Plot stems with plot boundary
ggplot() +
geom_sf(data = census_2008_bw_sf, size = 0.25)
# }
Run the code above in your browser using DataLab