R语言 向分面数据添加公式与拟合曲线

数据

源代码：

library(ggplot2)
data<-read.csv("D:\\rwork\\tiaocan_xiugai.csv",header = T)  #header = Falese的意思
dt<-as.data.frame(data)
dt2<-dt


#函数
lm_labels<-function(dat){
  count1<-dat$count1
  degree<-dat$degree
  mod<-nls(count1 ~ a*degree^b, start = list(a=2, b= 1.5))
  formula<-sprintf("italic(y)==%.3f*italic(x)^%.3f",
                   round(coef(mod)[1],3),round(coef(mod)[2],3))
  a=0
  b=0
  for(i in 1:15){
    ymean=mean(count1)
    a=a+(fitted(mod)[i]-ymean)^2
    b=b+(count1[i]-ymean)^2
  }
  r=a/b
  #r<-cor(dat$count1,dat$degree)
  r2<-sprintf("italic(R^2)==%.3f",r)