基于ggplot2为分面行为数据图添加趋势线方程的技术咨询

要在每个分面子图中自动显示线性趋势线的方程（包含斜率作为相关指数），无需手动拆分数据，推荐使用ggpmisc包的stat_poly_eq函数，它能自动按分面分组拟合模型并生成标注。

步骤1：安装并加载依赖包 #

# 首次使用时执行安装
install.packages("ggpmisc")
# 加载所需包
library(ggplot2)
library(ggpmisc)

步骤2：修改原绘图代码 #

将stat_poly_eq层添加到你的ggplot代码中，配置参数以显示方程和相关统计量：

ggplot(data = water, mapping = aes(x = avg_time, y = time, color = as.factor(ID))) +
  facet_wrap(~action) +
  geom_point() +
  geom_smooth(method = lm, color = 'black') +
  # 添加回归方程标注
  stat_poly_eq(
    formula = y ~ x,  # 匹配线性回归模型
    aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")),
    parse = TRUE,  # 解析为规范的数学表达式
    color = "black",  # 标注颜色与趋势线统一
    size = 3  # 调整标注字体大小
  ) +
  labs(
    x = 'INDEX', 
    y = 'Time spent in each behavior (sec)', 
    title = 'FIGURE 12: A little bit of p hacking for the water trails', 
    color = 'Spider ID'
  )

参数说明 #

• formula = y ~ x：与geom_smooth(method = lm)使用的线性模型保持一致
• ..eq.label..：自动生成包含斜率、截距的回归方程
• ..rr.label..：自动生成R²值（可根据需求选择是否保留）
• parse = TRUE：将文本解析为LaTeX格式数学表达式，让方程显示更专业

无额外包的替代方案（不推荐） #

如果不想安装新包，可提前按action分组拟合线性模型，提取方程文本后用geom_text手动添加，但步骤繁琐且易出错，优先推荐ggpmisc的方法。

内容的提问来源于stack exchange，提问作者Katie Van Havel