如何在函数内为ggplot的stat_summary传递参数并添加观测数注释

你的问题大概率出在以下几个ggplot2自定义stat_summary函数的常见误区里，结合暴露-反应分析的分面箱线图场景，逐一拆解：

1. 自定义give.n函数的返回格式不符合要求 #

stat_summary要求自定义函数必须返回包含x、y、label列的data.frame，很多人会错误地直接打印文本或返回非数据框格式，导致无法生成注释。另外，需要确保函数能正确识别分面和分组变量，计算每个分面内各x类别下的观测数。

错误示例（常见问题）：

give.n <- function(x){
  cat("n =", length(x)) # 错误：用cat打印而非返回数据框
}

正确的give.n函数写法：

give.n <- function(x) {
  data.frame(
    y = max(x, na.rm = TRUE) * 1.1, # 把注释放在箱线图上方10%的位置
    label = paste0("n = ", length(na.omit(x))) # 统计非缺失值的观测数
  )
}

2. stat_summary的参数配置错误 #

• 必须指定geom = "text"，明确告诉ggplot要生成文本注释；
• 颜色设置color = "blue"要放在stat_summary的参数中（而非aes映射内），否则会被当作变量映射触发报错；
• 分面场景下，无需额外指定group参数，ggplot会自动传递分面的分组逻辑给stat_summary。

错误示例（常见问题）：

stat_summary(fun.data = give.n, aes(color = "blue")) # 错误：颜色放在aes映射里

正确的stat_summary调用：

stat_summary(fun.data = give.n, geom = "text", color = "blue", size = 3)

3. 分面变量的分组逻辑未被正确识别 #

如果你的分面变量是在facet_wrap或facet_grid中设置的，要确保数据中的分面变量是因子类型，give.n函数会自动按「分面+ x轴类别」的组合分组计算观测数（ggplot会自动传递分组信息）。

完整修复后的示例代码 #

构造示例患者暴露-反应数据 #

set.seed(123)
patient_data <- data.frame(
  exposure = rnorm(200, mean = 50, sd = 15),
  response = factor(sample(c("缓解", "稳定", "进展"), 200, replace = TRUE)),
  category = factor(sample(c("A组", "B组", "C组"), 200, replace = TRUE))
)

自定义绘图函数 #

plot_exposure_response <- function(data) {
  ggplot(data, aes(x = response, y = exposure)) +
    geom_boxplot(fill = "lightgray", color = "black") +
    # 正确添加蓝色观测数注释
    stat_summary(
      fun.data = give.n,
      geom = "text",
      color = "blue",
      size = 3
    ) +
    facet_wrap(~ category) +
    labs(x = "响应分类", y = "暴露值", title = "患者暴露-反应分面箱线图") +
    theme_bw()
}

# 定义符合要求的give.n函数
give.n <- function(x) {
  data.frame(
    y = max(x, na.rm = TRUE) * 1.1,
    label = paste0("n = ", length(na.omit(x)))
  )
}

# 运行函数
plot_exposure_response(patient_data)

额外注意事项 #

• 如果x轴变量是数值型，需先转换为因子factor(response)，否则stat_summary会按连续变量处理，无法正确计算每个分类的观测数；
• 调整注释位置可修改y = max(x, na.rm = TRUE) * 1.1中的倍数，比如1.2会把注释放得更高；
• 存在缺失值时，必须用na.omit(x)统计有效观测数，避免结果错误。

内容的提问来源于stack exchange，提问作者Ksenia Kasey Arzumanova