tidycensus_density.R

library(tidycensus)
library(tidyverse)

dc_density <- get_acs(
  geography = "tract",
  state = "DC",
  variables = "B01003_001",
  geometry = TRUE,
  keep_geo_vars = TRUE
) %>%
  mutate(density_km2 = estimate / (ALAND / 1000000))

ggplot(dc_density) + 
  geom_sf(aes(fill = density_km2)) + 
  theme_void() + 
  scale_fill_viridis_c() + 
  labs(
    title = "Population density in Washington, DC",
    fill = "Pop/km2"
  )