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)