if(!require(tidyverse))install.packages(“tidyverse”)

library(tidyverse)

Inc2022 <- read.csv(“https://drkblake.com/wp-content/uploads/2024/02/Inc2022.csv”)

Inc2017 <- read.csv(“https://drkblake.com/wp-content/uploads/2024/02/Inc2017.csv”)

Income <- left_join(Inc2022,Inc2017, by = join_by(GEOID == GEOID))

Income <- Income %>% mutate(Change = HHInc2022 - HHInc2017) head(Income, 10)

Income <- Income %>% mutate(Level = case_when(HHInc2022 > 99999 ~ “$100k+”, HHInc2022 < 100000 ~ “<$100k”, .default = “Error”))

LevelbyCounty <- Income %>% group_by(County, Level) %>% summarize(Count = n()) %>% pivot_wider(names_from = Level, values_from = Count)

RichDistricts <- Income %>% filter(Level == “$100k+”) head(RichDistricts, 10)

RichDistricts <- RichDistricts %>% arrange(desc(HHInc2022)) head(RichDistricts, 10)