tidycensus_example.r

library("tidyverse")
library("tidycensus")

census_api_key("INSERT_API_KEY")

wright <- get_acs(
  state = "MN",
  county = "Wright",
  geography = "county subdivision",
  variables = "B19013_001",
  geometry = TRUE,
  year = 2020
) %>%
# filter out the townships
filter(str_detect(NAME,  regex("city", ignore_case = TRUE)))

# get the county outline
wright_county <- get_acs(
  state = "MN",
  county = "Wright",
  geography = "county",
  variables = "B19013_001",
  geometry = TRUE,
  year = 2020
)

ggplot() + 
geom_sf(wright_county, mapping = aes(), color = "black") +
geom_sf(wright, mapping = aes(), color = NA, fill = "white")