West Indian Ancestry

Methods Section

Here we are making a choropleth map to show the West Indian ancestry in Brooklyn

library(tidyverse)
library(tidycensus)
library(sf)
library(scales)
library(viridisLite)
##library(plotly)

ggplot(data = west_indian, mapping = aes(fill = pct_west_indian))  + 
  geom_sf(color = "#ffffff") +
  theme_void() +
  scale_fill_distiller(breaks=c(0, .2, .4, .6, .8, 1),
                       direction = 1,
                       na.value = "transparent",
                       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")+ 
  geom_sf(data = nabes |> filter(BoroName == "Brooklyn"), 
          color = "gray", fill = NA, lwd = 0.25) + 
  geom_sf(data = boros |> filter(boro_name == "Brooklyn"), 
          color = "black", fill = NA, lwd = .5)

##I keep getting an error message for this part and I’m not quite sure why

west_indian_2263 <- st_transform(west_indian, 2263)

#remove unnecessary fields in the neighborhood shapefile
nabes_selected <- nabes |>
  select(BoroCode, BoroName, NTA2020, NTAName)

west_indian_nabes <- west_indian_2263 |>
  st_join(nabes_selected, 
          left = TRUE,
          join = st_intersects,
          largest = TRUE)

Warning: attribute variables are assumed to be spatially constant throughout all geometries

ggplot(data = west_indian_nabes |> 
         filter(NTAName == "Crown Heights (North)"),
       mapping = aes(fill = pct_west_indian))  + 
  geom_sf(color = "#ffffff",
          lwd = 0) +
  theme_void() +
  scale_fill_distiller(breaks=c(0, .2, .4, .6, .8, 1),
                       direction = 1,
                       na.value = "transparent",
                       name="Percent West Indian Ancestry (%)",
                       labels=percent_format(accuracy = 1L)) +
  labs(
    title = "Brooklyn, Crown Heights, West Indian Ancestry by Census Tract",
    caption = "Source: American Community Survey, 2016-20"
  ) + 
  geom_sf(data = nabes |> filter(NTAName == "Crown Heights (North)"), 
          color = "black", fill = NA, lwd = 0.5)