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')
 
     
    

