在ggplot2的geom_bar(stat="identity")中添加总计标签的问题

问题背景 #

使用R 4.3.1（R Studio），在ggplot2的geom_bar(stat="identity")柱状图中添加柱子总计数值标签时遇到以下问题：

• 使用stat="count"搭配label=after_stat(count)时触发报错
• 直接以APIDAE作为标签，每个柱子内显示分组内多个值而非总和

基础代码：

library(ggplot2)

ggplot(gold, aes(x=Site, y=APIDAE, fill=Site, colours))+
  geom_bar(stat = "identity")+
  scale_fill_manual(name = "Sites", values = gold_cols)+
  ylab(label = "Famille apidae")+
  ylim(0, 50)+
  ggtitle("Présence Apidae entre sites")

遇到的错误信息：

Error in geom_text():
! Problem while computing stat.
ℹ Error occurred in the 2nd layer.
Caused by error in setup_params():
! stat_bin() must only have an x or y aesthetic.
Run rlang::last_trace() to see where the error occurred.

Error in geom_bar():
! Problem while computing stat.
ℹ Error occurred in the 1st layer.
Caused by error in setup_params():
! stat_count() must only have an x or y aesthetic.
Run rlang::last_trace() to see where the error occurred.

错误原因 #

1. stat="count"使用场景错误：stat_count()/stat_bin()用于对数据计数统计，要求仅指定x或y其中一个美学映射。但你的代码已指定y=APIDAE并使用stat="identity"（直接调用y轴数值），两者逻辑冲突导致报错。
2. 直接用APIDAE作标签的问题：原始数据中每个Site对应多条APIDAE记录，直接映射会显示所有单条记录值，而非分组总和。

解决方案 #

方案1：提前汇总数据（推荐） #

先用dplyr按Site分组计算APIDAE总和，再基于汇总后的数据绘图：

library(ggplot2)
library(dplyr)

# 汇总数据
gold_summary <- gold %>%
  group_by(Site) %>%
  summarize(total_APIDAE = sum(APIDAE, na.rm = TRUE))

# 绘图并添加标签
ggplot(gold_summary, aes(x=Site, y=total_APIDAE, fill=Site))+
  geom_bar(stat = "identity")+
  geom_text(aes(label = total_APIDAE), vjust = -0.5, size = 4)+ # vjust控制标签垂直位置
  scale_fill_manual(name = "Sites", values = gold_cols)+
  ylab(label = "Famille apidae")+
  ylim(0, 55)+ # 扩大y轴范围避免标签被截断
  ggtitle("Présence Apidae entre sites")

方案2：绘图时实时计算总和 #

无需提前汇总，用geom_text配合stat="summary"直接计算分组总和：

library(ggplot2)

ggplot(gold, aes(x=Site, y=APIDAE, fill=Site))+
  geom_bar(stat = "identity")+
  # 用stat="summary"计算总和并添加标签
  geom_text(stat = "summary", fun = sum, aes(label = after_stat(y)), vjust = -0.5, size = 4)+
  scale_fill_manual(name = "Sites", values = gold_cols)+
  ylab(label = "Famille apidae")+
  ylim(0, 55)+
  ggtitle("Présence Apidae entre sites")

额外注意点 #

• 原代码aes()中的colours是无效参数，已移除，避免不必要的警告或错误。
• 调整ylim上限值确保标签不超出绘图区域；vjust参数可调整标签相对柱子的位置（负数向上，正数向下）。

内容的提问来源于stack exchange，提问作者Shawna O'Flaherty