Modeling the Post-Containment Elimination of Transmission of COVID-19 ===================================================================== * Flávio Codeço Coelho * Luiz Max Carvalho * Raquel M Lana * Oswaldo G Cruz * Leonardo S Bastos * Claudia T Codeço * Marcelo F C Gomes * Daniel Villela ## Abstract Roughly six months into the COVID-19 pandemic, many countries have managed to contain the spread of the virus by means of strict containment measures including quarantine, tracing and isolation of patients as well strong restrictions on population mobility. Here we propose an extended SEIR model to explore the dynamics of containment and then explore scenarios for the local extinction of the disease. We present both the deterministic and stochastic version fo the model and derive the ℛ and the probability of local extinction after relaxation (elimination of transmission) of containment, ℙ. We show that local extinctions are possible without further interventions, with reasonable probability, as long as the number of active cases is driven to single digits and strict control of case importation is maintained. The maintenance of defensive behaviors, such as using masks and avoiding agglomerations are also important factors. We also explore the importance of population immunity even when above the herd immunity threshold. Keywords * COVID-19 * Epidemiology * Mathematical modelling * pandemic * SARS-COV-2 * probability of extinction ## 1. Introduction The COVID-19 pandemic is among the top three biggest in the last on hundred years, reaching levels only previously seen in Influenza pandemics[1]. At the time of writing, more than 8 million confirmed cases have been reported globally, with more than 500 thousand deaths [2]. Most affected countries still observe transmission, even if number of new cases show signs of reduction. Even with few cases, if cities (or countries) decide to lift quarantine measures, transmission may increase again due to the presence of sizeable portion of susceptible individuals suddenly at greater risk of infection. The initial containment response varied considerably across countries, with some countries displaying more success than others in avoiding infections and subsequent deaths[3, 4, 5]. The economic impact of the containment efforts in the form of quarantines, lock-downs, suspension of international travelling and other drastic measures, has been remarkable [6]. This economic strain has forced many countries towards an early suspension of many of the most severe containment measures such as quarantines and mobility restrictions [7]. Naturally, the re-normalization of social interaction brings with it may concerns about the potential for a second wave of transmissions, which could potentially grow out of control [8]. Thus, an important question that emerges after a period of isolation is: *What’s the probability that local infections will be eliminated once the isolation is lifted?*. Evaluating scenarios for incidence evolution after these initial containment efforts requires models which accommodate both biological and population-level dynamics. In particular, models that can represent properly the immunological aspects of COVID-19 progression as well as the impact of containment mechanisms, such as quarantine and social distancing. Many models have been proposed recently to deal with the temporal evolution of the epidemic[9, 10, 11], but one key aspect that must be considered is the contribution of stochastic fluctuations to the interruption of local transmission after the number of active cases is brought close to zero. Kucharsky et al. (2020, [12]), used an stochastic transmission model to estimate the daily reproduction number, *ℛ**t*, which is often used to predict disease extinction (ℛ*t* *<* 1), but empirical estimates of basic reproduction numbers are very sensitive to noise in the testing rates as well as to changes in case definition. Moreover, in real populations, *ℛ**t* can move back above one quite easily in response to changes in the population protective behavior. In this paper, we approach the issue of local disease extinction by calculating the probability of extinction of the disease as a stochastic epidemic process. In the context of epidemics, the correct epidemiological terminology for the stochastic extinction is the local “elimination of infections”. In this paper, however, we shall continue to use the term extinction throughout, as it is shorter and more in line with the literature on stochastic processes. The probability of extinction is greater when the number of infected is low, so we calculate this probability assuming post-containment scenarios where cases have been dropped to very low numbers. We start by presenting the deterministic version of the model and derive its basic reproduction number. Then we derive an stochastic version of the same model and use it to calculate an analytical expression for the extinction probability. We conclude by looking at scenarios of local extinction and discussing how it applies to real scenarios, including also the impact of the fraction of the population already immunized upon the lifting of containment. ## 2. Methods ### 2.1. The SEIAHR model First, we describe the (deterministic) model used to represent COVID-19 dynamics. The Susceptible-Exposed-Infectious-Removed (SEIR) model is a classic model for diseases for which it is important to take into account an incubation period, and variants of it have been employed in numerous COVID-19 modelling studies [13, 14]. Here, we propose a variation of this model with added asymptomatic, hospitalized compartments and a quarantine mechanism (fig. 1). Ther is no explicit compartments for Quarantied an dead individuals as they are pure sink states that do not influence the main dynamics. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/17/2020.06.15.20132050/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/F1) Figure 1: Block diagram of the SEIAHR model. compartments *Q* and *death* are included for illustrative purposes only. The model is represented as system of ordinary differential equations: ![Formula][1] ![Formula][2] ![Formula][3] ![Formula][4] ![Formula][5] ![Formula][6] with *λ* = *β*(*I* + *A*) as the force of infection. State variables *S, E, I, A, R, H* represent the fraction of the population in each of the compartments, thus *S*(*t*) + *E*(*t*) + *I*(*t*) + *A*(*t*) + *H*(*t*) + *R*(*t*) = 1 at any time *t*. Quarantine enters the model through the parameter *χ* which can be taken as a constant or as a function of time, *χ*(*t*), that represents the modulation of the isolation policies. Quarantine works by blocking a fraction *χ* of the susceptibles from being exposed, i.e., taking part on disease transmission. Time-varying quarantine is achieved through multiplying *χ* by activation (eq. 2) and deactivation (eq. 3) functions: ![Formula][7] and ![Formula][8] where *s* and *e* are the start and end of the isolation period (*e > s*), respectively. A finite period *τ* = *e* − *s* of quarantine can be defined by the combined effect of both functions: ![Formula][9] Figure 2 illustrates the activation and deactivation of quarantine. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/17/2020.06.15.20132050/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/F2) Figure 2: *χ*(*t*) upon activation and deactivation of quarantine for *χ* = 0.3. Panel on the left shows each function separately(*χA*1(*t*) and *χD*1(*t*)) set to *t* = 1. The right panel shows a combination with activation on *t* = 1 and deactivation on *t* = 6 (*χ*(*t*) = *χA*1(*t*)*D*6 (*t*)). ### 2.2. Basic Reproduction Number The basic reproduction number for the SEIAHR model can be derived using the next generation matrix method[15], which we show in Appendix A. The expression for *ℛ* is ![Formula][10] From equation (5) we can obtain *ℛ**t* = *ℛ**S*(*t*) and another reproduction number, denoted by ℛ*c*, which we shall call *control reproduction number*, as it represents the average number of secondary cases infected by primary cases under some control scenario – i.e. whenever *χ >* 0. ### 2.3. Probability of extinction The question of the local extinction the disease can be more realistically addressed with a stochastic version of the SEIAHR model, where the combination of a discrete state and stochastic state transitions allow for actual extinctions to occur. The transition rates from the ODE model (Eqs 1) listed on table 1, can be used to build a continuous–time markov chain model where time is continuous but the state is discrete, allowing for a more realistic description of population changes over time. View this table: [Table 1:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/T1) Table 1: State transitions and rates for the stochastic SEIAHR model. In the state transition column, only the changing state variables are indicated. This model is a multivariate stochastic process *{S*(*t*), *E*(*t*), *I*(*t*), *A*(*t*), *H*(*t*) *}* where *R*(*t*) = *N* − (*S*(*t*) + *E*(*t*) + *I*(*t*) + *A*(*t*) + *H*(*t*)). We will leave the equation of R out, because it is decoupled from the rest of the system. A joint probability function is associated with the set of random state variables, *{S*(*t*), *E*(*t*), *I*(*t*), *A*(*t*), *H*(*t*)*}*, ![Formula][11] which leads to the Kolmogorov forward equation ![Formula][12] As a continuous–time branching process, the extinction threshold for the stochastic model is closely related to the corresponding one in the deterministic model but depends on the initial number of infectious individuals[16]. Based of the properties of this kind of stochastic processes, Whittle (1955) calculated the probability of extinction for the stochastic SIR model to be ![Graphic][13], where *i* is the initial number of infectious individuals[17]. We can apply the same technique described for the stochastic SEIR model by Allen and Lahodny [16] to derive the probability of extinction for the SEIAHR model. The analysis presented in Allen and Lahodny [16] assume proximity to the DFE with a large enough number of susceptibles and a small number of infectious individuals. Assuming a value for ℛ from other epidemics or estimated from initial exponential growth. For a realistic application to the COVID-19 epidemic at the moment *t* of the relaxation of population lockdown, we need to acknowledge that *S*(*t*) *< N*, whilst not knowing what the exact value of *S* at time *t* and thus the effective reproduction number at the time. Nevertheless, we know from equation (5) that ℛ*t* is a function of *S*(*t*). Therefore, we calculate ℙ for various values of *ℛ**t*. #### Analytical expression for ℙ The probability of extinction, ℙ, can be computed analytically following the derivation of probability-generating functions for the system of equations in (6), the details of which are given in Ap- pendix B. The probability of extinction is computed from the fixed points of the PGFs, *q*1, *q*2 and *q*3, which lie in (0, 1)3. With the fixed points in hand, we arrive at ![Formula][14] where *k**i* are the initial states *k*1 = *E*(0), *k*2 = *I*(0) and *k*3 = *A*(0). Here we denote the moment of relaxation of containment as *t* = 0. #### Numerical approximation to ℙ The expression in (7), while very useful, is derived from the assumption that the initial state is small compared to the size of the population and that *S*(0) ≈ *N*. Due to the possible deviations from the theoretical value of ℙ when *S*(0) *< N*, In the results, we always present estimates of ℙ , obtained by simulation as well. We can approximate ℙ for a given ℛ and an initial number of infected individuals (*E* + *I* + *A*). This can be accomplished by running a large number of simulations of the stochastic model and computing the fraction of simulations in which the virus is driven out of the population (*E* + *I* + *A* = 0) without first causing an outbreak. To facilitate setting up the simulation to specific *ℛ*, we can rewrite the force of infection by replacing *β* as a function of ℛ*R* (eq. A.8): ![Formula][15] ### 2.4. Effect of Increasing Seroprevalence After a first wave of infections a fraction of the population will become immunized against SARS-COV-2. This will have a protective effect on the population even if *S*(*t*) is still above the so-called “herd immunity” threshold of 1*/ℛ* . The resulting adjusted ℙ for different levels of population immunization can be determined through simulation of the stochastic SEIAHR model (eq. 6). ## 3. Results Figure 3 shows how imposing strong containment measures can reverse the growth of the number of cases. But a second wave starts immediately after the restrictions are relaxed. Figure 4 shows that in three out of ten post-containment simulations, the containment resulted in the extinction of the disease, with the other resulting on a second wave. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/17/2020.06.15.20132050/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/F3) Figure 3: Effect of containment, *χ* = .99 for 55 days starting on the 50th day of the epidemic. Simulation parameters were ϕ = 0.01, ℛ = 1.7, *ρ* = 0.21, *δ*= 0.04, *α* = 0.34, *µ* = 0.02, *p* = 0.76, *s* = 50, *e* = 105. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/17/2020.06.15.20132050/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/F4) Figure 4: Ten runs of the stochastic SEIAHR model with the same parameters used on fig. 3. In 3 out of 10 runs the containment could eliminate the disease but the remaining 7 gave origin to a second wave. All simulations had *I*(0) = 2 and ℛ = 1.7, on a population of 5000. Depending on the value of ℛ*t* at the moment the quarantine is lifted and the number of remaining active cases, the probability of extinction can favor the stochastic extinction of the disease (Fig 5). ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/17/2020.06.15.20132050/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/F5) Figure 5: Probability of local extinction, ℙ as a function of ℛ*t* and *I* .(a)for the SIR model and (b) for the SEIAHR model. Tables 2 and 3 contain the values of ℙ for different scenarios of number of infectious, with high and low ℛ respectively. View this table: [Table 2:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/T2) Table 2: Probabilities of extinction for ℛ*t* = 2.5. Approximate ℙ was calculated from a set of 10000 runs of the stochastic SEIAHR model with different number of initial infectious individuals(*I* ). View this table: [Table 3:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/T3) Table 3: Probabilities of extinction for ℛ*t* = 1.7. Approximate ℙ was calculated from a set of 10000 runs of the stochastic SEIAHR model with different number of initial infectious individuals(*I*). View this table: [Table 4:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/T4) Table 4: Probabilities of extinction for ℛ*t* = 1.1. Approximate ℙ was calculated from a set of 10000 runs of the stochastic SEIAHR model with different number of initial infectious individuals(*I*). We also estimated the ℙ adjusted for different fractions of immune individuals (*R*) in population from 10000 simulations of the stochastic model. The results are presented in table 5. View this table: [Table 5:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/T5) Table 5: Adjustedℙ for different levels seroconversion (*R*(0)) of the population. Results are the fraction of stochastic extinctions in 10000 simulations, with *ℛ* = 1.1. Table 6 shows the expected time in days until extinction under different levels of population immunity and post-containment number of infected. View this table: [Table 6:](http://medrxiv.org/content/early/2020/06/17/2020.06.15.20132050/T6) Table 6: Expected time to elimination of infections, in days, counting from the end of containment. Time is expressed as mean and 95% interval. Values where estimated from a 10000 runs of the stochastic SEIAHR model and rounded to the nearest integer number of days. ## 4. Discussion Most countries who managed to control the first wave of COVID-19 did it by means of imposed restrictions on human-to-human contact. Be it through quarantine, tracing and isolation, use of masks or a combination of social distancing tactics, they managed to minimize contacts and thus transmission. Here we have introduced a transmission model (SEIAHR) in both deterministic and stochastic formulation, which includes the isolation of susceptibles as a means of reducing transmission as well as different levels of infectious individuals. From these models (deterministic and stochastic versions) it is possible to explore the consequences of social distancing as well as to calculate the probability of a second wave upon the suspension of distancing behavior. Figure 3 shows how imposing strong social distancing can interrupt the transmission, bringing the prevalence to near zero. However, if steps are not taken to permanently change the way people interact in their daily routine, at work, school, public transport, etc. The reproduction number upon lifting of the restrictive regulations, will still be substantially higher than one and will drive a powerful second wave if the population is still far from herd immunity conditions. In the deterministic version of SEIAHR, a second wave will always happen due to the asymptotic way that the number of infectious approach zero. Treating the epidemic process as the stochastic process that it actually is, we can see that the probability of local extinction post-reactivation is quite substantial (fig 5). In the stochastic SEIAHR, extinction events will happen with slightly lower probabilities than those of a SIR model justifying using more detailed model to study this problem. Figure 4 shows stochastic extinction taking place in three out of ten runs, with ℛ*t* = 1.7 and *I* = 2. The Stochastic model allows us to compute the probabilities of eliminating local transmission under various scenarios, and can be a useful tool for planning when to lift restrictions to human mobility and interaction. Though the probabilities of elimination of local infections seem rather low on the basic scenarios described on tables 2, 3 and 4, these do not tell us the whole story. If one takes into account the immunity acquired by the population during the containment period, one can see from the results in table 5, that even if a location is still far from achieving herd immunity, any acquired immunity will greatly improve the chances local extinction substantially. Communities that managed to contain the disease and bring it to the brink of extinction with severe economic impact, need to know how likely they are in succeeding in their fight against the disease as they return to “normal” social and economic activities. Our results also reinforce the need to run seroprevalence surveys previous to the reopening so that the probability to eliminate local infections e properly adjusted for the context of each locality. Here we explored how to improve the chances for local transmission elimination but it must be kept in mind that the ℛ*t* must be kept low (through the use of masks and social distancing) in the post-containment period to guarantee elimination and even when properly executed the time to elimination can vary from weeks to months, as shown in table 6. ## Data Availability Upon publication ## Appendix A. Derivation of the basic reproduction number In order to use the next-generation method, we start by identifying all *m* compartments containing infected individuals. In this case they are *E, I* and *A* and *m* = 3. Let ℱ be the vector of the rates of appearance of new infections in the three infectious compartments: ![Formula][16] Next, we let *𝒱**i*(*x*)− be rate infectious individuals leave compartment *i*, and *𝒱**i*(*x*)+ be the rate with which they enter compartment *i*: ![Formula][17] and ![Formula][18] After defining these we can calculate *𝒱* (*x*) = *V𝒱**i*(*x*)− − *𝒱**i*(*x*)+: ![Formula][19] Now, if we let *x* = *{E, I, A}*and *x*0 be the Disease-free equilibrium (DFE), we can define ![Formula][20] and ![Formula][21] The next generation matrix is given by *FV* −1: ![Formula][22] The spectral radius of 𝕂 at the DFE (when *S*(0) ≈ 1), is the basic reproduction number of the model, *ℛ* = *ρ*(*FV* −1), ![Formula][23] which, after simplification, gives equation (5). ## Appendix B. Probability generating functions Following [16], we derive the probability-generating functions (PGF) for each infectious compartment. Starting with *I**i*(0), the probability of an infected individual in state *i* producing offspring of type *j* given that *I**j*(0) can be obtained from ![Formula][24] The desired probabilities can be obtained by differentiating the PGF with respect to *z**i* and setting all ***z*** to 1. Notice *f**i* has a fixed point at *z*1 = … = *z**k* = 1. Computing the relevant probabilities is straightforward by keeping track of the possible transitions (given in Table 1) and considering that only one transition may occur in a given time interval Δ*t*. First, when *E*(0) = 1, *I*(0) = 0 and *A*(0) = 0, we have ![Formula][25] and when *E*(0) = 0, *I*(0) = 1 and *A*(0) = 0 the PGF is ![Formula][26] Finally, for the case *E*(0) = 0, *I*(0) = 0 and *A*(0) = 1: ![Formula][27] From the PGFs, we can obtain a matrix 𝕄 whose entries ![Graphic][28] are the expected number of offspring in state *j* from an individual in state *i*. Given the PGFs above, we arrive at ![Formula][29] To obtain the probability of extinction, ℙ, we need to find the fixed points of the PGFs, i.e. solutions to equations of the form fi(q1, q2, q3) = qi, qi ∈ (0, 1). After some tedious algebra, we arrive at ![Formula][30] ![Formula][31] ![Formula][32] With ![Graphic][33]. * Received June 15, 2020. * Revision received June 15, 2020. * Accepted June 17, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].List of epidemics, 2020. URL: [https://en.wikipedia.org/wiki/List\_of\_epidemics](https://en.wikipedia.org/wiki/List_of_epidemics), page Version ID: 962716222. 2. [2]. E. Dong, H. Du, L. Gardner, An interactive web-based dashboard to track covid-19 in real time, The Lancet infectious diseases 20 (2020) 533–534. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30120-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F17%2F2020.06.15.20132050.atom) 3. [3]. B. F. Maier, D. Brockmann, Effective containment explains subexponential growth in recent confirmed covid-19 cases in china, Science 368 (2020) 742–746. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5Mi83NDIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNi8xNy8yMDIwLjA2LjE1LjIwMTMyMDUwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 4. [4]. S. Zhang, Z. Wang, R. Chang, H. Wang, C. Xu, X. Yu, L. Tsamlag, Y. Dong, H. Wang, Y. Cai, Covid-19 containment: China provides important lessons for global response, Frontiers of Medicine (2020) 1–5. 5. [5]. C. J. Wang, C. Y. Ng, R. H. Brook, Response to covid-19 in taiwan: big data analytics, new technology, and proactiv testing, Jama 323 (2020) 1341–1342. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F17%2F2020.06.15.20132050.atom) 6. [6]. N. Fernandes, Economic effects of coronavirus outbreak (covid-19) on the world economy, Available at SSRN 3557504 (2020). 7. [7]. S. Reed, J. M. Gonzalez, R. Johnson, Willingness to accept tradeoffs among covid-19 cases, social-distancing restrictions, and economic impact: A nationwide us study, medRxiv (2020). 8. [8]. S. Xu, Y. Li, Beware of the second wave of covid-19, The Lancet 395 (2020) 1321–1322. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)30845-x&link_type=DOI) 9. [9]. K. Leung, J. T. Wu, D. Liu, G. M. Leung, First-wave covid-19 transmissibility and severity in china outside hubei after control measures, and second-wave scenario planning: a modelling impact assessment, The Lancet (2020). 10. [10]. M. Gatto, E. Bertuzzo, L. Mari, S. Miccoli, L. Carraro, R. Casagrandi, A. Rinaldo, Spread and dynamics of the covid-19 epidemic in italy: Effects of emergency containment measures, Proceedings of the National Academy of Sciences 117 (2020) 10484–10491. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzE5LzEwNDg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDYvMTcvMjAyMC4wNi4xNS4yMDEzMjA1MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. [11]. Z. Yang, Z. Zeng, K. Wang, S.-S. Wong, W. Liang, M. Zanin, P. Liu, X. Cao, Z. Gao, Z. Mai, et al., Modified seir and ai prediction of the epidemics trend of covid-19 in china under public health interventions, Journal of Thoracic Disease 12 (2020) 165. 12. [12]. A. J. Kucharski, T. W. Russell, C. Diamond, Y. Liu, J. Edmunds, S. Funk, R. M. Eggo, F. Sun, M. Jit, J. D. Munday, et al., Early dynamics of transmission and control of covid-19: a mathematical modelling study, The lancet infectious diseases (2020). 13. [13]. S. M. Kissler, C. Tedijanto, E. Goldstein, Y. H. Grad, M. Lipsitch, Projecting the transmission dynamics of sars-cov-2 through the postpandemic period, Science 368 (2020) 860–868. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5My84NjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNi8xNy8yMDIwLjA2LjE1LjIwMTMyMDUwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 14. [14]. K. Prem, Y. Liu, T. W. Russell, A. J. Kucharski, R. M. Eggo, N. Davies, S. Flasche, S. Clifford, C. A. Pearson, J. D. Munday, et al., The effect of control strategies to reduce social mixing on outcomes of the covid-19 epidemic in wuhan, china: a modelling study, The Lancet Public Health (2020). 15. [15]. P. van den Driessche, J. Watmough, Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission, Mathematical Biosciences 180 (2002) 29–48. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0025-5564(02)00108-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12387915&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F17%2F2020.06.15.20132050.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179220600004&link_type=ISI) 16. [16]. L. J. Allen, G. E. Lahodny Jr, Extinction thresholds in deterministic and stochastic epidemic models, Journal of Biological Dynamics 6 (2012) 590–611. 17. [17]. P. Whittle, The outcome of a stochastic epidemic—a note on Bailey’s paper, Biometrika 42 (1955) 116–122. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/42.1-2.116&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1955XF69000009&link_type=ISI) [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/graphic-8.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-10.gif [10]: /embed/graphic-12.gif [11]: /embed/graphic-14.gif [12]: /embed/graphic-15.gif [13]: /embed/inline-graphic-1.gif [14]: /embed/graphic-16.gif [15]: /embed/graphic-17.gif [16]: /embed/graphic-26.gif [17]: /embed/graphic-27.gif [18]: /embed/graphic-28.gif [19]: /embed/graphic-29.gif [20]: /embed/graphic-30.gif [21]: /embed/graphic-31.gif [22]: /embed/graphic-32.gif [23]: /embed/graphic-33.gif [24]: /embed/graphic-34.gif [25]: /embed/graphic-35.gif [26]: /embed/graphic-36.gif [27]: /embed/graphic-37.gif [28]: /embed/inline-graphic-2.gif [29]: /embed/graphic-38.gif [30]: /embed/graphic-39.gif [31]: /embed/graphic-40.gif [32]: /embed/graphic-41.gif [33]: /embed/inline-graphic-3.gif