如何在GGally::ggpairs绘图中添加置信区间数值？

当然可以实现！你的思路方向完全正确——通过自定义上三角的continuous函数扩展ggpairs，但原代码里的注释位置问题，主要是因为用了固定坐标x=0,y=0，这会在数据范围不包含0时导致文本错位或消失。下面是优化后的解决方案，能完美实现类似参考图的效果：

改进后的自定义函数 #

这个函数会自动计算Pearson相关系数和95%置信区间，并且自适应每个单元格调整文本位置：

library(GGally)
library(rlang)

my_high <- function(data, mapping, method = "pearson", conf.level = 0.95, ...) {
  # 提取映射的x和y变量（避免硬编码列索引，适配任意数据组合）
  x_var <- rlang::get_expr(mapping$x)
  y_var <- rlang::get_expr(mapping$y)
  x <- data[[x_var]]
  y <- data[[y_var]]
  
  # 一次性计算相关系数和置信区间，提升效率
  cor_result <- cor.test(x, y, method = method, conf.level = conf.level)
  cor_val <- round(cor_result$estimate, 3)
  ci_lower <- round(cor_result$conf.int[1], 3)
  ci_upper <- round(cor_result$conf.int[2], 3)
  
  # 构建显示文本：分行展示相关系数和置信区间，提升可读性
  display_text <- paste0("r = ", cor_val, "\n95% CI: [", ci_lower, ", ", ci_upper, "]")
  
  # 基于ggally_cor的基础样式，添加自定义文本
  p <- ggally_cor(data, mapping = mapping, ...) +
    geom_text(
      aes(x = Inf, y = Inf),
      label = display_text,
      hjust = 1.1, vjust = 1.1, # 把文本固定在单元格右上角，避免被截断
      size = 3 # 可根据需求调整字体大小
    )
  
  return(p)
}

使用示例 #

以mtcars数据集的连续变量为例测试效果：

GGally::ggpairs(
  mtcars[, c("mpg", "disp", "hp", "wt")],
  upper = list(continuous = my_high)
)

关键优化点 #

1. 避免硬编码列索引：用rlang提取映射变量，适配任意数据列的组合，更符合ggplot的映射逻辑。
2. 高效计算：只调用一次cor.test，同时获取相关系数和置信区间，避免重复计算浪费资源。
3. 自适应位置：用x=Inf, y=Inf配合hjust/vjust，让文本固定在每个单元格的右上角，不会因为数据范围变化而错位。
4. 清晰排版：用换行符\n把相关系数和置信区间分行展示，和参考图的可读性一致。

更进一步：完全自定义中心文本 #

如果你想让文本居中显示（更贴近参考图的样式），可以直接创建空的ggplot图层，完全控制内容：

my_high_center <- function(data, mapping, method = "pearson", conf.level = 0.95, ...) {
  x <- data[[rlang::get_expr(mapping$x)]]
  y <- data[[rlang::get_expr(mapping$y)]]
  
  cor_result <- cor.test(x, y, method = method, conf.level = conf.level)
  cor_val <- round(cor_result$estimate, 3)
  ci_lower <- round(cor_result$conf.int[1], 3)
  ci_upper <- round(cor_result$conf.int[2], 3)
  
  display_text <- paste0("r = ", cor_val, "\n[", ci_lower, ", ", ci_upper, "]")
  
  # 空ggplot+居中文本+极简主题，复刻参考图核心样式
  p <- ggplot(data, mapping) +
    geom_blank() +
    geom_text(aes(x = mean(x), y = mean(y)), label = display_text, size = 3.5) +
    theme_void()
  
  return(p)
}

# 使用这个版本
GGally::ggpairs(
  mtcars[, c("mpg", "disp", "hp", "wt")],
  upper = list(continuous = my_high_center)
)

这个版本会把相关系数和置信区间放在每个单元格的正中心，完全匹配你想要的参考图效果。

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