R-paquet deSolve: Erreur dans checkFunc

Question

Je suis assez nouveau avec R et j'essaie de résoudre un système d'équations différentielles avec le paquet de deSolve.

J'ai un total de 10000 timesteps et tous les 100 timesteps je dois changer les valeurs des variables. Le code i essayé me donne l'erreur suivante:

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) 

Je ne suis pas tout à fait sûr où l'erreur vient. Voici mon code à ce jour:

library(deSolve) 

#--------variables and parameter-------------------------- 


parameters <- c(up = 0.0001, hm = 0.1, hp = 0,889, mm = 0.1 , nt = 100) 
yini <- c(u = 1, m = 0.001, h = 0.001) 
times <- seq(0,100, by=0.1) 
out<- NULL 

#--------functions--------------------------------------- 
DiffU <- function(t, yini, parameters){ 
    with(as.list(c(yini, parameters)),{ 

    #functions 
    du <- -(up*u) + hm*h + (1/2*nt)*h 
    dm <- hp*h - mm*m - (1/nt)*m 

    # return functions 
    list(c(du, dm)) 
    }) 
} 

#--------------------------------------------- 
repeat{ 
    out<- rbind(out, ode(yini, times, DiffU, parameters)) 
    yini<-c(u = u+0.5*h, m = 0.001, h = m+0.5*h) 
    x<- x+1 
    if (x==100){ 
    break 
    } 
} 

Avez-vous des suggestions pour se débarrasser de l'erreur ou d'améliorer le code? Merci beaucoup!

Answer 1

Votre problème est que vous avez deux dérivés 'du' et 'dm', mais trois valeurs pour yini ('u', 'm' et 'h'). Cela provoque l'erreur.

Une solution consiste à (1) ajouter une troisième dérivée ou (2) supprimer 'h' de yini et l'ajouter aux paramètres.