如何实现ggplot facet_wrap固定列宽并按高度重排堆叠柱状图

针对你提出的两个可视化优化目标，以下是具体的解决方案和代码修改：

1. 确保所有分面的柱状图列宽一致 #

默认使用facet_wrap(~Predictor, scales = "free")时，每个分面的宽度固定，但x轴离散变量数量不同会导致柱子宽度不一致。要解决这个问题，可以借助ggh4x包的facet_wrap2()函数，通过space = "free_x"参数让分面宽度与x轴类别数量成正比，从而保证柱子宽度统一。

首先安装并加载ggh4x包：

install.packages("ggh4x")
library(ggh4x)

2. 根据堆叠高度重排柱状图 #

要按堆叠后的总高度（即每个Variable的RI总和）对柱子排序，需要先计算每个Variable的RI总和，再将Variable转换为按总和排序的因子：

# 计算每个Variable的RI总和，按总和排序
RI_ph_sorted <- RI_pH %>%
  group_by(Variable) %>%
  mutate(total_RI = sum(RI)) %>%
  ungroup() %>%
  arrange(total_RI) %>%
  mutate(Variable = factor(Variable, levels = unique(Variable)))

修改后的完整代码 #

library(ggplot2)
library(dplyr)
library(ggh4x)

# 数据预处理：修正因子水平
RI_pH <- data.frame(RI_pH)
RI_pH$Landcover <- factor(RI_pH$Landcover, levels = c("Bare", "Vegetation", "Mixed"))
RI_pH$Predictor <- factor(RI_pH$Predictor, levels = c("Landsat8", "Landsat8 & Sentinel1",
                                                      "Sentinel1", "All Variables"))

# 按堆叠总高度排序Variable
RI_pH_sorted <- RI_pH %>%
  group_by(Variable) %>%
  mutate(total_RI = sum(RI)) %>%
  ungroup() %>%
  arrange(total_RI) %>%
  mutate(Variable = factor(Variable, levels = unique(Variable)))

# 绘制分面柱状图，保证列宽一致
bar <- ggplot(RI_pH_sorted, aes(x = Variable, y = RI, fill = Landcover)) +
  geom_col(width = 0.5) +
  # 使用facet_wrap2实现等宽柱子
  facet_wrap2(~Predictor, scales = "free_x", space = "free_x") +
  theme(axis.text.x = element_text(angle = 90, hjust = 0.95, vjust = 0.2, size = 32),
        axis.text.y = element_text(size = 32),
        strip.text = element_text(size = 52, face = "bold")) # 修正原代码的拼写错误

# 导出图片
jpeg(file = "bar%03d.jpeg", width = 1180, height = 849,
     pointsize = 52, bg = "white", units = "mm", res = 95)
plot(bar)
dev.off()

关键说明 #

• space = "free_x"让分面宽度随x轴类别数量自动调整，确保每个柱子的物理宽度完全一致；
• 通过计算total_RI并重新设置Variable的因子水平，实现了按堆叠总高度排序柱子，便于视觉对比不同Variable的整体大小；
• 修正了原代码中element text的拼写错误，改为标准的element_text。

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