This is the coding that I have right now that creates the graph that I want. This is only for well 'aa1' though, and I want to learn how to make a loop so that I can make this graph for all of my wells. Any ideas?
longer_raindata %>%
select(well, metal, level, smpl_mth) %>%
filter(well == 'aa1') %>%
ggplot(aes(metal, level, fill = smpl_mth))+
scale_fill_manual(values = c('plum', 'lightsteelblue', 'darkolivegreen', 'indianred'))+
geom_col(position = "dodge")+
labs(title = 'Metals in Well AA1',
x = 'Metals',
y = 'µg/l')