初始参数估计时的 nls 奇异梯度矩阵

Question

我试图将第一个零级动力学曲线拟合到一些数据，当我绘制数据并向模型添加一些初始参数时，我得到了非常好的视觉拟合，但是当我运行和nls我得到错误：

nlsModel(formula, mf, start, wts) 中的错误：初始参数估计时的奇异梯度矩阵

我已经尝试过参数但没有成功，所以我根据模型模拟了一些数据，但仍然出现相同的错误。 任何可能出问题的建议：

t <- 1:50
CrStart <- 850
CsStart <- 1100
KrStart <- 0.27
KsStart <- 0.0065
y3 <- runif(1,825,875)*(1 - exp(-runif(1,0.25,0.3)*t)) + runif(1,1050,1150)*runif(1,0.00625,0.00675)*t+runif(50,0,100)

plot(t,y3)

m <- nls (y3~Cr*(1 - exp(-Kr*t)) + Cs*Ks*t, start=list(Cr=CrStart,Cs=CsStart,Kr=KrStart,Ks=KsStart))

Answer 1

您的问题是参数Cs和Ks被混淆了； 您的模型中有太多参数。 有无数种Cs和Ks组合适合该模型。

解决方案是将Cs*Ks更改为K并改为求解K 那么K = Cs*Ks 。