如何让ggplot2绘制的IQR图中中位数标识横线铺满对应条的宽度

你现有代码里用shape="_"的point层绘制中位数线的方式，无法动态匹配不同分组条的宽度，所以会出现横线过短的问题，可参考以下两种方案解决：

方案1：复用现有crossbar图层（无需额外安装包） #

crossbar geom本身自带中位数横线，只需要调整参数即可实现宽度自动匹配：

library(ggplot2)
library(dplyr)

# 自定义中位数+IQR统计量函数
median_IQR <- function(x) {
  data.frame(y = median(x), # 中位数
             ymin = quantile(x)[2], # 下四分位数
             ymax = quantile(x)[4])  # 上四分位数
}

iris <- mutate(iris, let = rep(c("a","b"), nrow(iris)/2))

ggplot(iris, aes(x=Species, y=Sepal.Length, color = let, fill=let )) +
  stat_summary(geom = "crossbar",
               fun.data = median_IQR, 
               position = position_dodge(width = 0.5),
               width = 0.4, # 控制每个IQR条的总宽度
               fatten = 3) # 控制中间中位数横线的粗细

方案2：使用ggdist包实现更灵活的自定义 #

如果需要更多样式调整，可使用专门做分布区间可视化的ggdist包，同时支持原始数据和汇总数据输入：

library(ggplot2)
library(dplyr)
library(ggdist)

iris <- mutate(iris, let = rep(c("a","b"), nrow(iris)/2))

ggplot(iris, aes(x=Species, y=Sepal.Length, color = let, fill=let )) +
  # 绘制IQR区间
  stat_interval(
    .width = 0.5, # 0.5分位区间即IQR
    position = position_dodge(width = 0.5),
    size = 8 # 调整IQR条的粗细
  ) +
  # 绘制中位数横线，宽度和IQR条自动匹配
  stat_summary(
    fun = median,
    geom = "errorbarh",
    aes(xmin = after_stat(x) - 0.2, xmax = after_stat(x) + 0.2),
    position = position_dodge(width = 0.5),
    linewidth = 1,
    color = "black"
  )

两种方案最终生成的图中，中位数横线都会和对应IQR条的宽度完全一致，且不会展示最小值、最大值和离群点，符合敏感数据的展示要求。

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