library(tidycensus)
library(tidyverse)
vars_2020 <- c(total_housing_units = "H1_001N",
total_vacant = "H1_003N")
fl_housing_data <- get_decennial(geography = "county",
state = "FL",
variables = vars_2020,
geometry = TRUE,
year = 2020)
fl_housing_data <- fl_housing_data |>
pivot_wider(names_from = variable, values_from = value) |>
mutate(percent_vacant = (total_vacant / total_housing_units) * 100)
fl_housing_data |>
ggplot(aes(fill = percent_vacant)) +
geom_sf(color = NA) +
scale_fill_viridis_c(option = "magma")