median_age_mapdeck.R

library(tidycensus)
library(mapdeck)
library(sf)
options(tigris_use_cache = TRUE)

tract_median_age <- get_acs(
  geography = "tract",
  variables = "B01002_001",
  state = c(state.abb, "DC", "PR"),
  geometry = TRUE
) %>%
  dplyr::filter(!is.na(st_dimension(.)))

mapdeck(style = mapdeck_style("light"),
        token = Sys.getenv("MAPBOX_TOKEN")) %>%
  add_polygon(data = tract_median_age,
              fill_colour = "estimate",
              legend = TRUE,
              tooltip = "estimate",
              auto_highlight = TRUE, 
              legend_options = list(
                title = "Median age"
              ))