R nls:找不到函数“；a「；

我试图用nls（）函数来拟合这些信息。 公式是：

但他可能会生气：）

dk1 Q H 1 43.754 21 2 74.434 52 3 94.601 72 4 115.327 81 5 135.751 109 6 72.990 54 7 35.756 12 8 36.685 2 9 35.799 2 >k1_H0 [1] 0 >n_n您可以尝试：

nls(Q ~ a*((H-k1_HO)^b)*exp(c*(log10(H-k1_HO)^2)), data = dk1, start = list( a = 0.1, b = 0.1, c=0.1), trace = T)

它应该与适当的起始值一起给出结果（否则错误会提到“奇异梯度”）

---

R需要一个*符号才能知道它是一个多应用程序，否则a（）或c（）将被视为函数。对于指数，你应该检查一下？exp.

你在

a（（H-k1_h0）中缺少了
*
。

也应该是

exp（c*（log10（H-k1_h0））**2）

因为

exp（4）

和

exp（1）^4

和

exp（1）**4

是相同的

nls(Q ~ a*((H-k1_HO)^b)*exp(c*(log10(H-k1_HO)^2)), data = dk1, start = list( a = 0.1, b = 0.1, c=0.1), trace = T)