如何在分面直方图中按物种(sp)显示均值文本？

解决步骤： #

1. 预先计算各物种均值并保存为独立数据框
   把分组计算出的均值存成单独的数据框，方便后续ggplot调用：

   library(dplyr)
   library(ggplot2)
   
   set.seed(111)
   rand.data <- rnorm(1000,5,1)
   sp <- rep(c("A","B","C","D"), each =250)
   rand.df <- data.frame(rand.data,sp)
   rand.df$sp <- as.factor(rand.df$sp)
   
   # 计算各物种均值并存储
   mean_df <- rand.df %>% 
     group_by(sp) %>% 
     summarise(rand.df.mean = mean(rand.data))
   

2. 修正geom_vline显示分组均值
   原来的geom_vline调用了全局数据的均值，现在指定mean_df为数据源，让每个分面匹配对应物种的均值：

   ggplot(rand.df, aes(x=rand.data)) + 
     geom_histogram(fill = "gray", color = "black") +
     theme_classic() +
     # 基于分组均值数据绘制竖线
     geom_vline(data = mean_df, aes(xintercept=rand.df.mean), 
                color="blue", linetype="dashed", size=1) +
     facet_grid(.~sp)
   

3. 添加均值文本标签
   使用geom_text或geom_label，同样以mean_df为数据源，指定文本位置（x为均值，y设为直方图峰值附近，可根据实际数据调整），并格式化显示内容：

   ggplot(rand.df, aes(x=rand.data)) + 
     geom_histogram(fill = "gray", color = "black") +
     theme_classic() +
     geom_vline(data = mean_df, aes(xintercept=rand.df.mean), 
                color="blue", linetype="dashed", size=1) +
     # 添加格式化的均值文本
     geom_text(data = mean_df, aes(x=rand.df.mean, y=50, 
                                  label=paste0("均值: ", round(rand.df.mean, 2))),
              color="red", fontface="bold") +
     facet_grid(.~sp)
   

关键说明： #

• 必须为geom_vline、geom_text指定独立数据源mean_df，否则ggplot会默认使用主数据框rand.df，导致计算全局均值或无法匹配分组。
• y轴的数值（示例中为50）需要根据直方图实际高度调整，确保文本显示在合适位置。若想自动匹配每个分面的峰值，可先计算各分面的最大频数，再将y设为该值的一定比例（比如80%）。

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