Hi,
to diagnose your problem, I tried the model with your original
parameters, using a Runge-Kutta fixed step integrator and smaller time
steps. You may try to make them even smaller.
rk4 has no warranted accuracy, so it is less reliable than "lsoda" etc.
However, it can be useful for
analysis. I tried different set of parameters, but it
>> still gives me the same error message. Maybe this has to do with my
>> equations ? But I really doubt it.
>>
>> Thank you for your time
>> Alex
>>
>> ----------
>> From: abdel.ha
rs, but it
> still gives me the same error message. Maybe this has to do with my
> equations ? But I really doubt it.
>
> Thank you for your time
> Alex
>
> --
> From: abdel.hallo...@gmail.com
> Date: Sat, 20 Feb 2016 11:01:00 -0600
&
analysis. I tried different set of parameters, but it still
gives me the same error message. Maybe this has to do with my equations ? But I
really doubt it.
Thank you for your time
Alex
From: abdel.hallo...@gmail.com
Date: Sat, 20 Feb 2016 11:01:00 -0600
Subject: Re: [R] Problems
I think your parameters are off. If you look at the simul data frame, it
gives you a bunch of NaNs after the first initialization. If you put lower
the timesteps s.t.
> times<-seq(0,200, by=0.01)
it begins to run but soon your values diverge, i1 & i2 going negative while
i12 goes way high. Not
Hello R-users,
I'm trying to build a SIR-like model under R,
using the "deSolve" package. I'm trying to do simulations of its dynamic
over time, with three differential equations. I'm also looking to
calculate the equilibrium state.
So far, my code looks like this
library(deSolve)
#This is
6 matches
Mail list logo
