ggpairs对角线绘制带正态曲线的直方图报错求助

这个错误的根源在于ggpairs在调用自定义对角线函数时，会自动传递额外的参数（比如mapping），但你的myhist函数只定义了data参数，完全没处理这些自动传入的参数，所以才会抛出unused argument的报错。咱们来一步步修复它：

1. 理解ggpairs自定义函数的参数要求 #

GGally的ggpairs在调用自定义对角线函数时，会默认传入data（当前列的数据集）、mapping（x轴的映射关系，比如~x1）以及其他可选参数。你的自定义函数必须能接收这些参数，否则就会报错。

2. 修改自定义直方图函数 #

调整函数定义，让它接收mapping和额外参数，同时从映射中提取当前变量来计算均值和标准差：

library(GGally)
library(ggplot2)
library(rlang) # 用于处理表达式映射

myhist <- function(data, mapping, ...) {
  # 从mapping中提取当前要绘制的变量名
  x_var <- rlang::get_expr(mapping$x)
  # 计算该变量的均值和标准差（加入na.rm处理缺失值）
  x_mean <- mean(data[[x_var]], na.rm = TRUE)
  x_sd <- sd(data[[x_var]], na.rm = TRUE)
  
  ggplot(data = data, mapping = mapping) +
    # 直方图用密度作为y轴，保留默认的binwidth或自定义传入
    geom_histogram(aes(y = ..density..), colour = "black", fill = "white", ...) +
    # 叠加正态密度曲线，用当前变量的均值和标准差
    stat_function(fun = dnorm, args = list(mean = x_mean, sd = x_sd), colour = "red", size = 1)
}

3. 测试完整代码 #

用你的示例数据（或模拟数据）调用ggpairs：

# 模拟示例数据
sample.data <- data.frame(
  x1 = rnorm(100),
  x2 = rnorm(100),
  x3 = rnorm(100),
  x4 = rnorm(100),
  x5 = rnorm(100),
  x6 = rnorm(100),
  F1 = rnorm(100),
  F2 = rnorm(100)
)

# 调用ggpairs，指定自定义对角线函数
ggpairs(sample.data, diag = list(continuous = myhist))

额外说明 #

• 如果不想依赖rlang包，也可以用基础R的方式提取变量名：x_var <- as.character(mapping$x)[2]，不过rlang的处理方式更稳妥，能应对复杂的映射表达式。
• 加入na.rm = TRUE是为了让函数在数据有缺失值时也能正常运行，提升健壮性。

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