双因素方差分析柱状图重叠：如何分离并实现单Y轴同图展示

我明白你遇到的问题——当尝试把nozzle和adj两个因素整合到同一柱状图时，柱子容易出现重叠，而且之前用分面的方式没法满足统一Y轴的需求。下面是两种简便的解决方案，既能完美解决重叠问题，又能把双因素展示在同一图中：

方案一：用颜色区分喷嘴，填充区分助剂 #

这种方式保留trat作为x轴，通过颜色区分nozzle，填充区分adj，让每个处理下的所有组合柱子并排展示：

# 先加载所需包
library(ggplot2)
library(Rmisc)

# 你的数据和统计汇总（保持不变）
data <- structure(list(nozzle = c("XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "XR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR", "AIXR"), trat = c("Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 34", "Cle 34", "Cle 34", "Cle 34", "Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 34", "Cle 34", "Cle 34", "Cle 34", "Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 34", "Cle 34", "Cle 34", "Cle 34", "Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 12.8", "Cle 34", "Cle 34", "Cle 34", "Cle 34"), adj = c("Without", "Without", "Without", "Without", "Without", "Without", "Without", "Without", "With", "With", "With", "With", "With", "With", "With", "With", "Without", "Without", "Without", "Without", "Without", "Without", "Without", "Without", "With", "With", "With", "With", "With", "With", "With", "With"), dw1 = c(3.71, 5.87, 6.74, 1.65, 0.27, 0.4, 0.37, 0.34, 0.24, 0.28, 0.32, 0.38, 0.39, 0.36, 0.32, 0.28, 8.24, 10.18, 11.59, 6.18, 0.2, 0.23, 0.2, 0.31, 0.28, 0.25, 0.36, 0.27, 0.36, 0.37, 0.34, 0.19)), row.names = c(NA, -32L), class = c("tbl_df", "tbl", "data.frame"))
data_sum <- summarySE(data, measurevar="dw1", groupvars=c("nozzle", "trat","adj"))

# 修正后的绘图代码
ggplot(data_sum, aes(x = trat, y = dw1, 
                     group = interaction(nozzle, adj),  # 按喷嘴+助剂的组合分组
                     fill = adj,                        # 填充区分助剂
                     color = nozzle)) +                 # 颜色区分喷嘴
  geom_bar(stat = "identity", width = 0.7, position = position_dodge(width = 0.8)) +
  geom_errorbar(aes(ymax = dw1 + se, ymin = dw1 - se), 
                position = position_dodge(width = 0.8),  # 和柱状图统一dodge宽度，避免错位
                width = 0.2, color = "black") +
  scale_fill_manual(name = "Presence of adjuvants", 
                    values = c('grey80', 'grey30'), 
                    labels = c("Without", "With")) +
  scale_color_manual(name = "Nozzle Type", 
                     values = c("black", "darkred")) +
  labs(x = "Treatment g.i.a. ha⁻¹", y = "Dry Mass (g)") +
  theme_classic()

关键调整点： #

• 消除重叠：通过interaction(nozzle, adj)设置分组，确保每个独立组合都有自己的柱子；同时让geom_bar和geom_errorbar使用相同的position_dodge宽度，保证对齐且不重叠。
• 统一Y轴：移除了facet_grid，把两个因素直接整合到同一坐标系中。
• 可读性优化：调整了柱子和误差线的宽度，让图表更清晰；修正了x轴标签的上标格式。

方案二：将处理和喷嘴组合作为x轴 #

如果更偏好紧凑的布局，可以把trat和nozzle的组合作为x轴，用填充区分adj：

ggplot(data_sum, aes(x = interaction(trat, nozzle), y = dw1, 
                     fill = adj)) +
  geom_bar(stat = "identity", width = 0.7, position = position_dodge(width = 0.8)) +
  geom_errorbar(aes(ymax = dw1 + se, ymin = dw1 - se), 
                position = position_dodge(width = 0.8), 
                width = 0.2, color = "black") +
  scale_fill_manual(name = "Presence of adjuvants", 
                    values = c('grey80', 'grey30'), 
                    labels = c("Without", "With")) +
  labs(x = "Treatment + Nozzle", y = "Dry Mass (g)") +
  theme_classic() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))  # 旋转x轴标签避免重叠

这种方式更直观地展示每个处理+喷嘴组合下的助剂效果，同样不会出现柱子重叠的问题。

内容的提问来源于stack exchange，提问作者Vitor Muller Anunciato