Flussi cartografici con tigris e tidycensus

Analizzare i dati del Censimento USA in R

Kyle Walker

Instructor

Mappa a punti sull'istruzione

Generare punti casuali con sf

Funzione chiave per generare punti casuali in sf: st_sample()

dc_dots <- map(c("White", "Black", "Hispanic", "Asian"), function(group) {
  dc_race %>%
    filter(variable == group) %>%
    st_sample(., size = .$value / 100) %>%
    st_sf() %>%
    mutate(group = group) 
}) %>%
  reduce(rbind)

Considerazioni sulla generazione casuale di punti

Per disegnare più velocemente:

dc_dots <- dc_dots %>%
  group_by(group) %>%
  summarize()

Per visualizzazioni più accurate:

dc_dots_shuffle <- sample_frac(dc_dots, size = 1)

Mappa a densità di punti di base con sf

plot(dc_dots_shuffle, key.pos = 1)

Mappa a densità di punti di base

Dati ausiliari con tigris

options(tigris_class = "sf")

dc_roads <- roads("DC", "District of Columbia") %>%
  filter(RTTYP %in% c("I", "S", "U"))

dc_water <- area_water("DC", "District of Columbia")

dc_boundary <- counties("DC", cb = TRUE)

Dati ausiliari con tigris

plot(dc_water$geometry, col = "lightblue")

Grafico delle acque di DC

Mappa a densità di punti con ggplot2

ggplot() + 
  geom_sf(data = dc_boundary, color = NA, fill = "white") + 
  geom_sf(data = dc_dots, aes(color = group, fill = group), size = 0.1) + 
  geom_sf(data = dc_water, color = "lightblue", fill = "lightblue") + 
  geom_sf(data = dc_roads, color = "grey") + 
  coord_sf(crs = 26918, datum = NA) + 
  scale_color_brewer(palette = "Set1", guide = FALSE) +
  scale_fill_brewer(palette = "Set1") +
  labs(title = "The racial geography of Washington, DC", 
       subtitle = "2010 decennial U.S. Census", 
       fill = "", 
       caption = "1 dot = approximately 100 people.\nData acquired with 
                  the R tidycensus and tigris packages.")

Mappa finale a densità di punti

Considerazioni per le mappe a densità di punti

Attenzione ai modi in cui le mappe a densità di punti possono essere fraintese
Scegli con cura colori qualitativi
Valuta bene la scelta dei layer ausiliari

Esercitiamoci!

Analizzare i dati del Censimento USA in R