Computational uncertainty quantification for random time‐discrete epidemiological models using adaptive gPC

Abstract

Population dynamics models consisting of nonlinear difference equations allow us to get a better understanding of theprocesses involved in epidemiology. Usually, these mathematical models are studied under a deterministic approach.However, in order to take into account the uncertainties associated with the measurements of the model input param-eters, a more realistic approach would be to consider these inputs as random variables. In this paper, we study therandom time-discrete epidemiological models SIS, SIR, SIRS, and SEIR using a powerful unified approach basedupon the so-called adaptive generalized polynomial chaos (gPC) technique. The solution to these random differenceequations is a stochastic process in discrete time, which represents the number of susceptible, infected, recovered,etc individuals at each time step. We show, via numerical experiments, how adaptive gPC permits quantifying theuncertainty for the solution stochastic process of the aforementioned random time-discrete epidemiological modeland obtaining accurate results at a cheap computational expense. We also highlight how adaptive gPC can be appliedin practice, by means of an example using real data.