geom_histogram中ID标签与填充色不匹配问题求助

使用给定数据集和代码生成直方图时，出现**ID标签的颜色未对应其所属代谢组（gengroup）**的问题：ID的X轴位置正确，但标签颜色和对应条形的填充色不匹配（如ID72属于PM组却未显示深绿色）。

原数据集代码 #

set.seed(4)
df <- data.frame(ID = factor(1:72), gengroup = c("UM","NM" ,"IM", "PM"), value = 2 - rgamma(72, 3, 2))

注：原代码中gengroup仅4个值，R会自动重复填充至72行，建议显式用rep(c("UM","NM","IM","PM"), each=18)避免警告。

原直方图代码 #

p1 <- ggplot(df, aes(x = value, fill = gengroup)) +
  scale_fill_brewer(aes(name= "Metabolizer group"),  palette = "Paired", labels= c("UM","NM" ,"IM", "PM"))+
  geom_histogram(bins = 30) +
  stat_bin(geom = "text", bins = 30,size =2, na.rm = TRUE,
           aes(label = ifelse(after_stat(count) == 0, NA, after_stat(group)), 
               group = ID, y = after_stat(count)),
           position = position_stack(vjust = 0.5))   +
  labs(x = NULL)

show(p1)

问题根源 #

1. stat_bin是对数据按bin聚合统计，after_stat(group)指向的是bin的分组，而非每个ID所属的gengroup，导致标签无法关联正确的组颜色。
2. 原代码中scale_fill_brewer的name参数错误放在aes()内，导致图例名称不生效。

修正后的代码 #

library(ggplot2)
library(dplyr)

set.seed(4)
# 修正gengroup的生成方式，避免自动重复的警告
df <- data.frame(
  ID = factor(1:72),
  gengroup = rep(c("UM","NM","IM","PM"), each = 18),
  value = 2 - rgamma(72, 3, 2)
)

# 计算每个value对应的bin中点，用于标签的X轴定位
bin_info <- ggplot_build(ggplot(df, aes(x = value)) + geom_histogram(bins = 30))$data[[1]]
df$bin_mid <- findInterval(df$value, bin_info$xmin) %>%
  sapply(function(x) mean(c(bin_info$xmin[x], bin_info$xmax[x])))

# 计算每个bin内的ID堆叠位置，避免标签重叠
df <- df %>%
  group_by(bin_mid) %>%
  mutate(y_pos = row_number()) %>%
  ungroup()

p1 <- ggplot(df, aes(x = value, fill = gengroup)) +
  # 修正scale_fill_brewer的参数写法
  scale_fill_brewer(name = "Metabolizer group", palette = "Paired", labels = c("UM","NM","IM","PM")) +
  geom_histogram(bins = 30) +
  # 用geom_text直接映射每个ID的gengroup，确保颜色匹配
  geom_text(
    aes(x = bin_mid, y = y_pos, label = ID, color = gengroup),
    size = 2, na.rm = TRUE,
    position = position_nudge(y = 0.2) # 微调Y位置，让标签在条形上方/内部居中
  ) +
  # 让文字颜色与填充色保持一致
  scale_color_brewer(name = "Metabolizer group", palette = "Paired", labels = c("UM","NM","IM","PM")) +
  labs(x = NULL)

show(p1)

关键调整说明 #

• 提前计算每个ID对应的bin中点，确保标签X轴位置完全匹配条形位置
• 替换stat_bin为geom_text，直接基于原始数据的gengroup映射颜色，保证标签颜色与条形填充色一一对应
• 修正scale_fill_brewer的参数写法，让图例名称正常显示
• 计算每个bin内的堆叠位置y_pos，避免标签重叠，同时用position_nudge微调位置提升可读性

内容的提问来源于stack exchange，提问作者Sweilem B.