Based on my data, I can prepare a simple dodge bar plot
geo <- data.frame("year" = c(2018, 2018, 2018, 2019, 2019, 2019, 2020, 2020, 2020),
