For this project, I used Median Household Income by State Dollars in 2010 and in 2015 data from GeoFRED. The data files are stored in the same directory.

Maryland state has highest median house income in both 2010 and 2015. According to the color scale, in 2010, the states that have median household income greater than 60,000 dollars are MD, NJ, CT, MA, VA, AK, and HI. In 2015, the states that have median household income greater than 70,000 dollars are MD, NJ, CT, DC, AK, and HI. The overall median household income increase by 10,000 dollars in 2015 comparing to 2010.

library(pacman)
p_load(tidyverse, statebins, readxl)
Median_Household_Income_by_State_Dollars_2010 <- read_excel("Median_Household_Income_by_State_Dollars_2010.xls", 
    skip = 1)

med_income_2010 <- Median_Household_Income_by_State_Dollars_2010 %>% select(region = `Region Name`, value = '2010')

statebins_continuous(state_data = med_income_2010, state_col = "region",
                     text_color = "black", value_col = "value",
                     font_size = 3,
                     legend_title="Dollars",
                     legend_position = "bottom",
                     brewer_pal = "PuBu", 
                     plot_title = "Median Household Income 2010", 
                     title_position = "top")

med_income_2010 %>% arrange(desc(value)) %>% head()
Median_Household_Income_by_State_Dollars_2015 <- read_excel("Median_Household_Income_by_State_Dollars_2015.xls", 
    skip = 1)

med_income_2015 <- Median_Household_Income_by_State_Dollars_2015 %>% select(region = `Region Name`, value = '2015')

statebins_continuous(state_data = med_income_2015, state_col = "region",
                     text_color = "black", value_col = "value",
                     font_size = 3,
                     legend_title="Dollars",
                     legend_position = "bottom",
                     brewer_pal = "PuBu", 
                     plot_title = "Median Household Income 2015", 
                     title_position = "top")

med_income_2015 %>% arrange(desc(value)) %>% head()