如何在ggplot分面（facets）中添加plotmath格式的标签？

常规预构建标签再手动添加的思路依然适用，只需要额外适配plotmath的解析规则即可，具体实现步骤如下：

• 按分面分组计算回归参数，同时按照plotmath语法构造标签字符串：需要上标用^标记，下标用[]标记，等于符号用==，空格用~占位，特殊符号直接用对应plotmath语法书写即可。
• 生成独立的标签数据框：每一行对应一个分面，保留分面分组变量和构造完成的标签字符串，同时可以提前指定每个标签的x、y坐标，控制显示位置。
• 绘图时调用geom_text()或geom_label()加载标签数据框，额外设置parse = TRUE参数，ggplot会自动将字符串解析为plotmath格式渲染。

代码示例 #

library(ggplot2)
library(dplyr)

# 按分面变量cyl分组计算回归参数并构造标签
label_df <- mtcars %>% 
  group_by(cyl) %>% 
  do(lm_model = lm(mpg ~ wt, data = .)) %>% 
  mutate(
    intercept = round(coef(lm_model)[1], 2),
    slope = round(coef(lm_model)[2], 2),
    r2 = round(summary(lm_model)$r.squared, 2),
    # 构造带R²上标的plotmath格式标签
    label = paste0("y == ", intercept, " + ", slope, "*x ~~ R^2 == ", r2)
  )

# 绘图输出
ggplot(mtcars, aes(x = wt, y = mpg)) +
  geom_point() +
  geom_smooth(method = "lm", se = FALSE) +
  # 加载标签数据，开启parse解析
  geom_text(data = label_df, aes(x = 2.5, y = 30, label = label), parse = TRUE, hjust = 0) +
  facet_wrap(~cyl)

注意事项 #

• 标签字符串的语法必须严格符合plotmath规则，否则会出现解析报错，常见错误包括把==写为=、没有用~代替空格等。
• 如果需要渲染希腊字母、斜体等特殊格式，直接在标签字符串中添加对应语法即可，比如italic(beta) == 0.6会渲染为斜体的β值。

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