Pairwise Approximations of Non-markovian Network Epidemics

Joint work with Zsolt Vizi (Bolyai Institute, University of Szeged, Hungary), Istvan Kiss (Department
of Mathematics, University of Sussex, United Kingdom)

Pairwise models have been proven to be a flexible framework for analytical approximations
of stochastic epidemic processes on networks that are in many situations much more accurate
than mean field compartmental models. The non-Markovian aspects of disease transmission
are undoubtedly important, but very challenging to incorporate them into both numerical
stochastic simulations and analytical investigations. Here we present a generalization of
pairwise models to non-Markovian epidemics on networks. For the case of infectious periods
of fixed length, the resulting pairwise model is a system of delay differential equations, which
shows excellent agreement with results based on the explicit stochastic simulations. For more
general distribution classes (uniform, gamma, lognormal etc.) the resulting models are PDEs
that can be transformed into systems of integro-differential equations. We derive pairwise
reproduction numbers and relations for the final epidemic size, and initiate a systematic
study of the impact of the shape of the particular distributions of recovery times on how
the time evolution of the disease dynamics play out.