ggplot2：分面图注释中文本与数学表达式的结合方法

你遇到的核心问题是要在注释里同时结合文本、动态变量值和数学格式（上标单位），expression()因为是静态表达式没法嵌入变量，而**bquote()**就是专门解决这个场景的工具——它能把环境中的变量值和plotmath格式的表达式无缝结合。

下面用ggplot2和基础绘图两种常见场景给你演示具体实现：

场景1：ggplot2分面图 #

先准备数据并计算每个分面的均值，再用geom_text配合bquote()添加注释：

library(ggplot2)
library(dplyr)

# 1. 准备示例数据并计算各分面的均值
mtcars_mean <- mtcars %>%
  group_by(cyl) %>%
  summarise(avg_value = round(mean(mpg), 2)) # 这里替换成你的y值均值计算

# 2. 绘制分面图并添加带单位的注释
ggplot(mtcars, aes(x = wt, y = mpg)) +
  geom_point() +
  facet_wrap(~cyl) + # 你的分面变量替换这里
  # 用bquote构造注释文本：嵌入变量+上标单位
  geom_text(
    data = mtcars_mean,
    aes(x = 4, y = 30, # 注释的位置，根据你的图调整
        label = bquote(Avg == .(avg_value) ~ m^2/ha)),
    size = 4, color = "red"
  )

代码解释： #

• bquote(Avg == .(avg_value) ~ m^2/ha)：
  • .(avg_value)：把计算好的均值变量插入到表达式中
  • ~：表示空格分隔文本和单位
  • m^2：用plotmath语法实现上标效果，对应你需要的m²
  • 整个表达式会被解析成带动态值和正确单位格式的注释文本

场景2：基础绘图分面 #

如果用base R的分面（比如par(mfrow)），同样用bquote()配合text()函数：

# 1. 计算各分面的均值
mean_values <- tapply(mtcars$mpg, mtcars$cyl, function(x) round(mean(x), 2))

# 2. 绘制分面图并添加注释
par(mfrow = c(1, 3))
for(cyl in unique(mtcars$cyl)){
  sub_data <- mtcars[mtcars$cyl == cyl, ]
  plot(sub_data$wt, sub_data$mpg, main = paste("Cylinders:", cyl))
  # 添加带单位的注释
  text(
    x = 4, y = 30,
    labels = bquote(Avg == .(mean_values[as.character(cyl)]) ~ m^2/ha),
    col = "blue"
  )
}

为什么expression()不行？ #

expression()是静态的，它只会解析内部的plotmath语法，无法引用外部的变量值。比如你写expression(Avg == avg_value ~ m^2/ha)，avg_value不会被替换成实际的均值，只是作为文本显示。而bquote()通过.语法允许你把环境中的动态变量嵌入到表达式里，完美适配注释需要动态值的场景。

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