R语言ggplot分组柱状图：为组与单柱添加X轴标签

问题描述 #

需要制作3个分组（condition：normal、stress、Nitrogen）各含4个柱子（specie：sorgho、poacee、banana、triticum）的偏移柱状图，给每个柱子添加独立标签，同时保留分组的主标签。尝试使用scale_x_discrete后所有X轴标签消失，原代码如下：

# create a dataset 
specie <- c(rep("sorgho" , 3) , rep("poacee" , 3) , rep("banana" , 3) , rep("triticum" , 3) )
condition <- rep(c("normal" , "stress" , "Nitrogen") , 4)
value <- abs(rnorm(12 , 0 , 15))
df <- data.frame(specie,condition,value)
df


# Grouped
set.seed(123)
df

ggplot(df, aes(fill = specie, y = value, x = condition)) +
  geom_bar(stat="identity", position = position_dodge(width = 0.9), width = 0.6) +
  scale_fill_manual(
    "specie",
    values = c("black", "black", "black", "black"),
    
  ) +
  theme_classic() +
  theme(panel.border = element_rect(colour = "black", fill=NA, size=1),
        text = element_text(size = 8, family = "sans"),
  ) + labs(x = "Nitrogen                  normal          stress", y = expression(paste("Total Phytoplankton Abundance (", x10^{6}, "cells/L)"))) +
  scale_y_continuous(expand = c(0,0))  +   
  theme(legend.position="none") + scale_x_discrete("", breaks=c(-0.25,0.25,0.75,1.25,1.75,2.25,2.75,3.25,3.75,4.25,4.75,5.25),
                     labels=rep(c("sorgho","poacee", "banana", "triticum"),times=3)) 

问题原因 #

原代码中scale_x_discrete的breaks参数使用了数值，但X轴映射的是离散分类变量condition，数值型breaks无法与离散变量的水平匹配，导致标签全部无法显示。

解决方案 #

将X轴转换为连续型变量，精准控制每个偏移柱子的位置，同时添加分组主标签。修改后的代码如下：

set.seed(123)
# 构建数据集
specie <- c(rep("sorgho" , 3) , rep("poacee" , 3) , rep("banana" , 3) , rep("triticum" , 3) )
condition <- rep(c("normal" , "stress" , "Nitrogen") , 4)
value <- abs(rnorm(12 , 0 , 15))
df <- data.frame(specie,condition,value)

# 给每个condition分配连续型数值位置
cond_levels <- unique(df$condition)
df$x_pos <- as.numeric(factor(df$condition, levels = cond_levels))

# 定义偏移宽度（与geom_bar的position_dodge参数一致）
dodge_width <- 0.9

ggplot(df, aes(fill = specie, y = value, x = x_pos)) +
  geom_bar(stat="identity", position = position_dodge(width = dodge_width), width = 0.6) +
  scale_fill_manual(values = rep("black", 4)) +
  theme_classic() +
  theme(
    panel.border = element_rect(colour = "black", fill=NA, size=1),
    text = element_text(size = 8, family = "sans"),
    legend.position = "none"
  ) +
  labs(y = expression(paste("Total Phytoplankton Abundance (", x10^{6}, "cells/L)"))) +
  scale_y_continuous(expand = c(0,0)) +
  # 设置每个柱子的X轴标签
  scale_x_continuous(
    name = "",
    # 计算每个偏移柱子的位置：分组中心位置 + 偏移量
    breaks = rep(1:length(cond_levels), each=4) + rep(seq(-dodge_width/2, dodge_width/2, length.out=4), times=3),
    labels = rep(c("sorgho","poacee", "banana", "triticum"), times=3)
  ) +
  # 添加分组主标签，放在每个分组的中心下方
  annotate("text", x = 1:length(cond_levels), y = -2, label = cond_levels, size = 3) +
  # 调整Y轴范围，确保分组标签不被裁剪
  coord_cartesian(ylim = c(0, max(df$value) + 2), clip = "off")

关键修改点 #

• 将离散的condition转换为连续型数值x_pos，便于精准控制柱子位置
• 用scale_x_continuous设置每个偏移柱子的breaks和对应标签，匹配柱子的实际位置
• 通过annotate添加分组主标签，结合coord_cartesian调整显示范围避免标签被裁剪
• 移除原手动编写的X轴标题，改用自动生成的分组标签

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