ggplot分组箱线图修改汇总统计量（自定义分位数）的实现方法

我想要修改如下箱线图中展示的汇总统计量：

我已通过如下代码生成箱线图：

ggplot(as.data.frame(beta2), aes(y=var1,x=as.factor(Year))) + 
  geom_boxplot(outlier.shape = NA)+
  ylab(expression(beta[1]))+
  xlab("\nYear")+
  theme_bw()

默认箱线图的箱体对应第一和第三分位数，我想要让箱体展示2.5%和97.5%分位数。我了解单组箱线图可以直接在geom_boxplot中添加如下参数修改展示内容：

aes(
    ymin= min(var1),
    lower = quantile(var1,0.025),
    middle = mean(var1),
    upper = quantile(var1,0.975),
    ymax=max(var1))

但该方法在生成分组箱线图时不生效，请问如何实现分组场景下的自定义汇总统计量展示？可使用鸢尾花（Iris）数据集示例演示：

ggplot(iris, aes(y=Sepal.Length,x=Species)) + 
  geom_boxplot(outlier.shape = NA)

补充说明 #

已采纳的答案确实可以实现需求，但我的数据框体量极大，该方法运行耗时较长，我找到了另一个适配大数据集、满足我特定需求的解决方案，可在RStudio社区搜索关键词「quantile box plot which is not an outlier box plot」查询对应内容。

内容的提问来源于stack exchange，提问作者username97