如何修改ggplot面板图中指定区域直方图的绘制层级

问题说明 #

制作了对比捕捞（Harvest）与非捕捞区域牡蛎壳高（Shell Height）频率的ggplot面板直方图，其中Matanzas River面板内的小直方图被较大的直方图遮挡，需要仅调整该面板的绘制顺序，其余面板保持不变，且不更改原有配色方案。

解决方案 #

核心是针对Matanzas River区域单独反转Harvest变量的绘制层级，其他区域保持原有顺序。通过新增一个分组变量，在该区域反转Harvest的水平，再用这个变量映射填充和颜色，同时保留原配色。

修改后的完整代码 #

# 加载必要包
library(ggplot2)
library(dplyr)

# 读取数据（假设Sh.csv已在工作目录中）
oyster_data <- read.csv("Sh.csv")

# 定义原配色
harvest_colors = c("Yes" = 'black', "No" = 'gray49')

# 新增分组变量，仅Matanzas River区域反转Harvest绘制顺序
oyster_data <- oyster_data %>%
  mutate(Harvest_order = ifelse(Region == "Matanzas River", 
                                factor(Harvest, levels = c("No", "Yes")),
                                factor(Harvest, levels = c("Yes", "No"))))

# 绘制直方图
ggplot(oyster_data, aes(x=Shell, fill=Harvest_order, color = Harvest_order)) +
  geom_histogram(alpha=0.5, position = 'identity') + 
  # 保持原配色与图例显示一致
  scale_fill_manual(values = harvest_colors, breaks = c("Yes", "No"), labels = c("Yes", "No")) + 
  scale_color_manual(values = harvest_colors, breaks = c("Yes", "No"), labels = c("Yes", "No")) +
  labs(x = "Shell Height (mm)", y = "Frequency") + 
  facet_wrap(~Region) +
  theme(panel.background = element_blank(), axis.line=element_line(size=1), 
        plot.title = element_text(hjust = 0.5)) +
  geom_vline(xintercept = 75, linetype = "dashed") + 
  ggtitle("Oyster Measurement Frequency by Region")

代码解释 #

• 新增Harvest_order变量：通过dplyr::mutate创建，当区域为Matanzas River时，将Harvest的水平设为c("No", "Yes")（让非捕捞组先绘制，捕捞组后绘制，避免小直方图被遮挡），其他区域保持原顺序。
• 映射新变量：将fill和color映射到Harvest_order，而非原Harvest变量，实现单独调整目标面板的绘制层级。
• 保持配色与图例一致性：在scale_fill_manual和scale_color_manual中设置breaks和labels，确保图例显示和原配色完全一致，不会因变量顺序改变产生混乱。

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