Question

我试图将deSolve包用于一组ODE，其中方程为辅助变量。我不断得到误差，其中导数的数量与初始条件向量的长度不同。我应该改变什么？

# rm(list=ls())  
library(deSolve)

exponential=function(t,state,parameters){ with(as.list( c(state,parameters)), {

   #Aux. Var.
   fX2 = pmax(0,1-(1-(d2/r12)*(X2/K2)))
   fX1 = X1/(X1+k1); 

   # equations (ODE)
   dX1 = C-((d1)*(X1))-(r12)*(X2)*fX2*fX1 # differential equaion
   dX2 = r12*(X2)*fX2*fX1-((d2)*(X2))

   return(list(c(dX1, dX2)))
   })
 }

# -- RUN INFORMATION

# Set Initial Values and Simulation Time
state = c(X1=2,X2=0.01,K2= 10)
times=0:100

# Assign Parameter Values
parameters = c(d1=0.001, d2=0.008, r12=0.3,C=0.5,k1= 0.001)

for (i in 1:length(times)){
  out= ode(y=state,times=times,func=exponential,parms=parameters)
  }

Error in checkFunc(Func2, times, y, rho) : 
  The number of derivatives returned by func() (2) must equal the length of
 the initial conditions vector (3)**

Answer 1

错误来自您定义的函数中的return：您的输入参数y的长度为3，但您只返回2个值，这是错误。您可以使用

解决问题

return(list(c(X1, X2, K2)))

另一种可能性是将K2带到参数，然后您的旧return是正确的。您必须确定K2是变量还是参数。

而BTW：为什么一个for循环与时间？在我看来，这是不必要的，因为ODE是在你提交给ode函数的时间间隔内解决的。

out= ode(y=state,times=times,func=exponential,parms=parameters)

使用deSolve在R-个导数误差中求解ODE

1 个答案: