R nls：找不到函數“ a”

Question

我正在嘗試使用nls（）函數適應信息。 公式為：

但是R會生氣:)

> 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_H0)^b)*(exp(1)^(c(log10(H-k1_H0)^2))), data = dk1, start = list( a = 0.1, b = 0.1, c=0.1), trace = T)
Error in eval(expr, envir, enclos) : could not find function "a"

Answer 1

您可以嘗試：

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。