Estimating the distribution of time to extinction of infectious diseases in mean-field approaches ================================================================================================= * Maryam Aliee * Kat S. Rock * Matt J. Keeling ## Abstract 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 recognise 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 these predictions agree very well with the results of stochastic models by analysing the simplified 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*). ## 1 Introduction For many infectious diseases the eventual aim of control measures is eradication completely removing the pathogen from host populations and the environment. This has been accomplished for just two infections – smallpox and rinderpest – but is the target for a number of other diseases including polio, Guinea Worm and Yaws [1–3]. A key question is, therefore, predicting the time to eradication, which has clear implications for the likely duration and hence cost of any control program. Markov-chain models are often used to study the spread of infectious diseases in a population [4, 5]. In such approaches the whole population is partitioned into compartments and system dynamics are described by rates of exchange between compartments. Both stochastic and deterministic versions of such models have been developed to study dynamics of a range of infectious diseases. The results of deterministic and stochastic versions of the same model generally agree well for large compartment populations, but start to deviate when there are either low numbers of infections, or for small population sizes [4, 6–9]. Deterministic models provide a mean-field approximation of infection dynamics described by ordinary differential equations (ODEs). The steady states and time evolution of these models are tractable and typically well-behaved in the whole phase space. Various analytical and numerical methods are available to solve a set of ODEs; this makes deterministic models computationally efficient and therefore allow the exploration of large regions of parameter space relatively quickly. This is particularly important for fitting such models to available data in order to estimate the underlying mechanistic parameters. Deterministic descriptions are not very practical, however, when it comes to studying rare events or very small populations. Importantly, the endpoint of an infection (when infection is eradicated) is by definition ambiguous in these models since the populations are represented by continuous variables that never reach zero. In practice, a proxy threshold is needed to determine when extinction may occur, however, it is not clear how such a threshold should be defined. Infection dropping below one may seem a natural measure, but in the literature different threshold values are applied [10, 11] generating shifts in the predicted extinction times. In such circumstances, stochastic models are more favourable. In these models, transitions between compartments are generated by stochastic events, whose average is well captured by the mean-field image, but express divergent behaviours individually. These events and accordingly system state are specified by probability functions allowing for uncommon incidents. Stochastic models are generally computationally expensive; however, they come with a good deal of flexibility and capabilities that we would like to take the advantage of to understand the extinction of diseases. In this manuscript, we will develop a basic framework to study the peri-elimination phase for an infectious disease using a simplified SIS (Susceptible-Infected-Susceptible) model that takes into account susceptible and infected compartments. We use solutions of the corresponding forward Kolmogorov equation to predict the point of infection extinction in deterministic solutions. We then extend this framework to estimate the extinction time of more complicated infection dynamics. Whilst other studies have indeed addressed similar questions around computation of mean extinction times of epidemiological or other physical systems [12–20], the present study focuses on the computation of the distribution in predicted extinction times which has not been done elsewhere. This aspect is crucially important as the policy implications of extinction occurring within shorter or longer periods could lead to different decision making. Indeed uncertainty in model prediction is not accounted for often enough in epidemiological modelling studies [21]. Representation of uncertainty allows the modeller to answer policy-driven questions such as “at what time would we expect there to be a *>*90% likelihood that elimination would be met?” or “what is the earliest we could hope to achieve zero infections?”, which the mean value does not address. For this study, we are motivated by a desire to analyse the endpoint of African sleeping sickness (*Gambiense* human African trypanosomiasis; gHAT), which is a vector-borne disease transmitted to humans by tsetse [22–25]. Humans transition through different stages of the disease and can either be diagnosed and recover (through treatment) or typically die without it. Medical intervention campaigns have had great success in bringing down the burden of this disease - from 37,385 in the most recent peak in 1998 to a record low of 953 cases in 2018 [26]. To continue the drive forwards, the World Health Organization (WHO) has a goal to achieve global elimination of transmission (EOT) of sleeping sickness by 2030 [27]. Quantitative methods provide a way to fore-cast the expected time to meet the goal under a variety of strategies, however there is a strong need for refinement of the methodology to analyse the distribution of extinction times in a more rigorous mathematical framework. Previous studies have used different threshold criteria for the endpoint for gHAT, by considering when the incidence of new infections falls below one of: 1 per 10,000, 1 per 100,000, less than 0.5 new infections across the entire health zone (*∼*100,000) per year, even to a more severe threshold of *<*1 per 1,000,000 [23, 28–30]. ## 2 Model Structure The simple SIS model is a one-dimensional characterisation of an infectious disease, in which recovery (or by treatment or otherwise) quickly leaves the individual susceptible to further infection [4, 5]. For a population of size *N*, the number of infected individuals can be expressed as: ![Formula][1] Here *β* is the transmission rate, *γ* is the recovery rate and *S* and *I* are the numbers of susceptible and infectious individuals respectively. Another important metric, *R* = *β/γ*, is the basic reproduction number. When *R* *>* 1 this model exhibits sigmoidal growth to the non-zero equilibrium (*I** = 1 *−* 1*/R*); whereas when *R* *<* 1 we observe exponential decay of infection to zero. Through changes to the transmission or recovery rates it is possible to modify the *effective reproduction number, R*eff and push it below this critical threshold. It is this latter case that is of greatest interest here, as we are most concerned with diseases that are bound for eventual extinction (assuming controls remain in place and *R*eff continues to be *<*1). The non-linear *SI* term in equation (1) precludes many analytical results; however, if we are purely interested in extinction it is acceptable to assume that the number of infected individuals is relatively small and hence that *S ≈ N*. In this limit, we obtain a simple birth-death process for infections, where birth corresponds to infecting someone and death refers to recovery [31] ![Formula][2] Here *β* plays the role of the *per capita* birth rate (*b*) and *γ* is equivalent to the *per capita* death rate (*d*). Throughout the rest of this paper we will use the parameters *b* and *d* to stress that we are approximating the dynamics as a birth-death process. Our condition that *R*eff *<* 1 now corresponds to *d > b*. ### 2.1 Dynamics without births We first consider the simplest case where *b* = *β* = 0, and so each infected individual recovers (or is treated) at a Poisson rate *d*, but no new infections occur. Starting with *n* infected individuals, the expected time for the first individual to recover is 1*/*(*dn*). After this recovery event, the population then has *n −* 1 infected individuals and the process can be iterated giving an expected time to extinction of: ![Formula][3] Figure 1A illustrates how the average extinction time *τ**n* depends on the initial number of infected people *n*, showing a logarithmic scaling for large *n*. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/12/2020.07.10.20150359/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/11/12/2020.07.10.20150359/F1) Figure 1: Extinction in a death process model. The average extinction time *τ**n* and the extinction threshold of infected population *I*ext are plotted as a function of initial number of infected population. Here, *d* = 1 per unit time, with no specific units considered. Starting with *I*(0) = *n*, the deterministic model predicts an exponential decay of the size of infected population: ![Formula][4] and hence at the expected extinction time we have: ![Formula][5] which we postulate gives us a threshold value of the deterministic model that equates to the mean extinction time of the associated stochastic model. Figure 1B shows how the threshold value changes with the initial number of people infected. For large *n*, the threshold plateaus to exp(*−E*), where *E* is Euler’s constant (0.5772) and hence the threshold is 0.5615. Therefore, we observe the useful property that the threshold is largely independent of *n*. ### 2.2 Birth-death dynamics Moving to the situation where *b* is non-zero, the dynamics are more complex. The expected time to extinction has to be formulated by considering the impact of birth and death from a given number of infected individuals (*n*). This generates an iterated set of equations for the extinction time starting with *n* infections: ![Formula][6] where the first term is the expected time until an event occurs. The other two terms take into account the possibility of a birth or death event transitioning to *n* + 1 or *n −* 1 infected individuals. There is a large relevant literature using similar approaches in epidemiology or in other contexts which lead to comparable equations for the mean extinction time [12–17, 32]. Equation 5 can be expressed as: ![Formula][7] which we solve as a matrix operation (**A***x* = *c*) for *x* = (*τ* … *τ**M*) with the boundary conditions *τ* = 0 and ![Graphic][8] for large *M*, the maximum number of initial people infected. We note that without loss of generality we can assume that *d* = 1 as this simply acts as a rescaling of time. Figure 2A shows the average extinction time calculated for different choices of *b*. The average extinction time scales with ![Graphic][9] and increases logarithmically with *n* for large values. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/12/2020.07.10.20150359/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/11/12/2020.07.10.20150359/F2) Figure 2: Extinction in a death-birth process model. The average extinction time, *τ**n*, the extinction threshold of infected population, *I*ext, and a rescaled version of this threshold are plotted as a function of initial number of infected population for three choices of birth rate *b*. Again, *d* = 1 per unit time. With the same approach as the previous section, we can calculate the threshold of number of people infected *I*ext = *n* exp((*b − d*)*τ**n*), using the predicted mean extinction time from equation 5. Figure 2B shows how the threshold depends on the initial size of the infected population, calculated for different choices of birth rate *b*. Our results confirm the useful property that the threshold is largely independent of *n* for moderate numbers of initially infected people (e.g. *n >* 500) and that a scaled threshold (by a factor of 1 *− b/d*) reaches the constant value of exp(*−E*) (Figure 2C). ### 2.3 Distribution of extinction times The above approach to formulate the mean extinction times can be extended to higher cumulative moments to estimate the variation and distribution of extinction time in the birth-death model. The second cumulative moment, ![Graphic][10], defined as the expected value of the square of the time to extinction, obeys: ![Formula][11] Rearranging the above form, and using our previous expression (Equation 5), we obtain the simplified sequential relation: ![Formula][12] This can be solved with a similar matrix operation as before (**A***x* = *c*′) with the boundary conditions *S* = 0 and ![Graphic][13] for large *M* to obtain ![Graphic][14] and hence the variance is given by ![Graphic][15] Knowing that the distribution of extinction times is highly asymmetric with a long right-hand tail, higher cumulative moments help provide a more accurate estimate of prediction intervals and the full distribution of extinction times. We find the average cube of extinction times, which is used to calculate the third cumulative moment, obeys: ![Formula][16] that can be simplified as ![Formula][17] This approach can be extended to calculate sequential relations of cumulative moments of higher orders *m* ![Formula][18] The full set of cumulative moments allows the probability distribution of extinction times to be calculated. Our numerical solutions suggest that the first four moments can give a good estimate of the distribution, approximated by a generalized F probability distribution function with four free parameters. This allows us to determine the prediction intervals for a given set of parameters (details in SI). As an alternative approach, we consider the Kolmogorov forward equations to describe the time-evolution of probabilities for the numbers infected [33]. These high-dimensional models are a numerically exact realisation of the full stochastic model, which can be computed for SIS model due to the relatively low number of states. The Kolmogorov forward equations are give by: ![Formula][19] Here *P**i* depicts the probability of having *i* individuals infected. The solution of this equation is given by **P** = exp (**A***t*) · **P****** for the vector of *P**i*s over time, where **A** is the coefficient matrix and **P****** is the initial state of the system. We can accordingly calculate the probability of extinction over time and the therefore the distribution. Figure 3 compares the estimated F distributions of extinction times from the cumulative moments with the corresponding prediction intervals with the solutions of the Kolmogorov forward equations. Again, without loss of generality we have taken *d* = 1, and varied *b < d*. The distribution shifts to the right and flattens as birth rate increases. These results demonstrate the very good agreement between the distributions, as well as the mean values and 95% prediction intervals. This analysis, therefore, provides confidence that our simple approach is a useful tool to estimate extinction times and prediction intervals for any choice of *b* and *d*. In the supplementary information (SI), we plot how variance, skewness, and kortosis change as a function initial population infected *n*. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/12/2020.07.10.20150359/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/11/12/2020.07.10.20150359/F3) Figure 3: Probability distribution of extinction time *τ**n* in the model of birth-death process starting with *n* = 1000. (A-D) Gray bars represent the solutions of Kolmogorov forward equation (light gray specifies 95% PI). The colored line shows the F distribution estimated from the first four moments of birth-death process. Corresponding mean values and 95% prediction intervals are indicated for each data set. Each panel corresponds to a birth rate *b*, when *d* = 1 per unit time for all. (E) Comparing the estimated F distribution functions plotted for all values of birth rate *b*. ## 3 Application to the complex dynamics of sleeping sickness We now investigate how well our theory holds against more realistic, and hence more complex, simulation models. Our previously developed model for gHAT [28, 34] is high dimensional, describing dynamics of vectors and two human risk groups. The model accounts for the natural history of disease progression; infected humans initially develop mild disease (“stage 1”) before progressing to more severe disease (“stage 2”) with neurological symptoms as the parasite crosses the blood-brain barrier. The infectious period ends with either death or successful diagnosis and treatment by screening programs [35]. This model is fitted to the human case incidence data recorded in different administrative regions (health zones of approximately 100,000 people) across the Democratic Republic of Congo to estimate the free parameters of the deterministic model [34]. Using a Markov Chain Monte Carlo approach, posterior parameter sets that match available case data have previously been generated (see [34] and SI for details). The deterministic equations can be numerically solved to generate the disease dynamics including the number of humans infected and new transmission each year for posterior parameter sets. We use these results and our knowledge of the underlying epidemiology to determine appropriate values of *d* and *b* that approximate the linear birth-death process (see SI for details). This allows us to generate the first four moments, and accordingly the distribution of extinction times for specific *d* and *b* rates following the framework developed in the previous sections. To verify these results, we run the stochastic discrete simulations of the gHAT model using the tau-leaping approach described elsewhere [25, 36] and select a 1-day timestep which approximates the full dynamics well due to the long timescale of gHAT infections (a Gillespie method could be used for more accurate simulations, however it is unlikely to modify the results we present here). We compute the extinction times for multiple stochastic realisations. These simulations assume that active screening (the fundamental driver of case identification and control) continues at a constant level, leading to decaying infection levels (Figures in the SI). Figure 4 compares the F probability distribution function derived from the estimated cumulants with the distribution of extinction times of one million stochastic simulations for two different parameter sets. The parameter sets represent the most likely of all posterior parameter sets found by fitting to data for two health zones named Mosango (a “moderate-risk” zone) and Kwamouth (a “high-risk” zone) [34] that reveal different infection dynamics over time, and therefore different predicted extinction times. Despite the complexity of the model for gHAT dynamics, there is still a very good agreement between the birth-death model with F distribution predictions and full stochastic gHAT infection dynamics (with differences in mean and 95% PI of less than a year) confirming the analysis based on birth-death model can deliver robust and accurate estimates of mean extinction time and the prediction intervals. Examples using other levels of screening coverage are shown in the SI. In these examples, for a fixed parameter set, we computed that the deterministic threshold, *I*ext, to compute the expected extinction time was 1.39 and 1.22 for Mosango and Kwamouth respectively. Other threshold values would need to be calculated for alternative parameter sets and alternative screening strategies for the same health zones. It should be noted that this study provides a simplified picture of eradication of a diseases taking into account instant feedback between transmission and the infected population. However in epidemiology, other measures such as elimination of transmission may be of interest. In the case of complex dynamics like gHAT, due to phase difference of humans and vectors as well as slow progression of the disease and recovering process, there is a delay between elimination of transmission and disease extinction. We speculate the lag period can be estimated based on the underlying dynamics. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/12/2020.07.10.20150359/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/11/12/2020.07.10.20150359/F4) Figure 4: Probability distribution of extinction time for the full (17-dimensional) gHAT model for two health zones, (A) Mosango and (B) Kwamouth with 30% active screening. Gray bars represent the solutions of one million stochastic simulations (light gray specifies 95% PI). The purple line shows the F distribution estimated from the first four moments of the simplified birth-death process approach (*b* = 0.38, *d* = 0.64 for Mosango and *b* = 0.25, *d* = 0.46 for Kwamouth). Corresponding mean values and 95% prediction intervals are indicated for each data set. Another challenge is to extend such analysis to allow for more general noise models, while our current stochastic models mainly account for Poisson events. In the case of gHAT dynamics, we are aware of possible over-dispersion terms that can be considered in the infection dynamics [25,34]. It is important to study how the predicted distribution of extinction time will be modified by such noise models. ## 4 Conclusion In this manuscript, we study the extinction times of diseases when *R*eff *<* 1 starting from a finite population infected, utilising results from a simplified birth-death model. We formulate the mean extinction time that is in agreement with the previous studies in the analogous limits [13–17]. Beyond that, our analysis provides a novel method to approximate the distribution of extinction times by calculating additional moments of extinction time. Our results agree well with the solutions of forward Kolmogorov equations, which offer a numerical exact prediction of the birth-death dynamics. Our findings are a vital extension to deterministic models, allowing a robust and rigorous assessment of extinction in a modelling framework that cannot naturally capture this phenomena. We considered the simplistic assumption of constant birth and death rates, however, the results can be easily extended by modifying the coefficient matrix when these rates are functions of infected population. Moreover, our analysis helps us to quantify a simple threshold that can be used to determine the mean extinction time from deterministic solutions. We showed that this threshold asymptotes to a universal plateau given by effective reproduction ratio for large initial numbers of infected individuals. This asymptotic threshold is likely to be a reasonable approximation for large populations where the number of infections is also likely to be large. Interestingly, the birth-death approach can be applied to diseases with more complex disease dynamics in order to predict their extinction times. We showed that it is valid for the example of local elimination of sleeping sickness, yet in a broad sense the dynamics are approximately SIS decaying to zero infection. Our study can be extended for SIR models, where the system state is two-dimensional by considering the recovered compartment. A non-trivial modification of the birth-death process would be helpful to consider to understand the richer dynamics of the SIR model [37]. ## Supporting information Supplementary information [[supplements/150359_file02.pdf]](pending:yes) ## Data Availability All the necessary data and details of methods are included in the supplementary information. ## Author contribution All authors contributed in developing the model and writing the manuscript. ## Data Accessibility This work does not include any original data. To describe the gHAT dynamics, we used the model and the fitted parameters from a published manuscript which is cited in the main text and discussed in details in the supplementary material. ## Acknowledgements The authors thank Dr Erick Mwamba Miaka and PNLTHA for original data collection, and WHO for data access (in the framework of the WHO HAT Atlas [27]). This work was supported by the Bill and Melinda Gates Foundation ([www.gatesfoundation.org](http://www.gatesfoundation.org)) through the NTD Modelling Consortium [OPP1184344] (M.A., K.S.R. and M.J.K.) and the Human African Trypanosomiasis Modelling and Economic Predictions for Policy (HAT MEPP) project [OPP1177824] (K.S.R. and M.K.). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. ## Footnotes * The title is changed to be more specific. There were a couple of paragraphs added to highlight the importance of our results. * Received July 10, 2020. * Revision received November 12, 2020. * Accepted November 12, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1].Alan R. Hinman. Global progress in infectious disease control. Vaccine, 16(11):1116–1121, 1998. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0264-410X(98)80107-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9682367&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000074539100008&link_type=ISI) 2. [2].Donald R. Hopkins. Disease eradication. New England Journal of Medicine, 368(1):54–63, 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMra1200391&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23281976&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312930600005&link_type=ISI) 3. [3]. Garrett T. Nieddu, Lora Billings, James H. Kaufman, Eric Forgoston, and Simone Bianco. Extinction pathways and outbreak vulnerability in a stochastic ebola model. Journal of The Royal Society Interface, 14(127):20160847, 2017. 4. [4]. Matt J Keeling and Pejman Rohani. Modeling infectious diseases in humans and animals. Princeton University Press, 2011. 5. [5]. R.M. Anderson and R.M. May. Infectious Diseases of Humans. Oxford Science Publications, Oxford, 1991. 6. [6].Ingemar Nåsell. Stochastic models of some endemic infections. Mathematical Biosciences, 179(1):1–19, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0025-5564(02)00098-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12047919&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176531300001&link_type=ISI) 7. [7]. M.J Keeling and J.V Ross. On methods for studying stochastic disease dynamics. Journal of The Royal Society Interface, 5(19):171–181, 2008. 8. [8]. Matt J. Keeling, Howard B. Wilson, and Steve W. Pacala. Reinterpreting space, time lags, and functional responses in ecological models. Science, 290(5497):1758–1761, 2000. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyOTAvNTQ5Ny8xNzU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTEvMTIvMjAyMC4wNy4xMC4yMDE1MDM1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 9. [9].Matt J. Keeling. Simple stochastic models and their power-law type behaviour. Theoretical Population Biology, 58(1):21–31, 2000. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/tpbi.2000.1475&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10944473&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089076600002&link_type=ISI) 10. [10]. T Déirdre Hollingsworth, Emily R Adams, Roy M Anderson, Katherine Atkins, Sarah Bartsch, María-gloria Basáñez, Matthew Behrend, David J Blok, Lloyd A C Chapman, Luc Coffeng, Orin Courtenay, Ron E Crump, Sake J De Vlas, Andy Dobson, Louise Dyson, Hajnal Farkas, Alison P Galvani, Manoj Gambhir, David Gurarie, Michael A Irvine, Sarah Jervis, Matt J Keeling, Louise Kelly-hope, Charles King, Bruce Y Lee, Epke A Le Rutte, Thomas M Lietman, Martial Ndeffo-mbah, Graham F Medley, Edwin Michael, Abhishek Pandey, Jennifer K Peterson, Amy Pinsent, Travis C Porco, Jan Hendrik Richardus, Lisa Reimer, Kat S Rock, Brajendra K Singh, Wilma Stolk, Subramanian Swaminathan, Steve J Torr, Jeffrey Townsend, James Truscott, Martin Walker, Alexandra Zoueva, and N T D Modelling Consortium. Quantitative analyses and modelling to support achievement of the 2020 goals for nine neglected tropical diseases. Parasites & Vectors, 8(1):630, 2015. 11. [11]. James E. Truscott, Marleen Werkman, James E. Wright, Sam H. Farrell, Rajiv Sarkar, Kristjana Ásbjörnsdóttir, and Roy M. Anderson. Identifying optimal threshold statistics for elimination of hookworm using a stochastic simulation model. Parasites Vectors, 10, 2017. 12. [12].Yasushi Masuda Ushio Sumita. On first passage time structure of random walks. Stochastic Processes and their Applications, 20(1):133–147, 1985. 13. [13].Ingemar Nåsell. Extiction and quasi-stationarity in the stochastic logistic SIS model, vol 2022. Lecture Notes in Mathematics. Springer, Berlin, 2011. 14. [14]. Tom Britton and Peter Neal. The time to extinction for a stochastic sis-householdepidemic model. Journal of Mathematical Biology, 61, 2010. 15. [15]. Damian Clancy and Elliott Tjia. Approximating time to extinction for endemic infection models. Methodology and Computing in Applied Probability, 20:1043–1067, 2018. 16. [16]. Håkan Andersson and Tom Britton. Stochastic epidemics in dynamic populations: quasi-stationarity and extinction. Journal of Mathematical Biology, 41:559–580, 2000. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s002850000060&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11196585&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000166066200005&link_type=ISI) 17. [17]. Charles R. Doering, Khachik V. Sargsyan, and Leonard M. Sander. Extinction times for birth-death processes: Exact results, continuum asymptotics, and the failure of the fokker–planck approximation. Multiscale Modeling & Simulation, 3(2):283–299, 2005. 18. [18]. Mark I Dykman, Ira B Schwartz, and Alexandra S Landsman. Disease extinction in the presence of random vaccination. Phys. Rev. Lett., 101(7):078101, August 2008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevLett.101.078101&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18764580&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) 19. [19]. Eric Forgoston and Richard O. Moore. A primer on noise-induced transitions in applied dynamical systems. SIAM Review, 60(4):969–1009, 2018. 20. [20]. Michael Assaf and Baruch Meerson. Extinction of metastable stochastic populations. Phys. Rev. E, 81:021116, Feb 2010. 21. [21].Probabilistic uncertainty analysis of epidemiological modeling to guide public health intervention policy. Epidemics, 6:37–45, 2014. 22. [22]. Kat S Rock, Steve J Torr, Crispin Lumbala, and Matt J Keeling. Predicting the impact of intervention strategies for sleeping sickness in two high-endemicity health zones of the Democratic Republic of Congo. PLoS neglected tropical diseases, 11(1):e0005162, 2017. 23. [23]. K.S. Rock, A. Pandey, M.L. Ndeffo-Mbah, K.E. Atkins, C. Lumbala, A. Galvani, and M.J. Keeling. Data-driven models to predict the elimination of sleeping sickness in former Equateur province of DRC. Epidemics, 18:101–112, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2017.01.006&link_type=DOI) 24. [24]. K S Rock, M L Ndeffo-Mbah, Soledad Castaño, Cody Palmer, A Pandey, K E Atkins, Joseph M Ndung’u, T D Hollingsworth, A P Galvani, Caitlin Bever, N Chitnis, and M J Keeling. Assessing Strategies Against Gambiense Sleeping Sickness Through Mathematical Modeling. Clinical Infectious Diseases, 66(suppl. 4):S286–S292, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciy018&link_type=DOI) 25. [25]. M Soledad Castaño, Maryam Aliee, Erick Mwamba Miaka, Matt J Keeling, Nakul Chitnis, and Kat S Rock. Screening strategies for a sustainable endpoint for gambiense sleeping sickness. The Journal of Infectious Diseases, page jiz588, 2019. 26. [26].World Health Organization. Global health observatory data. 2019. 27. [27]. José R. Franco, Giuliano Cecchi, Gerardo Priotto, Massimo Paone, Abdoulaye Diarra, Lise Grout, Pere P. Simarro, Weining Zhao, and Daniel Argaw. Monitoring the elimination of human african trypanosomiasis: Update to 2016. PLoS Neglected Tropical Diseases, 12(12):e0006890, 2018. 28. [28]. Kat S Rock, Steve J Torr, Crispin Lumbala, and Matt J Keeling. Quantitative evaluation of the strategy to eliminate human African trypanosomiasis in the DRC. Parasites & Vectors, 8(1):532, 2015. 29. [29]. M Soledad Castaño, Martial L Ndeffo-Mbah, Kat S Rock, Cody Palmer, Edward Knock, Erick Mwamba Miaka, Joseph M Ndung’u, Steve Torr, Paul Verlé, Simon E F Spencer, and Others. Assessing the impact of data aggregation in model predictions of HAT transmission and control activities. medRxiv, page 19005991, 2019. 30. [30]. C Huang, R E Crump, P Brown, S E F Spencer, E Mwamba Miaka, S Chancy, M J Keeling, and K S Rock. Shrinking the ghat map: identifying target regions for enhanced control in drc. medRxiv, 2020. 31. [31]. Tom Britton. Stochastic epidemic models: A survey. Mathematical Biosciences, 225(1):24–35, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mbs.2010.01.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20102724&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) 32. [32]. Ushio Sumita. On conditional passage time structure of birth-death processes. Journal of Applied Probability, 21(1):10–21, 1984. 33. [33]. A. N. Kolmogorov. Foundations of the theory of probability. New York: Chelsea Publishing, 1993. 34. [34]. Ronald E Crump, Ching-I Huang, Ed Knock, Simon E.F. Spencer, Paul E. Brown, Erick Mwamba Miaka, Shampa Chancy, Matt J Keeling, and Kat S. Rock. Quantifying epidemiological drivers of gambiense human african trypanosomiasis across the democratic republic of congo. 2020. 35. [35].World Health Organization. Control and surveillance of human african trypanoso-miasis: report of a who expert committee. Technical report, 2013. 36. [36]. Christopher N Davis, Kat S Rock, Erick Mwamba Miaka, and Matt J Keeling. Village-scale persistence and elimination of gambiense human African trypanosomiasis. medRxiv, page 19006502, 2019. 37. [37]. Matt J Keeling and Bryan T Grenfell. Understanding the persistence of measles: reconciling theory, simulation and observation. Proceedings. Biological sciences, 269(1489):335–343, 02 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2001.1898&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11886620&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F12%2F2020.07.10.20150359.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174327000002&link_type=ISI) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/graphic-8.gif [8]: /embed/inline-graphic-1.gif [9]: /embed/inline-graphic-2.gif [10]: /embed/inline-graphic-3.gif [11]: /embed/graphic-10.gif [12]: /embed/graphic-11.gif [13]: /embed/inline-graphic-4.gif [14]: /embed/inline-graphic-5.gif [15]: /embed/inline-graphic-6.gif [16]: /embed/graphic-12.gif [17]: /embed/graphic-13.gif [18]: /embed/graphic-14.gif [19]: /embed/graphic-15.gif