R Notebook

# Replicate BIA for Q1 income
for (i in unique(df$income_var)) {
  p <- ggplot(filter(df, income_var == i), aes(quints, fill = program, y = share_in_quant)) +
  geom_col(position = "dodge") +
    labs(title = i)+
    scale_fill_manual(values = c("sap" = "lightblue", "gov_relief" = "darkblue", "any_4ps" = "darkgreen"))
  print(p)
}