stat_compare_means won't work when using ggpubr box-plots where Y values are a list and facet.by is used

Question

I have 94 rows of samples with 26 metal concentrations per sample in columns. I also have a column for year the sample was taken (a factor with 6 levels), and one for sample