Mike Beddo
2010-May-03 06:20 UTC
[R] 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]]
Thomas Petzoldt
2010-May-12 06:31 UTC
[R] How to generate Mackey-Glass time series with "ddesolve" package?
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.