leerer Argumentfehler im rootSolve-Paket in R

Question

Ich verwende rootSolve-Paket in R, um ein System von 6 nichtlinearen Gleichungen mit 6 unbekannten Variablen zu lösen. Hier ist mein Modell

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])), 
         ) 

Aber wenn ich rufe

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

wo ich Wert übergeben initial.value und Parm, halte ich immer den Fehler

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

Warum ist das passiert ? Warum sollte es Argument 7 geben?

Ich habe auch versucht, Parameter im Modell explizit anzugeben, wie dies, aber immer noch den gleichen Fehler.

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

Sie haben ein Komma. Zum Beispiel:

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