Approximation of epidemic models by diffusion processes and their statistical inference

TL;DR

This paper develops a new statistical inference framework for epidemic models by approximating Markov jump processes with diffusion processes, enabling efficient parameter estimation from incomplete data.

Contribution

It generalizes previous diffusion approximations to non-autonomous systems and extends inference results to time-dependent diffusions, providing a practical analytical method.

Findings

Estimators are consistent and asymptotically Gaussian.

Performance is robust across various parameters like R0, N, and n.

Method performs well with realistic data sizes and population scales.

Abstract

Multidimensional continuous-time Markov jump processes $(Z(t))$ on $\mathbb{Z}^p$ form a usual set-up for modeling $SIR$-like epidemics. However, when facing incomplete epidemic data, inference based on $(Z(t))$ is not easy to be achieved. Here, we start building a new framework for the estimation of key parameters of epidemic models based on statistics of diffusion processes approximating $(Z(t))$. First, \previous results on the approximation of density-dependent $SIR$-like models by diffusion processes with small diffusion coefficient $\frac{1}{\sqrt{N}}$, where $N$ is the population size, are generalized to non-autonomous systems. Second, our previous inference results on discretely observed diffusion processes with small diffusion coefficient are extended to time-dependent diffusions. Consistent and asymptotically Gaussian estimates are obtained for a fixed number $n$ of observations, which corresponds to the epidemic context, and for $N\rightarrow \infty$. A correction term, which yields better estimates non asymptotically, is also included. Finally, performances and robustness of our estimators with respect to various parameters such as $R_0$ (the basic reproduction number), $N$, $n$ are investigated on simulations. Two models, $SIR$ and $SIRS$, corresponding to single and recurrent outbreaks, respectively, are used to simulate data. The findings indicate that our estimators have good asymptotic properties and behave noticeably well for realistic numbers of observations and population sizes. This study lays the foundations of a generic inference method currently under extension to incompletely observed epidemic data. Indeed, contrary to the majority of current inference techniques for partially observed processes, which necessitates computer intensive simulations, our method being mostly an analytical approach requires only the classical optimization steps.