如何在R语言中创建共享X轴的多面板直方图

嘿，我来帮你搞定这个共享X轴的多面板直方图问题！你现在用基础R的hist()函数逐个处理每个Site的数据，不仅繁琐还很难保证坐标轴对齐，其实用ggplot2的分面功能就能轻松实现需求，代码更简洁，输出的图也更规整。

步骤1：整合数据过滤逻辑 #

你之前针对每个Site做的分位数过滤异常值的操作，可以用分组处理一次性完成，不用重复写代码。这里我推荐用tidyverse工具集（包含ggplot2和dplyr）来简化数据处理：

# 加载所需包
library(tidyverse)

# 读取数据（保留你原来的读取逻辑）
ifile <- read.table("C:/files/text.lu", skip = 2, header = TRUE, sep="\t")
ifileVf2 <- data.frame(ifile["Vf2"], ifile["Site"])

# 按Site分组，过滤掉Vf2的0.001和0.999分位数之外的异常值
filtered_data <- ifileVf2 %>%
  group_by(Site) %>%
  filter(Vf2 > quantile(Vf2, 0.001) & Vf2 < quantile(Vf2, 0.999))

步骤2：绘制共享X轴的多面板直方图 #

用ggplot2的facet_wrap()或facet_grid()就能实现分面，默认就会共享坐标轴。如果你想让Y轴根据每个面板的频率自由调整（避免不同Site频率差异大导致图不好看），可以加上scales = "free_y"参数，X轴依然保持统一：

# 绘制纵向排列的多面板直方图（共享X轴）
ggplot(filtered_data, aes(x = Vf2)) +
  # 直方图参数：设置 bins 数量，填充色和边框色
  geom_histogram(bins = 150, fill = "steelblue", color = "white") +
  # 按Site分面，纵向排列（ncol=1），共享X轴，Y轴自由调整
  facet_wrap(~Site, ncol = 1, scales = "free_y") +
  # 使用简洁的主题
  theme_bw() +
  # 设置坐标轴标题
  labs(x = "Vf2 数值", y = "频数") +
  # 调整文字大小，提升可读性
  theme(
    strip.text = element_text(size = 12), # 面板标题文字大小
    axis.title = element_text(size = 14)  # 坐标轴标题文字大小
  )

可选调整：横向排列面板 #

如果想让两个面板横向并排展示，只需要把facet_wrap()里的ncol=1改成ncol=2即可：

facet_wrap(~Site, ncol = 2, scales = "free_y")

为什么不用基础R的hist()？ #

你原来的方式需要手动创建每个直方图，还要手动调整布局对齐坐标轴，不仅耗时还容易出错。ggplot2的分面功能会自动帮你完成坐标轴对齐、面板排版，而且样式调整更灵活，能快速生成符合学术或报告要求的图表。

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