如何用after_stat/stat_summary在ggplot直方图中动态添加均值竖线及标签

要实现无需硬编码的动态直方图（带均值竖线和自动定位的均值标签），可以结合stat_summary和ggplot_build()动态提取直方图统计数据，以下是两种可靠实现方式：

方式一：直接在stat_summary中动态计算坐标 #

library(tidyverse)

ggplot(mtcars, aes(x = mpg)) +
  # 自动分箱生成直方图
  geom_histogram(bins = 30, fill = "lightgray", color = "black") +
  # 均值竖线：用stat_summary计算x变量均值
  stat_summary(
    aes(xintercept = after_stat(x), y = 0),
    fun = mean, 
    geom = "vline", 
    color = "red",
    linetype = "dashed"
  ) +
  # 动态均值标签：x为均值，y为直方图最大频数的95%
  stat_summary(
    aes(
      x = after_stat(x),
      y = max(ggplot_build(plot)$data[[1]]$count) * 0.95,
      label = paste("mean:", round(after_stat(x), 1))
    ),
    fun = mean, 
    geom = "text",
    hjust = -0.1, # 调整标签位置避免与竖线重叠
    color = "red",
    size = 4
  )

核心逻辑： #

• 均值竖线：通过stat_summary的fun=mean直接计算x变量均值，用after_stat(x)引用结果作为竖线的x截距。
• 动态y坐标：ggplot_build(plot)$data[[1]]$count提取直方图的频数数据，取最大值后乘以0.95，确保标签始终位于直方图顶部下方的合理位置。
• 批量适配：无论数据集或分箱数调整，代码都会自动重新计算所有统计值，完全无需硬编码。

方式二：自定义统计函数封装逻辑 #

如果需要更清晰的代码结构，可将统计计算封装为自定义函数：

library(tidyverse)

# 自定义函数：同时计算x均值和最大频数的95%
get_hist_stats <- function(x) {
  # 提取当前直方图的频数数据
  hist_counts <- ggplot_build(current_plot())$data[[1]]$count
  tibble(
    x = mean(x, na.rm = TRUE),
    y = max(hist_counts) * 0.95,
    label = paste("mean:", round(mean(x, na.rm = TRUE), 1))
  )
}

ggplot(mtcars, aes(x = mpg)) +
  geom_histogram(bins = 30) +
  stat_summary(aes(xintercept = x, y = 0), fun = mean, geom = "vline", color = "darkgreen") +
  stat_summary(fun.data = get_hist_stats, geom = "text", hjust = -0.1, color = "darkgreen")

这种方式将统计逻辑独立封装，更便于批量复用和维护。

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