Estimating the distribution of time to extinction of infectious diseases in mean-field approaches

A key challenge for many infectious diseases is to predict the time to extinction under specific interventions. In general, this question requires the use of stochastic models which recognize the inherent individual-based, chance-driven nature of the dynamics; yet stochastic models are inherently computationally expensive, especially when parameter uncertainty also needs to be incorporated. Deterministic models are often used for prediction as they are more tractable; however, their inability to precisely reach zero infections makes forecasting extinction times problematic. Here, we study the extinction problem in deterministic models with the help of an effective ‘birth–death’ description of infection and recovery processes. We present a practical method to estimate the distribution, and therefore robust means and prediction intervals, of extinction times by calculating their different moments within the birth–death framework. We show that these predictions agree very well with the results of stochastic models by analysing the simplified susceptible–infected–susceptible (SIS) dynamics as well as studying an example of more complex and realistic dynamics accounting for the infection and control of African sleeping sickness (Trypanosoma brucei gambiense).