线性拟合

画图+拟合+拟合线

plot(x = wp$Europe,y = wp$Asia,pch=19)

m <- lm(wp$Asia~wp$Europe)

abline(m,col='blue',lty=1)

coef(m)

summary(m)$r.squared

加标注

legend('bottomright',

      legend = c('Data','Linear fit'),

      pch = c(19,NA),

      lty = c(NA,1),

      col = c('black','blue'),

      bty = 'n')

加公式

eq <- substitute(italic(y) == b+a*italic(x)*","~italic(R)^2==r2,

                list(b = format(coef(m)[1], digits = 2),

                      a = format(coef(m)[2], digits = 2),

                      r2 = format(summary(m)$r.squared, digits = 3)))

text(x=30000,y = 8500,labels = as.expression(eq),col = 'blue')