首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >r模型nls奇异梯度

r模型nls奇异梯度
EN

Stack Overflow用户
提问于 2018-12-03 08:29:14
回答 1查看 74关注 0票数 1

我尝试使用r中的nls函数将以下内容拟合到ADBUG模型中,但是奇异矩阵错误不断重复,我真的不知道如何继续这样做。

代码语言:javascript
复制
     nprice       nlv2
[1,] 0.6666667 1.91666667
[2,] 0.7500000 1.91666667
[3,] 0.8333333 1.91666667
[4,] 0.9166667 1.44444444
[5,] 1.0000000 1.00000000
[6,] 1.0833333 0.58333333
[7,] 1.1666667 0.22222222
[8,] 1.2500000 0.08333333
[9,] 1.3333333 0.02777778

代码:

代码语言:javascript
复制
fit <- nls(f=nprice~a+b*nlv2^c/(nlv2^c+d),start=list(a=0.083,b=1.89,c=-10.95,d=0.94))

nls中的错误(f= nprice ~a+b* nlv2^c/(nlv2^c + d),start = list(a = 0.083,:奇异梯度)

EN

回答 1

Stack Overflow用户

发布于 2018-12-04 07:54:49

Package通过函数nls提供了一个更新版本的nlxb,在大多数情况下避免了“奇异梯度”错误。

代码语言:javascript
复制
library(nlsr)
fit <- nlxb(f = nprice~a+b*nlv2^c/(nlv2^c+d),
       data = df,
       start = list(a=0.083,b=1.89,c=-10.95,d=0.94))
## vn:[1] "nprice" "a"      "b"      "nlv2"   "c"      "d"     
## no weights

fit$coefficients
##           a           b           c           d 
## -2.1207e+04  2.1208e+04 -7.4083e-01  1.6236e-05 

拟合系数离起始值很远,且很大,说明问题的根据不充分。

票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/53589971

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档