Skip to content

Instantly share code, notes, and snippets.

@walkerke
Created August 25, 2018 14:18
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 1 You must be signed in to fork a gist
  • Save walkerke/302e938ee6096fd181e199668b1b57bc to your computer and use it in GitHub Desktop.
Save walkerke/302e938ee6096fd181e199668b1b57bc to your computer and use it in GitHub Desktop.
library(tigris)
library(tidycensus)
library(tidyverse)
library(sf)
ma_income <- get_acs(geography = "county subdivision",
variables = "B19013_001",
state = "MA")
ma_subs <- county_subdivisions(state = "MA", cb = TRUE, class = "sf") %>%
select(GEOID)
ma_joined <- left_join(ma_subs, ma_income, by = "GEOID")
plot(ma_joined["estimate"])
@walkerke
Copy link
Author

image

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment