按Program D中位数重排ggplot分组箱线图的Y轴（Site）

核心思路 #

先单独计算每个Site对应的Program D的Result中位数（包含仅含D的Site），再用该中位数的降序结果重新指定Site因子的水平顺序，最后绘图时确保同Site的两个箱线图相邻。

具体实现步骤 #

1. 加载依赖包

library(dplyr)
library(ggplot2)

2. 计算每个Site的Program D中位数
   筛选出Program为D的记录，按Site分组计算中位数，确保仅含D的Site也能被纳入统计：

site_d_median <- TP %>%
  filter(Program == "D") %>%
  group_by(Site) %>%
  summarise(d_median = median(Result, na.rm = TRUE))

3. 重新排序Site因子
   将原数据的Site转换为有序因子，水平顺序按Program D中位数的降序排列：

TP$Site <- factor(TP$Site, 
                  levels = site_d_median %>% arrange(desc(d_median)) %>% pull(Site),
                  ordered = TRUE)

4. 绘制分组箱线图
   通过position_dodge参数保证同Site的D、H箱线图相邻：

ggplot(TP, aes(x = Result, y = Site)) +
  geom_boxplot(aes(fill = Program), outliers = FALSE, position = position_dodge(width = 0.75)) +
  theme_bw()

效果说明 #

• 示例数据中，Coal的Program D中位数为(0.65+0.58)/2=0.615，STP的为(0.55+0.53)/2=0.54，因此Coal会排在Y轴顶部，符合需求。
• 对于仅含Program D数据的Site，会自动按自身D组中位数参与排序，且仅显示对应箱线图。

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