Question

我试图将ODE模型拟合到某些数据并求解模型中参数的值。

我知道R中有一个名为FME的软件包，旨在解决这类问题。但是，当我尝试编写此程序包的手册之类的代码时，程序无法使用以下回溯信息运行：

lsoda（y，times，func，parms，...）出错：在采取任何集成步骤之前检测到非法输入 - 请参阅书面信息

代码如下：

x <- c(0.1257,0.2586,0.5091,0.7826,1.311,1.8636,2.7898,3.8773)
y <- c(11.3573,13.0092,15.1907,17.6093,19.7197,22.4207,24.3998,26.2158)

time <- 0:7

# Initial Values of the Parameters
parms <- c(r = 4, b11 = 1, b12 = 0.2, a111 = 0.5, a112 = 0.1, a122 = 0.1)

# Definition of the Derivative Functions
# Parameters in pars; Initial Values in y
derivs <- function(time, y, pars){
     with(as.list(c(pars, y)),{
         dx <- r + b11*x + b12*y - a111*x^2 - a122*y^2 - a112*x*y
         dy <- r + b12*x + b11*y - a122*x^2 - a111*y^2 - a112*x*y
         list(c(dx,dy))
     })
}

initial <- c(x = x[1], y = y[1])

data <- data.frame(time = time, x = x, y = y)

 # Cost Computation, the Object Function to be Minimized
 model_cost <- function(pars){
     out <- ode(y = initial, time = time, func = derivs, parms = pars)
     cost <- modCost(model = out, obs = data, x = "time")
     return(cost)
 }

 # Model Fitting
 model_fit <- modFit(f = model_cost, p = parms, lower = c(-Inf,rep(0,5)))

是否有人知道如何使用FME包并解决问题？

Answer 1

您的代码语法是正确的，它可以工作到最后一行。

您可以使用

检查您的代码

model_cost(parms)

这很好用，你可以看到

model_cost(parms)$model

你的初步猜测＆＃34;远离观察到的数据（比较＆＃34; obs＆＃34;和＆＃34; mod＆＃34;）。也许这是一个失败，因此拟合程序将无法达到观察到的数据。

一时间......我还用参数＆＃34; methods = ...＆＃34;检查了不同的方法。但仍然无效。

祝福，约翰内斯

编辑：如果您使用：

model_fit <- modFit(f = model_cost, p = parms)

没有任何下限，那么你将获得一个结果（即使有警告），但是a112是负面的，你想要省略。

使用FME软件包在R中拟合ODE

1 个答案: