R help - May 2010 - How to generate Mackey-Glass time series with "ddesolve" package?

I could use some help generating a time series for the Mackey-Glass equation:
dx/dt = 0.2 x(t-tau)/(1 + x(t-tau)^10) - 0.1 x(t) for the case there tau = 17. I
tried the "ddesolve" package but the dde(...) function seems to hang,
not producing anything. Can someone show me the R script how to do this?


-          Mike Beddo

	[[alternative HTML version deleted]]

Hello Mike,

please provide a reproducible example, so that we can see why and how 
ddesolve "hangs". In addition you may consider using function dede
from
package deSolve which allows to choose between different solver methods 
for dealing with delay differential equations.

Thomas Petzoldt


Am 03.05.2010 08:20, schrieb Mike Beddo:
> I could use some help generating a time series for the Mackey-Glass
> equation: dx/dt = 0.2 x(t-tau)/(1 + x(t-tau)^10) - 0.1 x(t) for the
> case there tau = 17. I tried the "ddesolve" package but the
dde(...)
> function seems to hang, not producing anything. Can someone show me
> the R script how to do this?
>
>
> -          Mike Beddo
>
> [[alternative HTML version deleted]]
>
> ______________________________________________ R-help at r-project.org
> mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do
> read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.