解决ggarrange合并子图时公共图例缺失全类别问题

核心思路 #

问题根源是每个子图的Mineral因子水平不完整，导致颜色映射和图例范围不一致。解决分两步：统一所有子图的颜色映射规则，生成包含全部Mineral类别的独立图例，再和子图组合。

具体步骤 #

1. 统一Mineral的因子水平 #

先将Mineral列转换为包含所有9个取值的因子，确保所有子图使用相同的颜色映射逻辑：

# 提取所有Mineral的唯一取值（或直接手动指定9个值的向量）
all_minerals <- unique(results$Mineral)
# 转换为因子，固定水平顺序
results$Mineral <- factor(results$Mineral, levels = all_minerals)

2. 生成完整的独立图例 #

绘制一个仅包含完整图例的空图，再提取图例对象：

library(ggpubr)

# 生成带完整图例的空图
full_legend_plot <- ggplot(results, aes(x = B, y = N, color = Mineral)) +
  geom_point(alpha = 0) + # 透明点，仅保留图例
  theme(legend.position = "bottom") + # 可调整图例位置
  guides(color = guide_legend(nrow = 2)) # 可选：调整图例行数，避免过宽

# 提取图例对象
full_legend <- get_legend(full_legend_plot)

3. 生成子图并组合 #

批量生成子图时关闭各自的图例，再用ggarrange组合子图和独立图例：

# 按mount分组生成子图列表
plot_list <- lapply(unique(results$mount), function(current_mount) {
  ggplot(results[results$mount == current_mount, ], aes(x = B, y = N, color = Mineral)) +
    geom_point() +
    labs(title = paste("Mount:", current_mount)) +
    theme(legend.position = "none") # 关闭子图自身图例
})

# 组合16张子图（示例为4行4列）
subplots_combined <- ggarrange(plotlist = plot_list, nrow = 4, ncol = 4)

# 将子图和图例组合，调整高度比例让图例更紧凑
final_plot <- ggarrange(subplots_combined, full_legend, nrow = 2, heights = c(10, 1))

# 输出最终图
print(final_plot)

关键优势 #

• 所有子图的Mineral颜色完全一致，不会因为子图子集缺失而改变颜色映射
• 公共图例包含全部9个Mineral类别，不会仅显示第一个子图的子集
• 基于ggpubr工具链，和原有的ggarrange使用逻辑一致，无需切换其他包

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