1

我正在尝试将我的数据拟合到用户定义的 Gompertz 方程。但是，我收到如下错误：

nls (Eq, data = GPdata, start = list(K = 1), algorithm = "default", : 奇异梯度中的错误

我知道这可能意味着我的起始值不足以收敛，但是除了猜测之外还有什么办法吗？我还尝试使用 nls2 函数（蛮力），扩展可能的 K 值的巨大向量以尝试但收到相同的错误。

这是我的代码，包括示例数据。

t = c(0, 2.4, 4, 5.9, 8.3, 10.2, 12, 14.5, 16.1, 17.8, 19, 21.1, 23, 24.7, 26.4, 28.3,
    30.4, 32.3, 34.2, 36)
Cells = c(14619994708, 18945074477, 19999450160, 23461507263, 22803949639, 19400861775,
        19806214226, 17390826226, 17390826226, 18186961935, 14066488593,
        13285017004, 11058253613, 8085853098, 9168523935, 8435042020, 
        8560592341, 7734294268, 7720887693, 7503367578)

GPdata = data.frame(t, Cells)

Ninf = (1 * 10 ^ 6) * (1 / 0.001)
Eq <- Cells ~ exp(log(Ninf) * (1 - exp(-K * t)))
Gompstart = nls(Eq, data = GPdata, start = list(K = 1), 
        algorithm = "default", trace = TRUE)

任何帮助或建议将不胜感激。

谢谢！

4

1 回答 1

3

1) 请注意，对于所有t和K：

LHS >= min(LSH) = min(Cells) = 7503367578 > 1e9 = Ninf = max(RHS) >= RHS

因此所有的 LHS > RHS t，K因此该模型不起作用。

2）指数不会模拟上升然后下降的现象。

于 2013-04-30T19:43:02.017 回答