如何用ggplot2绘制保持柱形尺寸位置一致的分组子柱状图

首先给出数据框定义：

df1 <- data.frame(
  trait = c("CRP", "CRP", "CRP", "TotChol", "TotChol", "TotChol", "RBC", "RBC", "RBC", "glucose", "glucose", "glucose"),
  Variable = c("PRSeur_r2_PGScatalog", "PRSeur_r2_mol", "moliPRS_eurPRS", 
               "PRSeur_r2_PGScatalog", "PRSeur_r2_mol", "moliPRS_eurPRS", 
               "PRSeur_r2_PGScatalog", "PRSeur_r2_mol", "moliPRS_eurPRS", 
               "PRSeur_r2_PGScatalog", "PRSeur_r2_mol", "moliPRS_eurPRS"),
  Value = c(14, 3.5, 3.9, 11, 9, 9.4, 15, 7, 8.11, 3, 1.3, 1.4)
)

假设你原整体分组柱状图的代码类似这样（核心参数需保留）：

library(ggplot2)

# 原整体分组柱状图
p_original <- ggplot(df1, aes(x = trait, y = Value, fill = Variable)) +
  geom_col(position = position_dodge(width = 0.8), width = 0.7) +
  theme_classic()

要生成每个trait的独立子图且柱形尺寸、位置与原图完全一致，关键是保留绘图核心参数并固定坐标轴范围，以下两种方法可选：

方法一：用facet_wrap拆分并固定轴范围 #

这种方法快速生成子图面板，所有子图共享绘图参数：

p_subplots <- ggplot(df1, aes(x = trait, y = Value, fill = Variable)) +
  # 完全复用原图的柱形宽度和分组间距参数
  geom_col(position = position_dodge(width = 0.8), width = 0.7) +
  # 按trait拆分面板
  facet_wrap(~trait) +
  # 固定x轴范围为所有trait，保证柱形位置与原图对齐
  scale_x_discrete(limits = unique(df1$trait)) +
  # 可选：固定y轴范围，让所有子图y轴统一
  scale_y_continuous(limits = c(0, max(df1$Value) * 1.1)) +
  theme_classic() +
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1),
    strip.text = element_text(face = "bold")
  )

p_subplots

方法二：批量生成独立子图（适合单独导出） #

如果需要每个trait对应独立的图片文件，用purrr批量生成：

library(purrr)
library(gridExtra)

# 定义复用原图参数的绘图函数
plot_single <- function(target_trait) {
  df_sub <- df1[df1$trait == target_trait, ]
  
  ggplot(df_sub, aes(x = trait, y = Value, fill = Variable)) +
    geom_col(position = position_dodge(width = 0.8), width = 0.7) +
    # 固定x轴范围，对齐原图位置
    scale_x_discrete(limits = unique(df1$trait)) +
    scale_y_continuous(limits = c(0, max(df1$Value) * 1.1)) +
    theme_classic() +
    labs(title = target_trait) +
    theme(
      plot.title = element_text(hjust = 0.5),
      # 可选：隐藏x轴多余刻度，只保留当前trait的显示
      axis.text.x = element_text(color = ifelse(unique(df_sub$trait) == unique(df1$trait), "black", "white"))
    )
}

# 生成所有trait的子图列表
trait_list <- unique(df1$trait)
subplot_list <- map(trait_list, plot_single)

# 排列显示所有子图
grid.arrange(grobs = subplot_list, ncol = 2)

# 可选：批量导出子图
# walk2(subplot_list, trait_list, ~ggsave(paste0(.y, ".png"), plot = .x, width = 6, height = 4))

核心注意点 #

• 必须完全复用原整体图中geom_col()的width和position_dodge()的width参数，这是保证柱形尺寸一致的关键
• 固定x轴范围为所有trait的集合，确保每个子图的柱形在x轴上的位置与原图对应trait的柱形位置完全对齐
• 统一y轴范围可让子图视觉风格更一致，按需选择

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