ggplot2點圖+線性趨勢+公式+R2+p值
原創weixin_43948357 最后發布於2020-04-06 00:28:58 閱讀數 13 收藏
展開
正文
先看效果。
R語言代碼如下，
data("faithful")
library(ggplot2)
p <- ggplot(faithful,aes(x=eruptions,y=waiting)) + geom_point() + stat_smooth(method='lm',formula = y~x,colour='red')
model.lm<-lm(formula = waiting ~ eruptions, data = faithful)
#summary(model.lm)
#對於一元線性回歸方程y=ax+b，Intercept是指的截距，x對應的是系數。
l <- list(a = as.numeric(format(coef(model.lm)[1], digits = 4)),
b = as.numeric(format(coef(model.lm)[2], digits = 4)),
r2 = format(summary(model.lm)$r.squared, digits = 4),
p = format(summary(model.lm)$coefficients[2,4], digits = 4))
eq <- substitute(italic(y) == a + b %.% italic(x)~","~italic(R)^2~"="~r2~","~italic(P)~"="~p, l)
p + geom_text(aes(x = 4, y = 50, label = as.character(as.expression(eq))), parse = TRUE)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
參考文獻
【1】https://www.jianshu.com/p/d6a5eb41ecec
【2】https://www.zhihu.com/question/54776335/answer/141091529
————————————————
版權聲明：本文為CSDN博主「weixin_43948357」的原創文章，遵循 CC 4.0 BY-SA 版權協議，轉載請附上原文出處鏈接及本聲明。
原文鏈接：https://blog.csdn.net/weixin_43948357/java/article/details/105336901