ggplot分组柱状图按物种排序标注及多元素展示优化问询

1. 将ggplot柱状图按物种层级（Sp1、Sp2、Sp3）排序
2. 为分组柱状图标注所属的Specie名称
3. 优化Genotype和Treatment的分组柱状图展示，需区分不同物种的柱状图，且不使用facet或wrap，构建同时展示Genotype、Treatment、Specie三个元素的分组柱状图

ggplot(df, aes(x=Genotype, y=var, fill=Treatment)) +
  geom_bar(stat="identity", position=position_dodge(), color="black")+
  theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1, colour = "black", face = "bold"))+
  geom_text(aes(label=Specie), position = position_dodge(width = 0.9), col='blue', angle=90)

1. 数据排序处理 #

先将Genotype按照Specie的顺序重新排序，让柱状图按物种分组展示：

# 获取每个基因型对应的物种并排序
geno_specie <- unique(df[, c("Genotype", "Specie")])
geno_specie <- geno_specie[order(geno_specie$Specie, geno_specie$Genotype), ]
# 将Genotype转为有序因子，指定排序规则
df$Genotype <- factor(df$Genotype, levels = geno_specie$Genotype)

2. 优化图形展示（区分物种+标注+分组） #

通过物种设置柱状图边框颜色、Treatment设置填充色来同时区分两个维度；优化标注位置避免重叠，添加分隔线增强物种分组的视觉区分：

library(ggplot2)
library(dplyr)

# 数据预处理：计算每个基因型-处理组的均值（避免重复数据导致柱状图叠加）
df_summary <- df %>%
  group_by(Genotype, Treatment, Specie) %>%
  summarise(var_mean = mean(var), .groups = "drop")

# 按物种顺序确定基因型的展示顺序
geno_order <- df_summary %>%
  distinct(Genotype, Specie) %>%
  arrange(Specie, Genotype) %>%
  pull(Genotype)
df_summary$Genotype <- factor(df_summary$Genotype, levels = geno_order)

# 绘制优化后的柱状图
ggplot(df_summary, aes(x=Genotype, y=var_mean, fill=Treatment)) +
  # 柱状图：用物种设置边框颜色，处理组设置填充色
  geom_bar(stat="identity", position=position_dodge(width=0.9), 
           aes(color=Specie), size=1) +
  # 物种标注：放在柱状图顶部，居中展示
  geom_text(aes(label=Specie), position=position_dodge(width=0.9), 
            vjust=-0.5, size=4, fontface="bold") +
  # 添加物种分组的垂直分隔线
  geom_vline(xintercept = c(2.5, 4.5), linetype="dashed", color="gray50", size=0.8) +
  # 自定义颜色：给每个物种设置独特边框色，处理组设置填充色
  scale_color_manual(values = c("Sp1"="darkred", "Sp2"="darkblue", "Sp3"="darkgreen")) +
  scale_fill_manual(values = c("CT"="#999999", "HS"="#E69F00")) +
  # 主题优化
  theme_bw() +
  theme(
    axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1, face="bold"),
    axis.title = element_text(face="bold", size=12),
    legend.title = element_text(face="bold"),
    panel.grid.major.x = element_blank()
  ) +
  labs(x="基因型", y="变量均值", fill="处理", color="物种")

关键优化点说明 #

• 排序逻辑：将Genotype转为有序因子，强制按Sp1→Sp2→Sp3的顺序排列，实现物种层级排序
• 维度区分：用fill=Treatment区分处理组，color=Specie区分物种，无需facet即可同时展示三个元素
• 标注优化：将物种标注放在柱状图顶部，避免与x轴标签重叠，更清晰直观
• 视觉分组：添加垂直虚线分隔不同物种的基因型组，强化视觉上的物种分组
• 数据处理：先计算均值，避免原始重复数据导致柱状图叠加，符合统计展示逻辑

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