#“line2P” # y=a*x+b

#“line3P” # y=a*x^2+b*x+c

#“log2P” # y=a*ln(x)+b

#“exp3P” # y=a*exp(b*x)+c

#“power3P” # y=a*x^b+c)

#install.packages("basicTrendline")

#Pan-genome,“power3P”

#Core-genome,"exp3P"

library(basicTrendline)

data1<-read.csv("pan.csv",header = TRUE)

x <- data1$x

y <- data1$y

p1<-trendline(x,y,model="power3P", ePos.x= "top")

data2<-read.csv("Core.csv",header = TRUE)

x <- data2$x

y <- data2$y

p2<-trendline(x,y,model="exp3P", ePos.x= "top")

参考来自：http://blog.sina.com.cn/s/blog_6ac25ea00102yfpc.html

转载本文请联系原作者获取授权，同时请注明本文来自林国鹏科学网博客。

链接地址：https://wap.sciencenet.cn/blog-3448646-1278095.html?mobile=1

收藏

分享