How to automatically adjust ylim for boxplots in for loop

Question

I have a for loop to create box plots for a couple dozen analytes grouped by location:

for(parameter in analyte_list) {
  print(parameter)
  df_analytes <-