我无法让nls函数与r一起使用。我究竟做错了什么?这是我的代码
data = read.table("http://math.univ-lyon1.fr/~gciuperca/lactt.txt", header = FALSE);
plot(data$V1, data$V2)
model <- nls(V2 ~ (a*(V1 + 1)^b)*exp(-c*V1), start = c(a = 0, b = 2, c = 0), data=data)
这是我得到的错误(以法语表示抱歉)
Erreur dans nlsModel(formula, mf, start, wts) :
matrice de gradient singulière pour les estimations initiales des paramètres
您是否尝试过其他起始值?
nls(V2 ~ (a*(V1 + 1)^b)*exp(-c*V1), start = c(a = 1, b = 2, c = .1), data=data)
# Nonlinear regression model
# model: V2 ~ (a * (V1 + 1)^b) * exp(-c * V1)
# data: data
# a b c
# 15.001 0.180 0.003
# residual sum-of-squares: 0.001354
#
# Number of iterations to convergence: 9
# Achieved convergence tolerance: 1.39e-08
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句