我的第一帖----关于nls的问题

初学R语言

各位巨巨们好，刚接触R语言，下面的代码看的不是很懂，就看懂了前面一段，后面的4段看不懂，请大家帮忙看看，解释下，不胜感激。=。-

# Here is the program that would do the optimisation of ap = f(c) to retrieve
# a*_1, a*_2 and S.
f <- function(as1=as1,as2=as2,Sf=Sf,chl=chl)
{(as1-as2)/Sf*(1.-exp(-Sf*chl))+as2*chl }


res <- nls(y ~ f(as1,as2,S,chl=chl),start=list(as1=0.1,as2=0.05,S=3.))
Sf = as.numeric(coefficients(res)[3])
astar1=as.numeric(coefficients(res)[1])
astar2=as.numeric(coefficients(res)[2])

初学R语言

nls是解非线性模型的参数，能不能解线性超定方程组呢？那解线性超定方程组怎么弄？有边界条件怎么加？

初学R语言

nls是解非线性模型的参数，能不能解线性超定方程组呢？那解线性超定方程组怎么弄？有边界条件怎么加？

南楠

res <- nls(y ~ f(as1,as2,S,chl=chl),start=list(as1=0.1,as2=0.05,S=3.))

就是一个 模拟的非线性模型吧
Sf = as.numeric(coefficients(res)[3]) 这个是找出他的系数
是可以解线性的模型的 我认为啊，因为非线性模型就是通过转变 成了线性的这个意思吧

初学R语言

南楠 发表于 2015-12-6 12:49
res

Error in numericDeriv(form[[3L]], names(ind), env) :
这是什么错误呢？

初学R语言

南楠 发表于 2015-12-6 12:49
res

Error in nlsModel(formula, mf, start, wts) :
  参数初始估计值里有奇异梯度矩阵

初学R语言

南楠 发表于 2015-12-6 12:49
res

Error in nlsModel(formula, mf, start, wts) :
  参数初始估计值里有奇异梯度矩阵

Ploter

nls对 参数的初始值很敏感  很容易出现  什么奇异梯度矩阵，我建议 先用1st opt 来确定初始值  或者先转化为线性的 然后反推 非线性的初始值 。