解决R中的非线性方程组

时间:2018-08-02 13:21:55

标签: r

我想求解一个非线性方程组,但是我陷入了R代码中。

install.packages("nleqslv")
library(nleqslv)
se_function<-function(A,b,c){
  n<-nrow(Se_data)
  k1<-Se_data$HTA.
  k2<-Se_data$RMA.
  k3<-Se_data$CMA.
  F<-Se_data$rc_mean
  Area<-Se_data$MTA
  F[n]<-A+b[1]*exp(k1[n]*Area[n]/c[1])+b[2]*exp(k2[n]*Area[n]/c[2])+b[3]*exp(k3[n]*Area[n]/c[3])
  F
}
startX<-rnorm(50)
se_result<-dfsane(startX,se_function,control = list(maxit=500,trace=FALSE))

我要求解的方程:

在这个等式中,我已经知道x,y和k [1,2,3],并且我想知道A,b [1,2,3]和c [1,2, 3]。 x,y和k有25个小组,但我收到了以下消息

Error in fn(par, ...) : argument "b" is missing, with no default
Error in dfsane(startX, se_function, control = list(maxit = 500, trace = FALSE)) : 
  Failure in initial functional evaluation.

请帮助我,您能解释一下如何使用function()吗? 非常感谢。

0 个答案:

没有答案