R中rootSolve包中的空参数错误

Question

我正在使用 R 中的 rootSolve 包来求解具有 6 个未知变量的 6 个非线性方程组。这是我的模型

model <- function(x, parms) c(F1 = x[1] - parms[1] - 1 / ((-parms[7]) * (1 - x[4])),
                       F2 = x[2] - parms[2] - 1 / ((-parms[7]) * (1 - x[5])),
                       F3 = x[3] - parms[3] - 1 / ((-parms[7]) * (1 - x[6])),
                       F4 = x[4] - exp(parms[4] + parms[7] * x[1]) / (1 + exp(parms[4] + parms[7] * x[1]) + exp(parms[5] + parms[7] * x[2]) + exp(parms[6] + parms[7] * x[3])),
                       F5 = x[5] - exp(parms[5] + parms[7] * x[2]) / (1 + exp(parms[4] + parms[7] * x[1]) + exp(parms[5] + parms[7] * x[2]) + exp(parms[6] + parms[7] * x[3])),
                       F6 = x[6] - exp(parms[6] + parms[7] * x[3]) / (1 + exp(parms[4] + parms[7] * x[1]) + exp(parms[5] + parms[7] * x[2]) + exp(parms[6] + parms[7] * x[3])),
                       )

但是当我打电话给

new.equi = multiroot(model, start = initial.value, parms = parm)

在我将值传递给 initial.value 和参数的地方，我不断收到

的错误

Error in c(F1 = x[1] - parms[1] - 1/((-parms[7]) * (1 - x[4])), F2 = x[2] -  : argument 7 is empty

为什么会这样？为什么要有参数7？

我也试过像这样在模型中显式指定参数，但仍然得到同样的错误。

model <- function(x) c(F1 = x[1] - 1.265436 - 1 / (2.443700 * (1 - x[4])),
                          F2 = x[2] - 1.195844 - 1 / (2.443700 * (1 - x[5])),
                          F3 = x[3] - 1.288660 - 1 / (2.443700 * (1 - x[6])),
                          F4 = x[4] - exp(4.600528 - 2.443700 * x[1]) / (1 + exp(4.600528 - 2.443700 * x[1]) + exp(3.924360 - 2.443700 * x[2]) + exp(4.643808 - 2.443700 * x[3])),
                          F5 = x[5] - exp(3.924360 - 2.443700 * x[2]) / (1 + exp(4.600528 - 2.443700 * x[1]) + exp(3.924360 - 2.443700 * x[2]) + exp(4.643808 - 2.443700 * x[3])),
                          F6 = x[6] - exp(4.643808 - 2.443700 * x[3]) / (1 + exp(4.600528 - 2.443700 * x[1]) + exp(3.924360 - 2.443700 * x[2]) + exp(4.643808 - 2.443700 * x[3])),

)

Answer 1

你有一个尾随逗号。例如：

> c(1,2,)
Error in c(1, 2, ) : argument 3 is empty

R中rootSolve包中的空参数错误

empty argument error in rootSolve package in R

r

equation-solving

nonlinear-optimization