使用ggpubr箱形图其中Y值是列表且使用facet.by时，stat_compare_means无法使用

如何解决使用ggpubr箱形图其中Y值是列表且使用facet.by时，stat_compare_means无法使用

我有94行样品，每列样品中有26种金属浓度。我还有一列用于说明年份的采样（一个因子有6个级别），还有一个用于采样位置（总计3个因子）。 我想按年份将每个元素按位置（共26个图）进行分组，每组包含一组箱形图，其中每个图面均使用stats_compare_means显示均值比较。 我一次可以做一种金属

#example data with 3 metals
my.df <- data.frame(Location = sample(c("1","2","3"),100,replace = TRUE),Year = as.factor(sample(c("2011","2012","2015","2016","2017","2018"),replace = TRUE)),Co=sample(1:100),Cd=sample(1:100),Ni=sample(1:100))

#define comparisons between all three sites in a list
my_comparisons <- list( c("1","2"),c("2",c("1","3") )

#plot cobalt data grouped by locations faceted by sample year + mean comparisons based on list
Co<-ggboxplot(my.df,x = "Location",y = "Co",color = "Location",palette = "jco",add="jitter",facet.by = "my.df$Year")+stat_compare_means(comparisons = my_comparisons)

与其对每个元素执行此操作，不如我可以创建一个包含所有元素的列标题列表，并使用该列表定义Y。我确实获得了所有图表，但是stat_compare_means的结果为NULL而不是地块列表。

#list of columns to feed Y variable
metals<-colnames(my.df[,3:5])
#feed this to the box plot
metals.all<-ggboxplot(my.df,y = metals,facet.by = "my.df$Year")+stat_compare_means(comparisons = my_comparisons)

这行不通吗？ 如果没有stat_compare_means，我会得到每种金属的图

metals.all<-ggboxplot(my.df,facet.by = "my.df$Year")