# NOT RUN {
data(covid19sf_vaccine_geo)
head(covid19sf_vaccine_geo)
library(sf)
library(dplyr)
df <- covid19sf_vaccine_geo %>% filter(area_type == "Analysis Neighborhood")
plot(df[, c("percent_pop_series_completed", "geometry")],
main = "San Francisco - Percentage of Fully Vaccinated Population")
# }
Run the code above in your browser using DataLab