我尝试使用r中的nls函数将以下内容拟合到ADBUG模型中,但是奇异矩阵错误不断重复,我真的不知道如何继续这样做。
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代码:
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,:奇异梯度)
发布于 2018-12-04 07:54:49
Package通过函数nls提供了一个更新版本的nlxb,在大多数情况下避免了“奇异梯度”错误。
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 拟合系数离起始值很远,且很大,说明问题的根据不充分。
https://stackoverflow.com/questions/53589971
复制相似问题