使用未知参数在R中执行fsolve(打包实用程序)

时间:2018-07-26 12:51:47

标签: r

我目前正在尝试在软件包pracma中使用功能fsolve。我想知道您是否可以帮助我解决一个小问题。

我有N个方程式的系统,带有一个未知参数(超定系统)。

让y是内生变量,x是解释变量,eta是未知参数。

等式为:Y = exp (X+eta)/(1+exp(X+eta))

代码:

Length_Vector <- 8

功能

F <- function(eta) {

y <- c(1,1,0,0,0,1,0,1)

x <- c(5,4,6,8,4,6,15,14)

x_modified <- c(x) + eta

return(c(exp(x_modified)/(1+exp(x_modified)) - y))
}

矩阵初始化

eta0 <- matrix(0, Length_Vector , 1)

计算解决方案

Solution  <- matrix(fsolve(F, eta0, maxiter=1000, tol=0.01)x, nrow = Length_Vector , ncol = 1)

当我计算此程序时,该程序运行,但是我获得了eta值不同的向量。该函数似乎可以解决由N个方程和N个未知参数(而不是1个未知参数)组成的系统。

如果您能帮助我解决这个问题,我将不胜感激。

0 个答案:

没有答案