🔍 检索绿色通道
     0.ESI查询      1.JCR-IF查询      2.CiteScore-IF查询     3.SCI检索     4.中科院分区查询     5.SSCI检索     6.AHCI检索     7.ESCI检索(新增)     8.EI检索     9.EI检索备用     10.ISTP/CPCI检索备用(新增)      11.梅斯期刊智能查询     12.CSCD检索     13.CSSCI检索     14.北大核心检索      15.国自然官方查询     16.国自然梅斯查询     17.国自然科学网查询     18.论文被检索情况查询方法

R中y=a*exp(b*x)+c的selfStart函数


R中y=aexp(bx)+c的selfStart函数

Weiping Mei < meiweipingg@gmail.com >

# selfStart method for
# y = a *exp(b*x)+ c

expModel<-function(predictor,a,b,c){a*exp(b*predictor)+c}  #定义对应模型的函数。

  expModelInit <- function(mCall,LHS, data)
  {
xy <- sortedXyData(mCall[["predictor"]],LHS, data)
y=xy[,"y"]
x=xy[,"x"]
adjy=y-min(y)+1
xadjy=data.frame(x,adjy)

lmFit <- lm(log(adjy) ~ x)
coefs <- coef(lmFit)
get.b <- coefs[2]   #slope

nlsFit<-nls(adjy~cbind(1+exp(b*x),exp(b*x)),
            start = list(b=get.b),data = xadjy,algorithm = "plinear",
            nls.control(maxiter = 5000000,minFactor = 10^(-10)))

coef<-coef(nlsFit)
b<-coef[1]
c<-coef[2]+min(y)-1
a<-coef[3]+coef[2]

value <- c(a,b,c)
names(value) <- mCall[c("a","b","c")]
value
}

SSexp3P<-selfStart(expModel,expModelInit,c("a","b","c"))

# getInitial(y~SSexp3P(x,a,b,c),data = z)

您能在R中通过github安装 “fitting”函数包,其中的”SSexp3P”的函数即可得到 “y=a*exp(b*x)+c”回归的各个参数的值(迭代结果):

x<-1:5
y<-c(2,14,16,18,19)
z<-data.frame(x,y)

fitexp3Pss<-nls(y~SSexp3P(x,a,b,c),data=z)
fitsum <- summary(fitexp3Pss)   # 输出各参数最终结果
fitsum


<已有 次阅读>


由于本文作者水平有限,文中如有错误之处,欢迎大家批评指正!

① 本文仅代表作者个人观点,不代表任何其它立场,欢迎交流合作!

② 转载与分享请注明:本文源于为学为研网 http://meiweiping.cn

(>看完记得五星好评哦亲<)

分享到:
    友情链接
    e教学: 蓝墨云班课     传课网     慕课网     作业帮
    BBS: 小木虫     丁香园     科学网bbs     零点花园     数据狗论坛     ResearchGate
    预印本: 中国科技论文在线(教育部)     ChinaXiv(中科院)     bioRxiv(冷泉港)     arXiv(康奈尔大学)
    科研资讯: 科学网   科学之家     51science     Publons     J-STAGE     日本の研究
    信息查询: 中国研究生招生信息网     学信网(学历查询)     中国人事考试网     更多>>
    网站建设: Hexo     极简图床     sm.ms     html-unicode     UTF-8转码     颜色代码
    关于本站: 关于我们      Logo含义      广告合作      免责声明      联系我们      RSS订阅