Mapping proportion of those with West Indian Ancestry in BK, NY

This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.

Try executing this chunk by clicking the Run button within the chunk or by placing your cursor inside it and pressing Cmd+Shift+Enter.

library(tidyverse)
library(tidycensus)
library(sf)
library(scales)
library(RColorBrewer)
library(plotly)
library(knitr)

Methods

Borough boundaries are represented with a shapefile that was downloaded from NYC Open Data. Census tract boundaries and data are from the 2016-2020 5-year American Community Survey, accessed with the tidyverse R package. Census data includes: -the number of people with West Indian ancestry -total population that reports ancestry The proportion of people with West Indian ancestry for each census tract was calculated and the census tracts with no residents were removed.

boros <- st_read("data/raw/geo/Borough Boundaries.geojson")

raw_ancestry <- suppressWarnings(get_acs(geography = "tract",
  variables = c(ancestry_pop = "B04006_001",
  west_indian = "B04006_094"),
  state='NY',
  county = 'Kings',
  geometry = T,
  year = 2020,
  output = "wide"))
west_indian <- raw_ancestry |>
mutate(pct_west_indian = west_indianE/ancestry_popE) |>
filter (ancestry_popE > 0)

Results

west_indian_map <- ggplot(data = west_indian, mapping = aes(fill = pct_west_indian)) +
  suppressWarnings(geom_sf(data = west_indian,
          mapping = aes(fill = pct_west_indian,
                        text = paste0(NAME, ":",
                                      "<br>Perncet West Indian Ancestry",
                                      scales::percent(pct_west_indian,
                                                      accuracy=1))),
    color = "transparent")) +
  theme_void() +
  scale_fill_distiller (breaks=c(0, .2, .4, .6, .8, 1),
  direction = 1,
  na.value = "#fafafa",
  name="Percent West Indian Ancestry (%)",
  labels=percent_format (accuracy = 1L)) +
  labs (title = "Brooklyn, West Indian Ancestry by Census Tract",
caption = "Source: American Community Survey, 2016-20") +
suppressWarnings(geom_sf (data = boros |> filter(boro_name == "Brooklyn"),
color = "black", fill = NA, Iwd = .5))
ggplotly(west_indian_map, tooltip = "text")