Log In Sign Up

Numerical methods and hypoexponential approximations for gamma distributed delay differential equations

by   Tyler Cassidy, et al.

Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic Gamma distributed DDEs are not currently available. Accordingly, modellers often resort to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations. In this work, we develop a functionally continuous Runge-Kutta method to numerically integrate the gamma distributed DDE and perform numerical tests to confirm the accuracy of the numerical method. As the functionally continuous Runge-Kutta method is not available in most scientific software packages, we then derive hypoexponential approximations of the gamma distributed DDE. Using our numerical method, we show that while using the common Erlang approximation can produce solutions that are qualitatively different from the underlying gamma distributed DDE, our hypoexponential approximations do not have this limitation. Finally, we implement our hypoexponential approximations to perform statistical inference on synthetic epidemiological data.


page 1

page 2

page 3

page 4


Magnus integrators for linear and quasilinear delay differential equations

A procedure to numerically integrate non-autonomous linear delay differe...

Efficiently and easily integrating differential equations with JiTCODE, JiTCDDE, and JiTCSDE

We present a family of Python modules for the numerical integration of o...

A New family of methods for solving delay differential equations

In the present paper, we introduce a new family of θ-methods for solving...

Periodic Pólya Urns, the Density Method, and Asymptotics of Young Tableaux

Pólya urns are urns where at each unit of time a ball is drawn and repla...

Bifurcation analysis of Bogdanov-Takens bifurcations in delay differential equations

In this paper, we will perform the parameter-dependent center manifold r...