R、R²、p 值和回归方程

R, R², p-value and regression equation

此代码为我提供了回归方程和 R2 的图:(但我需要提及方程将在其中的 x 和 y 中(手动)

CORRELATIONP3 <-CORRELATIONP2[product=='a',]

x<-CORRELATIONP3$b
y<-CORRELATIONP3$p


df <- data.frame(x = x)
m <- lm(y ~ x, data = df)
p <- ggplot(data = df, aes(x = x, y = y)) +
  scale_x_continuous("b (%)") +
  scale_y_continuous("p (%)")+
  geom_smooth(method = "lm", formula = y ~ x) +
  geom_point()
p

eq <- substitute(italic(y) == a + b %.% italic(x)*","~~italic(r)^2~"="~r2,
             list(        a = format(coef(m)[1], digits = 4),
                          b = format(coef(m)[2], digits = 4),
                          r2 = format(summary(m)$r.squared, digits = 3)))

dftext <- data.frame(x = 3, y = 0.2, eq = as.character(as.expression(eq)))

p + geom_text(aes(label = eq), data = dftext, parse = TRUE)

但是,通过这段代码,我得到了 R 和 p 值:这里关于 R 和 p 值的信息自动适合绘图,为什么?我也想要第一个。

  CORRELATIONP3 <-CORRELATIONP2[product=='a',]

x<-CORRELATIONP3$b
y<-CORRELATIONP3$p


df <- data.frame(x = x)
m <- lm(y ~ x, data = df)
p <- ggplot(data = df, aes(x = x, y = y)) +
  scale_x_continuous("b (%)") +
  scale_y_continuous("p (%)")+
  geom_smooth(method = "lm", formula = y ~ x) +
  geom_point()
p

eq <- substitute(italic(r)~"="~rvalue*","~italic(p)~"="~pvalue, list(rvalue = sprintf("%.2f",sign(coef(m)[2])*sqrt(summary(m)$r.squared)), pvalue = format(summary(m)$coefficients[2,4], digits = 3)))


dftext <- data.frame(x = 30, y = 0.4, eq = as.character(as.expression(eq)))
p + geom_text(aes(label = eq), data = dftext, parse = TRUE)

你能告诉我如何将所有 4 条信息连接到一个图中吗? (R、R2、方程和 p 值)

除此之外,我希望这些信息可以自动安装到绘图中,而不是手动安装。

好的,我不确定这是否有效,因为您没有提供可重现的数据示例,但我想您只需要重命名其中一个变量即可,例如:

eq2 <- substitute(italic(r)~"="~rvalue*","~italic(p)~"="~pvalue, 
list(rvalue = sprintf("%.2f",sign(coef(m)[2])*sqrt(summary(m)$r.squared)), 
pvalue = format(summary(m)$coefficients[2,4], digits = 3)))

然后你改变你把它放在你的地块上的点,就在第一个地块的另一个块下面一点。这里的 x 和 y 指的是文本标签的位置,因此请尝试使用它们直到您的文本看起来不错。

dftext2 <- data.frame(x = 30, y = 0.12, eq2 = as.character(as.expression(eq2)))
p + geom_text(aes(label = eq2), data = dftext2, parse = TRUE)

请让我知道这是否有效以及这是否是您的意思。