属于比较常用的技巧。下面我用ggplot2内置数据集进行演示：
先来看一下，仅有散点和拟合线的图
df=diamonds[sample(1:dim(diamonds)[1],40),]
df$price=df$price / 10000
p=df%>%ggplot(aes(carat,price))+
geom_point(size=4,alpha=0.3,color="#6baed6")+
geom_smooth(method = "lm", formula = y~x, color = "#756bb1", fill = "#cbc9e2")+ #颜色选自https://colorbrewer2.org/
theme_bw()+
theme(
panel.grid.major = element_blank(),panel.grid.minor = element_blank()
)
p
ggsave("tmp1.pdf",width = 8,height = 8,units = "cm")
如果想添加拟合的方程，可以使用ggpmisc包的stat_poly_eq()函数，该函数可以拟合多项式方程并生成多种标签，比如p值、决定系数R方，这些由函数计算出来的量通常名称前后有两个点
..eq.label..
..adj.rr.label..
..p.value.label..
对应的代码只有一行，其余代码和上面一样：
library(ggpmisc)
p+stat_poly_eq(aes(label=paste(..eq.label..,..adj.rr.label..,..p.value.label..,sep = "~~~~")),formula = y~x,parse=T,size=2.5)
#这一行paste()中的内容是一个类似Latex的表达，“~”表示空格，parse=T表示将这个语句翻译成可读形式；size=2.5表示字体大小
ggsave("tmp2.pdf",width = 8,height = 8,units = "cm")
如果需要添加相关系数，可以按照相同的思路，先求出对应的值，再加到stat_poly_eq()的label参数里
cor=round(cor(df$carat,df$price),2)
p+stat_poly_eq(aes(label=paste("italic(r)~`=`~",cor,sep = "")),formula = y~x,parse=T,size=4)
ggsave("tmp3.pdf",width = 8,height = 8,units = "cm")
另外也可以使用ggpubr包的stat_cor()函数，添加相关系数以及p值
library(ggpubr)
p+stat_cor()
ggsave("tmp4.pdf",width = 8,height = 8,units = "cm")
将stat_poly_eq()和stat_cor()结合起来也是可以的，当然两个文本的位置需要调整一下
p+stat_poly_eq(aes(label=..eq.label..),formula = y~x,parse=T,size=4)+
stat_cor(label.y=0.9,size=4)
因水平有限，有错误的地方，欢迎批评指正！
