将线性模型方程添加到绘图中

Question

我想用ggplot2将线性模型的方程式添加到绘图中。这个问题已经被问过很多次了，但是我的以下代码有问题-以前没有出现过问题-所以我认为这不是重复的问题。

我从stackoverflow的post中多次使用了以下代码

library(ggplot2)
df <- data.frame(x = c(1:100))
df$y <- 2 + 3 * df$x + rnorm(100, sd = 40)
p <- ggplot(data = df, aes(x = x, y = y)) +
  geom_smooth(method = "lm", se=FALSE, color="black", formula = y ~ x) +
  geom_point()
p

lm_eqn = function(m) {

  l <- list(a = format(coef(m)[1], digits = 2),
            b = format(abs(coef(m)[2]), digits = 2),
            r2 = format(summary(m)$r.squared, digits = 3));

  if (coef(m)[2] >= 0)  {
    eq <- substitute(italic(y) == a + b %.% italic(x)*","~~italic(r)^2~"="~r2,l)
  } else {
    eq <- substitute(italic(y) == a - b %.% italic(x)*","~~italic(r)^2~"="~r2,l)    
  }

  as.character(as.expression(eq));                 
}
p1 = p + geom_text(aes(x = 25, y = 300, label = lm_eqn(lm(y ~ x, df))), parse = TRUE)
p1

问题在于，在图中，方程的系数被c（）包围：

Answer 1

似乎代码无法正确地从系数向量中删除名称。尝试改用此版本

lm_eqn <- function(m) {

  l <- list(a = format(unname(coef(m)[1]), digits = 2),
            b = format(unname(abs(coef(m)[2])), digits = 2),
            r2 = format(summary(m)$r.squared, digits = 3));

  if (coef(m)[2] >= 0)  {
    eq <- substitute(italic(y) == a + b %.% italic(x)*","~~italic(r)^2~"="~r2,l)
  } else {
    eq <- substitute(italic(y) == a - b %.% italic(x)*","~~italic(r)^2~"="~r2,l)    
  }

  as.character(as.expression(eq));                 
}

我刚刚加入了几个unname()通话。