如何在ggplot2中随线性模型线条宽度同步调整置信带宽度？

当你增大stat_smooth的线条size后，默认的置信带（基于模型统计置信区间）会被粗线条遮挡。要让置信带边缘避开粗线条，需要手动调整置信区间的上下限，把线条半宽对应的数值偏移量加入计算。

方法一：手动计算并绘制 #

先拟合模型生成预测数据，再根据线条厚度调整置信区间：

library(ggplot2)
library(dplyr)

# 1. 拟合线性模型并生成默认预测置信区间
lm_model <- lm(Sepal.Length ~ Petal.Width, data = iris)
pred_data <- predict(lm_model, 
                     newdata = data.frame(Petal.Width = seq(min(iris$Petal.Width), max(iris$Petal.Width), length.out = 100)),
                     interval = "confidence") %>%
  as.data.frame() %>%
  mutate(Petal.Width = seq(min(iris$Petal.Width), max(iris$Petal.Width), length.out = 100))

# 2. 计算线条厚度对应的y轴数据偏移量
line_size <- 5  # 你设置的线条size
y_range <- diff(range(iris$Sepal.Length))
plot_height_mm <- 177.8  # 默认ggplot绘图高度7英寸=177.8mm，可按需调整
mm_to_y_unit <- y_range / plot_height_mm
offset <- (line_size / 2) * mm_to_y_unit  # 线条半宽对应的y轴偏移

# 3. 调整置信区间上下限
pred_data <- pred_data %>%
  mutate(lwr_adj = lwr - offset,
         upr_adj = upr + offset)

# 4. 绘图：先画置信带，再画粗线条
ggplot(iris, aes(x = Petal.Width, y = Sepal.Length)) +
  geom_point() +
  geom_ribbon(data = pred_data, aes(x = Petal.Width, ymin = lwr_adj, ymax = upr_adj),
              fill = "red", alpha = 0.2) +
  geom_line(data = pred_data, aes(y = fit), color = "red", size = line_size)

方法二：自定义stat实现复用 #

如果需要多次使用，可以自定义一个stat，自动根据线条size调整置信带：

library(ggplot2)

# 自定义Stat
StatSmoothAdj <- ggproto("StatSmoothAdj", StatSmooth,
                         compute_group = function(data, scales, ...) {
                           result <- super$compute_group(data, scales, ...)
                           # 获取线条size，默认1
                           line_size <- data$size[1] %||% 1
                           # 计算y轴mm到数据单位的转换
                           y_range <- diff(scales$y$range$range)
                           plot_height_mm <- 177.8
                           mm_to_y_unit <- y_range / plot_height_mm
                           offset <- (line_size / 2) * mm_to_y_unit
                           # 调整置信带上下限
                           result$ymin <- result$ymin - offset
                           result$ymax <- result$ymax + offset
                           result
                         })

# 封装成可调用的函数
stat_smooth_adj <- function(mapping = NULL, data = NULL, geom = "smooth",
                            position = "identity", ..., method = "auto",
                            formula = y ~ x, se = TRUE, na.rm = FALSE,
                            show.legend = NA, inherit.aes = TRUE) {
  layer(
    data = data,
    mapping = mapping,
    stat = StatSmoothAdj,
    geom = geom,
    position = position,
    show.legend = show.legend,
    inherit.aes = inherit.aes,
    params = list(
      method = method,
      formula = formula,
      se = se,
      na.rm = na.rm,
      ...
    )
  )
}

# 使用示例
ggplot(iris, aes(x = Petal.Width, y = Sepal.Length)) +
  geom_point() +
  stat_smooth_adj(method = "lm", col = "red", size = 5, fill = "red", alpha = 0.2)

注意事项 #

• 偏移量的计算依赖于绘图区域的实际高度，如果修改了图的尺寸（比如用ggsave指定高度），需要同步调整plot_height_mm的值。
• 这种方法保证线条半宽在y轴方向上刚好对齐置信带边缘，解决粗线条遮挡置信带的问题。

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