Dynamical SPQEIR model assesses the effectiveness of non-pharmaceutical interventions against COVID-19 epidemic outbreaks ========================================================================================================================= * Daniele Proverbio * Françoise Kemp * Stefano Magni * Andreas Husch * Atte Aalto * Laurent Mombaerts * Alexander Skupin * Jorge Gonçalves * Jose Ameijeiras-Alonso * Christophe Ley ## Abstract Against the current COVID-19 pandemic, governments worldwide have devised a variety of non-pharmaceutical interventions to suppress it, but the efficacy of distinct measures is not yet well quantified. In this paper, we propose a novel tool to achieve this quantification. In fact, this paper develops a new extended epidemic SEIR model, informed by a socio-political classification of different interventions, to assess the value of several suppression approaches. First, we inquire the conceptual effect of suppression parameters on the infection curve. Then, we illustrate the potential of our model on data from a number of countries, to perform cross-country comparisons. This gives information on the best synergies of interventions to control epidemic outbreaks while minimising impact on socio-economic needs. For instance, our results suggest that, while rapid and strong lock-down is an effective pandemic suppression measure, a combination of social distancing and contact tracing can achieve similar suppression synergistically. This quantitative understanding will support the establishment of mid- and long-term interventions, to prepare containment strategies against further outbreaks. This paper also provides an online tool that allows researchers and decision makers to interactively simulate diverse scenarios with our model. Keywords * COVID-19 * Cross-country comparison * Risk assessment * Epidemiological modelling * Non-pharmaceutical interventions * Public health ## 1 Introduction The current global COVID-19 epidemic has led to significant impairments of public life world-wide. To suppress the spread of the virus and to prevent dramatic situations in the healthcare systems, many countries have implemented a combination of rigorous measures like lock-down, isolation of symptomatic cases and the tracing, testing, and quarantine of their contacts. In order to gain information about the efficacy of such measures, a quantitative understanding of their impact is necessary. This can be based on statistical methods [1] and on epidemiological models [2]. While statistical methods allow for accurate characterization of the population’s health state, epidemiological modeling can provide more detailed mechanisms for the epidemic dynamics and allow investigating how epidemics will develop under different assumptions. Preliminary efforts have been made to quantify the contribution of different policy interventions [3], but these rely on complex models based on a number of assumptions. We base our study on a classical SEIR-like epidemiological model. SEIR models are minimal mechanistic models that consider individuals transitioning through Susceptible→ Exposed → Infectious → Removed state during the epidemics [4]. The essential control parameter is the basic reproduction number *R* [5], that worldwide non-pharmaceutical suppression strategies aim at reducing below the threshold value 1. Building on this, we incorporate additional compartments reflecting different categories of intervention strategies, identified by socio-political studies [6]. In particular, the model focuses on three main suppression programs: social distancing (lowering the rate of social contacts), active protection (decreasing the number of susceptible people), and active removal of latent asymptomatic carriers [7]. This study investigates how these programs achieve repression both individually and combined, first conceptually and then by cross-country comparison. This information can supply Government decisions, helping to avoid overloading the healthcare system and to minimise stressing the economic system (due to prolonged lock-down). We expect our model, together with its interactive online tool, to contribute to crucial tasks of decision making and to prepare containment strategies against further outbreaks. ## 2 Methods This study links policy measures to epidemiological modelling and uses the developed model to quantitatively assess the efficacy of different interventions in six countries. In this section, we illustrate the modelling choice and the use of data. ### 2.1 The classical SEIR model SEIR models are continuous-time, mass conservative compartment-based models of infectious diseases [4, 8]. They assume homogeneous propagation media (or fully connected graphs) and focus on the evolution of mean properties of the closed system. All of these models, from the more conceptual to more realistic versions, e.g. SEIR with delay [9], spatial coupling [10], or individual-based models [11], are classical and widely used tools to investigate the principal mechanisms governing the spread of infections and their dynamics. Main compartments of SEIR models (see Fig. 1, framed) are: susceptible S (the pool of individuals likely to be infected), exposed E (corresponding to latent carriers of the infection), infectious I (individuals having developed the disease and being contagious) and removed R (those that have processed the disease, being either recovered or dead). The model’s default parameters are the average contact rate *β*, the inverse of mean incubation period *α* and the inverse of mean contagious period *γ*. When focusing on infection dynamics rather than patients’ fate, the latter combines recovery and death rate [12]. From these parameters, epidemiologists calculate the “basic reproduction number” *R* = *β/γ* [13] at the epidemic beginning. During the epidemic progression, isolation after diagnosis, vaccination campaigns and active suppression measures are in action. Hence, we speak of “effective reproduction number” ![Graphic][1] [14]. ![Fig. 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F1.medium.gif) [Fig. 1](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F1) Fig. 1 Scheme of the SPQEIR model. The basic SEIR model (framed blue blocks) is extended by the red blocks to the SPQEIR model. Parameters that are linked to repression strategies are shown in red. Interpretation and values of parameters are given in Table 2. ### 2.2 Data and analyzed countries Governments worldwide have issued a number of social measures, including those for public health safeguard, economic support, movement restriction and non-pharmaceutical interventions to hamper disease spreading. Scholars from political sciences and sociology have recorded and classified such measures [15, 16]. Among the resources listed on the World Health Organization “Tracking Public Health and Policy Measures” [6], we used information from the ACAPS database [17] that contains a curated categorization of policy measures. ACAPS is an independent, non-profit information provider helping humanitarian actors respond more effectively to disasters. The ACAPS analysis team has aggregated and classified interventions from different sources (media, governments and international organizations), for all countries and in time. Suppression measures against the epidemic are classified under “Movement restrictions”, “Lock-down”, “Social Distancing” and “Monitoring and Surveillance”. Our modelling choice is based on these categories, which are reflected by additional compartments to the classical SEIR model (see next section). Epidemiological data for all selected countries and regions were obtained from the COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University [18]. The data are from 22 Jan 2020 to 08 July 2020. Lombardy data were obtained from the Protezione Civile Italiana data repository “Dati COVID-19 Italia” [19], from 22 Feb 2020 to 08 July 2020. This study analyses the effect of suppression measures in flattening the curve. Despite having a precise starting date, such measures take some days to be fully effective. We estimate an average delay using the Google Mobility Reports [20, 21] for the selected countries. Google provides changes in mobility with respect to a monthly baseline, w.r.t. 6 locations: Retail & Recreation, Grocery & Pharmacy, Transit stations, Workplaces, Residential, Parks. We average the decrease in mobility at the first four locations (corresponding to those where social mixing happens more frequently [22]) to get a proxy of the time needed for hard lock-down to be fully effective (cf. Fig. 4c). ### 2.3 The extended SPQEIR model to reflect suppression strategies SEIR models reproduce the typical bell-shaped epidemic curves for the number of infected (and still infectious) people. These quantify the main stressors for both the health system, i.e. the peak of the curve, and the economic system, i.e. the time 𝕋 passed until no new infections occur. Mainstream suppression measures against the epidemic aim at flattening the curve of new infections [7]. However, the classical SEIR model is not granular enough to investigate suppression measures when they need to be considered or should be sequentially reduced if already in place. Therefore, we extend the classical SEIR model as in Fig. 1 (red insertions) into the SPQEIR model, to reflect the intervention categories described above. The model can be summarized as follows: * The classical blocks S, E, I, R are maintained; * A social distancing parameter *ρ* is included to tune the contact rate *β*; * Two new compartments are introduced where: * – Protected P includes individuals that are isolated from the virus through lock-down, thus reducing the susceptible pool; * – Quarantined Q describes latent carriers that are identified and quarantined after monitoring and tracing. We do not introduce a second quarantined state for isolation of confirmed cases after the Infectious state [23, 2] but consider this together with the Removed state (see Liu et al. [24] and references therein). Quarantining infected symptomatic patients is a necessary first step in every epidemic [25]. An additional link from Q to R, even though realistic, is neglected as both compartments are already outside the “contagion system” and would therefore be redundant from the perspective of evolution of the infection. In general, protected individuals can get back to the pool of susceptible after a while, but here we neglect this transition, to focus on simulating repression programs alone. Long-term predictions could be modelled even more realistically by considering such link, that would lead to an additional parameter to be estimated and is beyond the scope of the present paper. The model has in total 6 parameters. Three of them (*β, α, γ* introduced in Fig. 1) are based on the classical SEIR model. The new parameters *ρ, µ, χ* account for alternative repression programs (see Table 2 for details). Commonly, social distancing is modeled by the parameter *ρ*. It tunes the contact rate parameter *β*, resulting in the effective reproduction number ![Graphic][2]. This occurs in a closed-system setting where all individuals belong to the susceptible pool, but interact less intensively with each other. The parameter *µ* stably decreases the susceptible population by introducing an active protection rate. This accounts for improvements of public health, e.g. stricter lock-down of communities, or physical reduction of a country’s population like reduced commuters’ activity. This changes the effective reproduction number into ![Graphic][3] with *T* being the number of days the measures are effective [26]. The parameter *χ* introduces an active removal rate of latent carriers. Intensive contact tracing and improved methods to detect asymptomatic latent carriers may enhance the removal of exposed subjects from the infectious network. Following earlier works [27, 28] and adjusting the current parameters, ![Graphic][4] can be then expressed as ![Graphic][5]. Parameter values that are not related to suppression strategies are set from COVID-19 epidemic literature [24, 29]. We use mean values as the main focus of the present model lies on sensitivity analysis of suppression parameters. Our model can be further extended by time dependent parameters [25]. Default values for suppression parameters are *{ρ, µ, χ}* = {1, 0, 0}, corresponding to the classical SEIR model. The dynamics of our SPQEIR model is described by the following system of differential equations: ![Formula][6] with conservation of the total number of individuals, meaning ![Graphic][7] with N = S + E + I + R + P + Q. As value for the qualitative study, we used N = 10,000. For the cross-country assessment, N is adjusted to true population values for each country. Overall, the effective reproductive number becomes ![Formula][8] with *T* being the number of days that the measures leading to compartment P are active. Suppression measures are initiated several days after the first infection case. Hence, we activate non-default parameter values after a delay *τ*. In conceptual simulations, we set it arbitrarily without loss of generality. For data fitting, we fit and compare it to the official date when measures are initialized (cf. Table 1). To integrate the model numerically, we use the *odeint* function from *scipy*.*integrate* Python library. View this table: [Table 1:](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/T1) Table 1: Test countries, with measures implemented (following the ACAPS database [17]), corresponding parameter in our SPQEIR model and starting date. For Lombardy, we used the Italian official date for lock-down. Differently from other countries, Ireland issued measures on two different dates; we use this case to compare social distancing and lock-down effect in a single country. We also report the (rounded) population of each country. ### 2.4 Model fitting We fit the model to the official number of currently infected (active) cases, for each considered country. Model fitting to the infectious curves is performed in two steps, using the parameters known to be active (cf. Table 1). First, we estimate the “model consistent” date of first infection, so that the simulated curve matches the reported data of active infections. This initial step corresponds to setting the time initial conditions of the SEIR model [2, 26]. The fitting is performed with default parameter values, on a subset of data corresponding to the first outbreak, from first case until when measures are implemented (cf. Table 1). We use a grid search method for least squares, sufficient to fit a single parameter: ![Formula][9] where *t* is the “model consistent” estimated date of first infection, *t**m* refers to the date measures are implemented, ![Graphic][10] and *x* are respectively reported and model-predicted data, and *n* is the number of points between *t* and *t**m*. The second step estimates the suppression parameters that yield the best fitting of the simulated SPQEIR curve on reported data, during the first phase with implemented measures. This period is identified between the starting date *t**m* (also included in the fitting) and the phase-out date *t**p*, cf. Table 1. Holding the epidemic parameters to literature values to achieve cross-country comparison on intervention parameters alone, the fitting is performed for a set of suppression parameters relative to each country, as reported by policy databases (cf. Table 1). The fit is performed with *lmfit* Python library. When more parameters are involved at once, we also perform a comparative information analysis between our extended model and the simplest SEIR that lumps parameters under a single “social distancing” *ρ*. This shows what we gain in distinguishing the intervention parameters, not only in terms of improved interpretation, but also in terms of model fitting. We employ a reduced *χ*2 metric to evaluate the goodness-of-fit for both models, considering the degrees of freedom [31]: ![Formula][11] where *n*′ is the number of data points until phase-out, *k* is the number of parameters in the model, *y**j* are estimated values and ![Graphic][12] the expected ones. The lower the ![Graphic][13] value, the better the SPQEIR model fits the data w.r.t. the simplest “social distancing” SEIR model. ## 3 Results We first focus on the conceptual analysis of the effect of suppression interventions, initially for single measures (social distancing, active protection and active quarantining) and subsequently for a number of synergistic approaches. In particular, we study how crucial quantities, namely the infectious peak height and time to zero infectious, depend on suppression parameters and affect ![Graphic][14]. We define 𝕋 as the time when there are less than 0.5 individuals in the I compartment. This because ODE models approximate discrete quantities with continuous variables. Finally, we perform model fitting and intervention assessment over a set of countries. This provides quantitative outputs about the effectiveness of real measures, informing about the synergies applied and enabling cross comparison. ### 3.1 Mathematical analysis of single suppression measures #### 3.1.1 Only social distancing The parameter *ρ* captures social distancing effects, taking values in the interval [0, 1], where 0 indicates no contacts among individuals while 1 is equivalent to no actions taken. Without loss of generality, simulations consider a delay of 10 days from the first infection to the time social distancing is initiated. Fig. 2 reports simulation results. The curve of infectious is progressively flattened by social distancing (2a) and its peak suppressed (2b). However, the eradication time gets delayed for decreasing *ρ*, until a threshold yielding a disease-free equilibrium rapidly (2c). In this case, the critical value for *ρ* is 0.4, leading to ![Graphic][15]. However, we notice that values of *ρ ≃* 0.3 or lower are more effective in suppressing the epidemic faster. This is in line with early findings, suggesting that epidemic control with social distancing alone should be done “well or not at all” [32]. ![Fig. 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F2.medium.gif) [Fig. 2](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F2) Fig. 2 (a) Effects of social distancing on the epidemic curve. The grey area indicates when measures are not yet in place. (b) The peak is progressively flattened until a suppression is reached for sufficiently small *ρ*. For these settings, the critical value for *ρ* is 0.4 (it pushes ![Graphic][16] below 1). (c) Unless *ρ* is small enough, stronger measures of this kind might delay the suppression time 𝕋 of the epidemic. #### 3.1.2 Only active protection As above, our simulations take into account 10 days delay from the first infection to the initiation of active protection. The range of *µ* is only up to values similar to those measured in China [26]. Higher values are considered for step-wise hard lock-down (see below). The results are reported in Fig. 3. We see that small precautions can make an initial difference, but then the effects saturate (Fig. 3a,b). The time to zero infectious is decreased with higher values of active protection (Fig. 3a,b). In particular, *µ* = 0.01 *d*−1 suppresses the epidemics in about 6 months by protecting 70% of the population. Higher values of *µ* achieve suppression faster, while protecting almost 100% of the population. If protection is mostly achieved through isolation, this is unrealistic. However, this models effective vaccination strategies. ![Fig. 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F3.medium.gif) [Fig. 3](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F3) Fig. 3 (a) Effects of active protection on the infectious curve. The grey area indicates when measures are not yet in place. *µ* is expressed in *d*−1. (b) Dependency of peak height on *µ*: the peak is rapidly flattened for increasing *µ*, then it is smoothly reduced for higher parameter values. (c) High *µ* values are effective in anticipating the complete eradication of the epidemic, but require protecting more than 90% of the population. We also consider hard lock-down strategies which isolate many people at once [33]. This corresponds to reducing S to a relatively small fraction rapidly. Since *µ* is a rate, we mimic a step-wise hard lockdown by setting a high value to *µ*, but its effect only lasts for a short period of time, see Fig. 4b. We thus use the notation *µ**ld*. In the figure, an example shows how to rapidly protect about 68% of the population with a step-wise *µ**ld* function. In particular, we use an average four-days long step-wise *µ**ld* function (Fig. 4b) to mimic the rapid, but not abrupt, change in mobility observed in many countries by Google Mobility Reports [21] (Fig. 4c). Lock-down effects are reported in Fig. 4a: a hard lock-down is effective in suppressing the epidemic curve and in lowering the eradication time. ![Fig. 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F4.medium.gif) [Fig. 4](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F4) Fig. 4 (a) Flattening the infectious curve by hard lock-down. Rapidly isolating a large population fraction is effective in suppressing the epidemic spreading. (b) Modeling hard lock-down: high *µ**ld* (orange) is active for four days to isolate and protect a large population fraction rapidly (blue). As an example, we show *µ**ld* = 0.28 *d*−1 if *t* ∈ [10, 4]. It results in protecting about 68% of the population in two days. Higher values, e.g. *µ**ld* = 0.65 *d*−1 would protect 93% of the population at once. (c) Google Mobility Report visualization [21] for analysed countries, around the date of measures setting. Each line reports the mean in mobility change across Retail & Recreation, Grocery & Pharmacy, Transit stations, and Workplaces, around the date of implementation of the measures. A minimum of 4 days (from top to bottom of steep decrease) is required for measures to be fully effective. Abbreviation explanations: AT = Austria, CH = Switzerland, DK = Denmark, IL = Israel, IR = Ireland, LO = Lombardy. #### 3.1.3 Only active quarantining The simulations in this part are based on realistic assumptions: testing a person is effective only after a few days that that person has been exposed (to have a viral charge that is detectable). This induces a maximal quarantining rate *θ*, which we set *θ* = 0.33 *d*−1 as testing is often considered effective after about 3 days from contagion [34]. Therefore, we get the active quarantining rate *χ* = *χ*′ *θ*, where *χ*′ is a tuning parameter associated e.g. to contact tracing. As *θ* is fixed, we focus our analysis on *χ*′. As above, we also assume that testing starts after the epidemic is seen in the population, i.e. some infectious are identified with the usual 10 days delay in the activation of measures. The corresponding results are reported in Fig. 5. The curve is progressively flattened by latent carriers quarantining and its peak suppressed, but the eradication time gets delayed for increasing *χ*′. This happens until a threshold value of ![Graphic][17] that pushes ![Graphic][18] below 1. This value holds if we accept a strategy based on testing, with *θ* = 0.33. If preventive quarantine of suspected cases does not need testing (for instance, it is achieved by contact tracing apps), the critical *χ*′ value could be drastically lower. In particular, ![Graphic][19] if *θ* = 1 *d*−1, i.e. latent carriers are quarantined the day after a contact. ![Fig. 5](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F5.medium.gif) [Fig. 5](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F5) Fig. 5 (a) Effects of active latent carriers quarantining on the epidemic curve. The grey area indicates when measures are not yet in place. (b) The peak is progressively flattened until a disease-free equilibrium is reached for sufficiently large *χ*. (c) Unless *χ*′ is large enough, stronger measures of this kind might delay the complete eradication of the epidemic. Note that the critical *χ*′ can be lowered for higher *θ*, e.g. if preventive quarantine does not wait for a positive test. The parameter *χ*′ tunes the rate of removing latent carriers. Hence, it combines tracing and testing capacities, i.e. probability of finding latent carriers (*P**find*) and probability that their tests are positive (*P*+). The latter depends on the false negative rate *δ*− as ![Formula][20] So, *χ*′ = *P**find* *P*+. Hence, suppressing the peak of infectious requires an adequate balance of accurate tests and good tracing success as reported in Fig. 6. Further quantifying the latter would drastically improve our understanding of the current capabilities and of bottlenecks, towards a more comprehensive feasibility analysis. ![Fig. 6](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F6.medium.gif) [Fig. 6](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F6) Fig. 6 Assessing the impact of *P**find* and *P*+ on the peak of infectious separately. This way, we separate the contribution of those factors to look at resources needed from different fields, e.g. network engineering or wet lab biology. Solutions to boost the testing capacity like [35] could impact both terms. ### 3.2 Synergistic scenarios Fully enhanced active quarantining and active protection might not be always feasible, e.g. because of limited resources, technological limitations or welfare restrictions. Therefore a synergistic approach is very attractive as it can flatten the curve. This section shows a number of possible synergies, concentrating as before on abstract scenarios to investigate the effect of combining different suppression programs. As case studies, we consider the 6 synergistic scenarios listed below. Parameters are set without being specific to real measures taken: their value is so far conceptual and meaningful when compared across scenarios. Just like above, we consider a 10 days delay from the first infection to issuing measures; as suggested in other studies [36], delaying action could worsen the situation. To differentiate between a rapid isolation and a constant protection, we introduce *µ**ld* (for “hard lock-down strategies”, see Section 3.1.2) separated from *µ*. To get ![Graphic][21], we follow Eq. 1, considering *χ* = *χ*′ *θ* as in Section 3.1.3 and *T* = 4 (along the steep decay after measures are in place, in a best case scenario). Our scenarios are the following: 1. Many European countries opted for a lock-down strategy. A quite large fraction of the population was isolated, individuals were recommended to self-quarantine in case of suspected positiveness, social distancing got mandatory but was sometimes not fully followed, masks and sprays were suggested for protection. So, we set an initial “hard lock-down” *µ**ld* = 0.12 *d*−1 to protect around 38% of the population quickly. Then we chose *ρ* = 0.7, *χ*′ = 0.15 and *µ* = 0.008 *d*−1. This yields ![Graphic][22]. 2. An alternative procedure is to rapidly protect only the population fraction at high risk (*µ**ld* = 0.06 *d*−1, driving 15% of initial S to P). Then, we assume an improvement in individual safety giving *µ* = 0.01 *d*−1. Social distancing is relaxed (*ρ* = 0.8) but latent carrier quarantine is enforced (*χ*′ = 0.5). This gives ![Graphic][23] 3. In case preventive quarantine of latent carriers is not greatly effective (*χ*′ = 0.1), and in case of low protection rate and scarce isolation (*µ* = 0.004 *d*−1, *µ**ld* = 0.08 *d*−1), we rise social distancing for all individuals doing business as usual (*ρ* = 0.5). In this case, ![Graphic][24]. 4. If there are no safety devices that provide an adequate protection (*µ* = 0 *d*−1), we set *ρ* = 0.45, *µ**ld* = 0.2 *d*−1 and *χ*′ = 0.2 to get ![Graphic][25]. 5. This case has higher ![Graphic][26] than the previous ones, namely ![Graphic][27]. The corresponding parameters are *µ**ld* = 0.1 *d*−1, *µ* = 0.002 *d*−1, *ρ* = 0.7, *χ*′ = 0.1. 6. Finally, we consider “draconian” measures such that ![Graphic][28] only through isolation and massive latent carriers quarantining. So, *µ**ld* = 0.6 *d*−1 and *χ*′ = 0.3 while *ρ* = 1 and *µ* = 0 *d*−1. Simulation results are reported in Fig. 7. Different synergies lead to different timing, even though the peak is contained similarly (Fig. 7a). This has an impact on the cumulative number of cases (Fig. 7b) that will be reflected on the death toll. This holds even when the ![Graphic][29] values are very close, as in scenarios 1 to 4. Focusing on scenarios 2 and 3, we notice that prevention measures and latent quarantine accelerate the suppression, even when isolating only vulnerable people. This achieves similar effects as strong social distancing. In addition, active protective measures with relatively low values further concur in suppressing the peak. This finding asks for rapid assessment of masks and sanitising routines. ![Fig. 7](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F7.medium.gif) [Fig. 7](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F7) Fig. 7 Simulations of the 6 synergistic scenarios. (a) Curves of infectious Individuals, (b) Cumulative cases. The grey area indicates when measures are not yet in place. It is evident that scenarios leading to the same ![Graphic][30] could show different patterns and suppression timing. (c) Distribution of times to zero infections 𝕋 for different scenarios. Overall, the strength of suppression measures influences how and how fast the epidemic is flattened. *µ**ld* mostly governs the peak height after measures are implemented, *ρ* mainly tunes the curve steepness together with *µ*, while *χ* shifts the decaying slope up and down. Overall, a ![Graphic][31] suffices to avoid breakdown of the health system, but its effects could be too slow for the economic system. Decreasing its value with synergistic interventions could speed up epidemic suppression. A careful assessment of measures’ strength is thus recommended for cross-country comparison. ### 3.3 Model fitting and interventions assessment #### 3.3.1 Model fitting As described in the Methods section, we first estimate the “model consistent” date of first infection *t*, that is, the temporal initial condition for the SPQEIR model. We do not claim that this is the true date of first infection in a country. On the contrary, it is the starting date of infections in case of homogeneous transmission, under the assumption of no superspreading events [37], and with the hypothesis of coherent *R* (cf. Table 2). During the second fitting step, we also estimate the date at which suppression measures start having effect on the infectious curve, *t**m*. Comparing *t**m* with official intervention dates from Table 1, we notice that about 8 days are necessary to register lock-down effects. This is consistent to early findings on lock-down effectiveness [38]. Estimated dates are reported in Table 3. View this table: [Table 2:](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/T2) Table 2: SPQEIR model parameters with their standard values for the COVID-19 pandemic from literature [24, 30]. Here “*d*” stands for days. View this table: [Table 3:](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/T3) Table 3: Dates of official detection of first COVID-19 case [18], estimated dates for first infection *t* (according to Eq. 2) and date at which measures start being effective *t**m*, per country. Although consistent with recent literature [26] that suggests the first infection happened weeks before the first official detection, these retrospective dates should be interpreted under the model assumptions. Then, we fit suppression parameters to data affected by policy measures, from their estimated starting date *t**m* to phase-out *t**p* (cf. Table 1). Results of the model fitting are reported in Fig. 8. The SPQEIR model, with appropriate parameters for each country, is fitted to reported infection curves and, overall, model fitting have good agreement with data. This supports the model structure as very simple yet realistic enough to capture the main dynamical behaviour of the infection curves in multiple countries. In addition, it allows for each country to obtain multiple sets of parameters, providing a good fit and representing different strategies. Finally, it allows a comparison between different countries through the corresponding best fit parameters. For Ireland, although initial social distancing advises were issued on 13th March (cf. Table 1), fitting the complete curve was only possible when considering the lock-down date (28th March) as the major driver of the suppression. ![Fig. 8](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.04.22.20075804/F8.medium.gif) [Fig. 8](http://medrxiv.org/content/early/2020/11/24/2020.04.22.20075804/F8) Fig. 8 Results of model fitting. Infection curves for the considered countries (dotted) are fitted with the SPQEIR model with appropriate parameters (red curves). We also show a comparison with the fitted curve obtained from the “basic” SEIR model with only social distancing (turquoise curves). Parameter values are reported for each country, as well as the corresponding ![Graphic][32] (for the gray area, following Eq. 1) and ![Graphic][33]. The period of measures enforcement, from *t**m* to *t**p*, is highlighted by the grey region. Time progresses from the estimated day of first infection *t* (cf. Table 3). Population fraction refers to country-specific populations (cf. Table 1). After phase-out, we prolong the fitted curve (parameter values unchanged) to compare observed data with what could have been if measures had not been lifted (dashed lines). From the data, we can observe a resurgence of cases that points to possible “second outbreaks” (particularly in Israel). The model matches very well those regions like Austria and Switzerland that avoided saturation of the healthcare system and thus reported reliable data (cf. Fig. 4c). Ireland reported intermittent data, while Lombardy is not perfectly represented, probably because of data reporting and larger heterogeneity in its spacial patterns. Finally, the reduced *χ*2 metric (Eq. 3) usually attains minimum values for the complete SPQEIR model rather than the simple social distancing one. Hence, the SPQEIR is confirmed to be informative, on top of being fully interpretable and linked to recognised social policy categories. #### 3.3.2 Cross-country interventions assessment Fitting a number of countries with the same model containing the same epidemiological parameters allows to perform a quantitative and consistent comparison on the efficacy of their interventions. In Fig. 8, parameter values providing the best fit of model to data are reported, together with the simulation results (mean values) calculated by *lmfit* algorithm. Different synergies yield similar values for ![Graphic][34], but the curve is different in its evolution as already observed in the previous sections. As expected from the model analysis above, the lower ![Graphic][35] is (below 1), the faster the suppression of the epidemic. In addition, different parameter combinations generate differing curves, which might well explain differences in reported total cases and deaths between various countries. Comparing Austria, Denmark and Lombardy we can see that contact tracing and monitoring might play a role in speeding up the curve decay, despite the fact that population-wide interventions played a major role. In general, combined isolation and tracing strategies would reduce transmission in addition to social distancing or self-isolation alone. In particular, social distancing alone is effective only if very stringent, as suggested by “only social distancing” fitting. In general, a strong, rapid lock-down seems the best option, as also suggested by the conceptual analysis. However, intervening with additional synergies is a viable option to suppress the epidemic faster and with lower lock-down values. Finally, we observe the value of timely interventions: we see that intervening earlier with respect to the date of first infection helps reducing the daily curves by almost a factor of 10. For instance, we can compare Denmark and Lombardy in Fig. 8: the first one got a peak corresponding to about 0.08% of the whole population, while the second region registered a number of active cases of about 0.5% of the whole population. This translates in more than 3800 infectious on the Danish peak, and on more than 37000 on Lombardy’s. ## 4 Discussion The model is fitted until phase-out dates, when measures are progressively lifted and therefore the model assumptions do not hold anymore. In Fig. 8 we extrapolate the model, with same parameter values, after phase-out (dashed lines), to compare observed data to the most optimistic scenario, where measures would not have been lifted. We observe that, up to July 8th, the infection curves mostly maintained an inertial decreasing trend: despite some fluctuations that make them generally higher than the best scenario, they kept on following a downward trend similar to that of the model. We speculate that this phenomenon is linked to changed behaviors, face masks [39] and improved sanitising practices that maintained social distancing values, as well as contact tracing practices issued by many countries along with the phase-out. However, some countries (Israel in particular, but also Austria) showed a worrisome upward trend, possibly associate to a second outbreak. As this is not a low probability event, we stress the usefulness of our analysis to prepare for future developments in pandemic progression. It has been asked whether the peak of infections was reached because of herd immunity or because of interventions [40]. An added value of this study is to confirm that the peak of infection, for the considered countries, was not reached because of herd immunity. On the contrary, it is the effect of a number of suppression measures that reduced the number of cases artificially. This should warn about the high numbers of people that are still susceptible. We acknowledge the limitations of our analysis. Due to its structure and the use of ordinary differential equations, the model only accounts for average trends. However, it cannot reproduce fluctuations in the data, being them intrinsic in the epidemic, or from testing and reporting protocols that might differ among countries. In addition, the constant nature of parameters used in this analysis allows good agreement between model and data when countries implemented rapid and strong measures point-wise in time, with little follow-ups. Further statistical studies, with time varying parameters, could obtain more precise values. In the same way, transferring models from country to country requires fulfilling the same assumptions on model structure and basic hypothesis. This is shown by the different fitting performances, that suggest that a transfer is not always possible. In general, this study is not intended to make a ranking of country responses, nor to suggest that different strategies could have led to better outcomes. Contrariwise, it should be used as a methodological step towards quantitatively inquiring the effect of different intervention categories. It examines possible abstract scenarios and compares quantitative, model-based outputs, but it is not intended to fully represent specific countries nor to reproduce the epidemic complexity within societies. In fact, the model does not provide fine-grained quantification of specific interventions, e.g. how effective masks are in protecting people, how much proximity tracing apps increase *P**find*, how changes in behavior are associated with epidemic decline [41] and so on. We acknowledge that the new compartments cannot perfectly match policy measures, but are a reasonable approximation. Some real measures might also affect multiple parameters at once, e.g. safety devices and lock-down could impact both *µ* and *ρ*. Comparing results of this macro-scale model with those of complex, micro-scale ones [3] could inform researchers and policy makers about the epidemic dynamics and effective synergies to hamper it. Any conclusion should be carefully interpreted by experts, and the feasibility of tested scenarios should be discussed before reaching consensus. ## 5 Conclusion We have developed a minimal model to link intervention categories against epidemic spread to epidemiological model compartments. This allows quantitative assessment of non-pharmaceutical suppression strategies on top of social distancing, for a number of countries. Strategies have different effects on epidemic evolution in terms of curve flattening and eradication timing. As with previous studies [22, 42], we have observed the need to enforce containment measures (i.e., detect and isolate cases, identify and quarantine contacts and at risk neighborhoods) along with mitigation (i.e., slow down viral spread in the community with social distancing). By extending the classic SEIR model into the SPQEIR model, we distinguished the impact of different suppression programs in flattening the peak and anticipating the eradication of the epidemic. Depending on their strength and synergy, non-pharmaceutical interventions can hamper the disease from spreading in a population. First, we performed a complete sensitivity analysis of their effects, both alone and in synergy scenarios. Then, we moved from idealised representations to fitting realistic contexts, allowing preliminary mapping of intervention categories to abstract programs. We verified that the model is informative in interpolating the infection curves for a number of countries, and performed cross-country comparison. We could then obtain model-based outputs on the strength of interventions, for a number of countries that respected the model assumptions. This provides better, quantitative insights on the effect of suppression measures and their timing, and allows improved comparison. Overall, this work could contribute to quantitative assessments of epidemic suppression strategies. To tackle current epidemic waves, and against possible resurgence of contagion [43] (also cf. Fig 8), better understanding the effect of different non-pharmaceutical interventions could help planning mid- and long-term measures and to prepare preventive plans, until a vaccine is available. ## 6 Shinyapp A user-friendly online shinyapp to interactively simulate different scenarios with the SPQEIR model is available on: [https://jose-ameijeiras.shinyapps.io/SPQEIR\_model/](https://jose-ameijeiras.shinyapps.io/SPQEIR_model/). It allows to reproduce the present outputs and to perform sensitivity analysis. ## Data Availability No data are used. Simulation code is publicly available. [https://github.com/daniele-proverbio/Covid-19](https://github.com/daniele-proverbio/Covid-19) ## Ethics The application of anonymized data for the purpose of epidemic modelling has been endorsed by the accessed databases. ## Data accessibility Databases of social measures can be accessed at [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/phsm](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/phsm) [6]. ACAPS database is at [https://www.acaps.org/covid19-government-measures-dataset](https://www.acaps.org/covid19-government-measures-dataset) [17]. Worldwide epidemiological data collection from John Hopkins University is at [https://github.com/CSSEGISandData/COVID-19](https://github.com/CSSEGISandData/COVID-19) [18]. Lombardy data were retrived from [https://github.com/pcm-dpc/COVID-19](https://github.com/pcm-dpc/COVID-19) [19]. Google mobility data [21] were accessed through [https://ourworldindata.org/covid-mobility-trends](https://ourworldindata.org/covid-mobility-trends) [20]. The code for analysis can be found at [https://github.com/daniele-proverbio/assessing\_strategies](https://github.com/daniele-proverbio/assessing_strategies). ## Authors’ Contributions DP, AH designed the study. DP, FK, SM developed the model. DP, FK, JAA implemented the model. DP, FK, SM, AH, AA, LM, AS, JG, CL analyzed and interpreted the results. SM, AH, AS, JG, CL supervised and coordinated the project. DP, FK, SM, CL wrote the first draft. All authors contributed to the final draft. All authors gave their final approval for publication. ## Competing interest The authors declare no competing interests. ## Funding Fundings: DP and SM’s work is supported by the FNR PRIDE DTU CriTiCS, ref 10907093. FK’s work is supported by the Luxembourg National Research Fund PRIDE17/12244779/PARK-QC. A.H. work was partially supported by the Fondation Cancer Luxembourg. JG is partly supported by the 111 Project on Computational Intelligence and Intelligent Control, ref B18024. AA is supported by the Luxembourg National Research Fund (FNR) (Project code: 13684479). JAA is supported by the FWO research project G.0826.15N (Flemish Science Foundation), GOA/12/014 project (Research Fund KU Leuven), Project MTM2016-76969-P from the Spanish State Research Agency (AEI) co–funded by the European Regional Development Fund (ERDF) and the Competitive Reference Groups 2017–2020 (ED431C 2017/38) from the Xunta de Galicia through the ERDF. ## Acknowledgments The authors thank the Research Luxembourg - COVID-19 Taskforce for mutual collaborations. ## Footnotes * The title was slightly modified. In addition, we show the results of fitting the computational model to empirical data. * Received April 22, 2020. * Revision received November 24, 2020. * Accepted November 24, 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]. Yixiang Ng, Zongbin Li, Yi Xian Chua, Wei Liang Chaw, Zheng Zhao, Benjamin Er, Rachael Pung, Calvin J Chiew, David C Lye, Derrick Heng, et al. Evaluation of the effectiveness of surveillance and containment measures for the first 100 patients with COVID-19 in Singapore–January 2–February 29, 2020. 2020. 2. [2]. Giulia Giordano, Franco Blanchini, Raffaele Bruno, Patrizio Colaneri, Alessandro Di Filippo, Angela Di Matteo, and Marta Colaneri. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nature Medicine, pages 1–6, 2020. doi:[https://doi.org/10.1038/s41591-020-0883-7](https://doi.org/10.1038/s41591-020-0883-7). 3. [3]. Adam J Kucharski, Petra Klepac, Andrew Conlan, Stephen M Kissler, Maria Tang, Hannah Fry, Julia Gog, John Edmunds, CMMID COVID-19 Working Group, et al. Effectiveness of isolation, testing, contact tracing and physical distancing on reducing transmission of SARS-CoV-2 in different settings. medRxiv, 2020. doi: 10.1016/S1473-3099(20)30457-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30457-6&link_type=DOI) 4. [4]. Roy M Anderson and Robert M May. Population biology of infectious diseases: Part I. Nature, 280 (5721):361–367, 1979. doi:[https://doi.org/10.1038/280361a0](https://doi.org/10.1038/280361a0). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/280361a0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=460412&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.04.22.20075804.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979HE95900033&link_type=ISI) 5. [5]. Judith Legrand, Rebecca Freeman Grais, Pierre-Yves Boelle, Alain-Jacques Valleron, and Antoine Flahault. Understanding the dynamics of Ebola epidemicss. Epidemiology & Infection, 135(4):610–621, 2007. doi:[https://doi.org/10.1017/S0950268806007217](https://doi.org/10.1017/S0950268806007217). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268806007217&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16999875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.04.22.20075804.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246856500011&link_type=ISI) 6. [6].WHO. Tracking public health and social measures - a global dataset. [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/phsm](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/phsm), 2020. 7. [7]. Roy M Anderson, Hans Heesterbeek, Don Klinkenberg, and T Déirdre Hollingsworth. How will country-based mitigation measures influence the course of the COVID-19 epidemic? The Lancet, 395(10228):931–934, 2020. doi:[https://doi.org/10.1016/S0140-6736(20)30567-5](https://doi.org/10.1016/S0140-6736(20)30567-5). 8. [8]. William Ogilvy Kermack and Anderson G McKendrick. A contribution to the mathematical theory of epidemics. Proceedings of the royal society of london. Series A, Containing papers of a mathematical and physical character, 115(772):700–721, 1927. doi:[https://doi.org/10.1098/rspa.1927.0118](https://doi.org/10.1098/rspa.1927.0118). 9. [9]. Ping Yan and Shengqiang Liu. SEIR epidemic model with delay. The ANZIAM Journal, 48(1):119– 134, 2006. doi:[https://doi.org/10.1017/S144618110000345X](https://doi.org/10.1017/S144618110000345X). 10. [10]. Julien Arino, Jonathan R Davis, David Hartley, Richard Jordan, Joy M Miller, and P Van Den Driessche. A multi-species epidemic model with spatial dynamics. Mathematical Medicine and Biology, 22 (2):129–142, 2005. doi:[https://doi.org/10.1093/imammb/dqi003](https://doi.org/10.1093/imammb/dqi003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/imammb/dqi003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15778332&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.04.22.20075804.atom) 11. [11]. Neil M Ferguson, Derek AT Cummings, Christophe Fraser, James C Cajka, Philip C Cooley, and Donald S Burke. Strategies for mitigating an influenza pandemic. Nature, 442(7101):448–452, 2006. doi:[https://doi.org/10.1038/nature04795](https://doi.org/10.1038/nature04795). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature04795&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16642006&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.04.22.20075804.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239278900042&link_type=ISI) 12. [12]. Side Syafruddin and MSM Noorani. Seir model for transmission of dengue fever in selangor malaysia. IJMPS, 9:380–389, 2012. doi: 10.1142/S2010194512005454. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1142/S2010194512005454&link_type=DOI) 13. [13]. Christophe Fraser, Christl A Donnelly, Simon Cauchemez, William P Hanage, Maria D Van Kerkhove, T Déirdre Hollingsworth, Jamie Griffin, Rebecca F Baggaley, Helen E Jenkins, Emily J Lyons, et al. Pandemic potential of a strain of influenza A (H1N1): early findings. Science, 324(5934):1557–1561, 2009. doi: 10.1126/science.1176062. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMjQvNTkzNC8xNTU3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTEvMjQvMjAyMC4wNC4yMi4yMDA3NTgwNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. [14]. Christian L Althaus. Estimating the reproduction number of Ebola virus (EBOV) during the 2014 outbreak in West Africa. PLoS Currents, 6, 2014. doi: 10.1371/currents.outbreaks.91afb5e0f279e7f29e7056095255b288. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/currents.outbreaks.91afb5e0f279e7f29e7056095255b288&link_type=DOI) 15. [15]. Thomas Hale, Anna Petherick, Toby Phillips, and Samuel Webster. Variation in government responses to COVID-19. Blavatnik school of government working paper, 31, 2020. 16. [16].World Health Organization and others. Coronavirus disease 2019 (covid-19): situation report, 88. 2020. 17. [17].ACAPS. Covid-19 government measures dataset. [https://www.acaps.org/covid19-government-measures-dataset](https://www.acaps.org/covid19-government-measures-dataset), 2020. 18. [18]. Ensheng Dong, Hongru Du, and Lauren Gardner. An interactive web-based dashboard to track COVID-19 in real time. The Lancet infectious diseases, 20(5):533–534, 2020. doi:[https://doi.org/10.1016/S1473-3099(20)30120-1](https://doi.org/10.1016/S1473-3099(20)30120-1). [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%2F11%2F24%2F2020.04.22.20075804.atom) 19. [19].Dipartimento della Protezione Civile - Emergenza Coronavirus. Dati covid-19 italia. [https://github.com/pcm-dpc/COVID-19/tree/master/dati-regioni](https://github.com/pcm-dpc/COVID-19/tree/master/dati-regioni), Accessed 14 Jun 2020. 20. [20]. Esteban Ortiz-Ospina Max Roser, Hannah Ritchie and Joe Hasell. Coronavirus pandemic (covid-19). Our World in Data, 2020.[https://ourworldindata.org/coronavirus](https://ourworldindata.org/coronavirus). 21. [21].Google. COVID-19 Community Mobility Reports. [https://www.google.com/covid19/mobility/](https://www.google.com/covid19/mobility/), Accessed 14 Jun 2020. 22. [22]. Neil M Ferguson, Daniel Laydon, Gemma Nedjati-Gilani, Natsuko Imai, Kylie Ainslie, Marc Baguelin, Sangeeta Bhatia, Adhiratha Boonyasiri, Zulma Cucunubá, Gina Cuomo-Dannenburg, et al. Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand. Imperial College, London. DOI: [https://doi.org/10.25561/77482](https://doi.org/10.25561/77482), 2020. xdoi: 10.25561/77482. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.25561/77482&link_type=DOI) 23. [23]. Zhilan Feng. Final and peak epidemic sizes for SEIR models with quarantine and isolation. Mathematical Biosciences & Engineering, 4(4):675, 2007. doi: 10.3934/mbe.2007.4.675. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3934/mbe.2007.4.675&link_type=DOI) 24. [24]. Ying Liu, Albert A Gayle, Annelies Wilder-Smith, and Joacim Rocklö v. The reproductive number of COVID-19 is higher compared to SARS coronavirus. Journal of Travel Medicine, 2020. doi: [https://doi.org/10.1093/jtm/taaa021](https://doi.org/10.1093/jtm/taaa021). 25. [25]. Helen J Wearing, Pejman Rohani, and Matt J Keeling. Appropriate models for the management of infectious diseases. PLoS medicine, 2(7), 2005. doi:[https://doi.org/10.1371/journal.pmed.0020174](https://doi.org/10.1371/journal.pmed.0020174). 26. [26]. Liangrong Peng, Wuyue Yang, Dongyan Zhang, Changjing Zhuge, and Liu Hong. Epidemic analysis of COVID-19 in China by dynamical modeling. arXiv preprint arXiv:2002.06563, 2020. 27. [27]. Michael Y Li, John R Graef, Liancheng Wang, and János Karsai. Global dynamics of a SEIR model with varying total population size. Mathematical biosciences, 160(2):191–213, 1999. doi:[https://doi.org/10.1016/S0025-5564(99)00030-9](https://doi.org/10.1016/S0025-5564(99)00030-9). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0025-5564(99)00030-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10472754&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.04.22.20075804.atom) 28. [28]. Jane M Heffernan, Robert J Smith, and Lindi M Wahl. Perspectives on the basic reproductive ratio. Journal of the Royal Society Interface, 2(4):281–293, 2005. doi:[https://doi.org/10.1098/rsif.2005.0042](https://doi.org/10.1098/rsif.2005.0042). 29. [29].A Kucharski, T Russell, C Diamond, and Y Liu. Analysis and projections of transmission dynamics of nCoV in Wuhan. CMMID repository, 2, 2020. 30. [30]. Joseph T Wu, Kathy Leung, and Gabriel M Leung. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. The Lancet, 395(10225):689–697, 2020. doi:[https://doi.org/10.1016/S0140-6736(20)30260-9](https://doi.org/10.1016/S0140-6736(20)30260-9). 31. [31]. JR Taylor. An introduction to error analysis-university science books. Mill Valley, California, 1997. 32. [32]. Savi Maharaj and Adam Kleczkowski. Controlling epidemic spread by social distancing: Do it well or not at all. BMC Public Health, 12(1):679, 2012. doi:[https://doi.org/10.1186/1471-2458-12-679](https://doi.org/10.1186/1471-2458-12-679). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2458-12-679&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22905965&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.04.22.20075804.atom) 33. [33]. Xiuli Liu, Geoffrey JD Hewings, Minghui Qin, Xin Xiang, Shan Zheng, Xuefeng Li, and Shouyang Wang. Modelling the situation of COVID-19 and effects of different containment strategies in China with dynamic differential equations and parameters estimation. Available at SSRN 3551359, 2020. doi: [http://dx.doi.org/10.2139/ssrn.3551359](http://dx.doi.org/10.2139/ssrn.3551359). 34. [34]. Victor M Corman, Olfert Landt, Marco Kaiser, Richard Molenkamp, Adam Meijer, Daniel Kw Chu, Tobias Bleicker, Sebastian Brünink, Julia Schneider, Marie Luisa Schmidt, et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Eurosurveillance, 25(3):2000045, 2020. 35. [35]. Rudolf Hanel and Stefan Thurner. Boosting test-efficiency by pooled testing strategies for SARS-CoV-2. arXiv preprint arXiv:2003.09944, 2020. 36. [36].Deutsche Gesellschaft fü r Epidemiologie. Stellungnahme der Deutschen Gesellschaft fü r Epidemi-ologie (DGEpi) zur Verbreitung des neuen Coronavirus (SARS-CoV-2). [www.dgepi.de](http://www.dgepi.de), 2020. 37. [37]. Yang Liu, Rosalind M Eggo, and Adam J Kucharski. Secondary attack rate and superspreading events for SARS-CoV-2. The Lancet, 395(10227):e47, 2020. doi:[https://doi.org/10.1016/S0140-6736(20)30462-1](https://doi.org/10.1016/S0140-6736(20)30462-1). 38. [38]. Marco Vinceti, Tommaso Filippini, Kenneth J Rothman, Fabrizio Ferrari, Alessia Goffi, Giuseppe Maffeis, and Nicola Orsini. Lockdown timing and efficacy in controlling covid-19 using mobile phone tracking. EClinicalMedicine, page 100457, 2020. 39. [39]. Richard OJH Stutt, Renata Retkute, Michael Bradley, Christopher A Gilligan, and John Colvin. A modelling framework to assess the likely effectiveness of facemasks in combination with ‘lock-down’in managing the COVID-19 pandemic. Proceedings of the Royal Society A, 476(2238):20200376, 2020. doi:[https://doi.org/10.1098/rspa.2020.0376](https://doi.org/10.1098/rspa.2020.0376). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.2020.0376&link_type=DOI) 40. [40]. Lucy C Okell, Robert Verity, Oliver J Watson, Swapnil Mishra, Patrick Walker, Charlie Whittaker, Aris Katzourakis, Christl A Donnelly, Steven Riley, Azra C Ghani, et al. Have deaths from COVID-19 in Europe plateaued due to herd immunity? The Lancet, 2020. doi: 10.1016/S0140-6736(20)31357-X. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)31357-X&link_type=DOI) 41. [41]. Benjamin J Cowling, Sheikh Taslim Ali, Tiffany WY Ng, Tim K Tsang, Julian CM Li, Min Whui Fong, Qiuyan Liao, Mike YW Kwan, So Lun Lee, Susan S Chiu, et al. Impact assessment of non-pharmaceutical interventions against coronavirus disease 2019 and influenza in hong kong: an observational study. The Lancet Public Health, 2020. doi:[https://doi.org/10.1016/S2468-2667(2030090-6](https://doi.org/10.1016/S2468-2667(2030090-6)). 42. [42]. Corey M Peak, Rebecca Kahn, Yonatan H Grad, Lauren M Childs, Ruoran Li, Marc Lipsitch, and Caroline O Buckee. Individual quarantine versus active monitoring of contacts for the mitigation of covid-19: a modelling study. The Lancet Infectious Diseases, 2020. doi:[https://doi.org/10.1016/S1473-3099(20)30361-3](https://doi.org/10.1016/S1473-3099(20)30361-3). 43. [43]. Stephen M Kissler, Christine Tedijanto, Edward Goldstein, Yonatan H Grad, and Marc Lipsitch. Projecting the transmission dynamics of sars-cov-2 through the postpandemic period. Science, 368 (6493):860–868, 2020. doi: 10.1126/science.abb5793. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5My84NjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMS8yNC8yMDIwLjA0LjIyLjIwMDc1ODA0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/inline-graphic-5.gif [6]: /embed/graphic-2.gif [7]: /embed/inline-graphic-6.gif [8]: /embed/graphic-3.gif [9]: /embed/graphic-5.gif [10]: /embed/inline-graphic-7.gif [11]: /embed/graphic-6.gif [12]: /embed/inline-graphic-8.gif [13]: /embed/inline-graphic-9.gif [14]: /embed/inline-graphic-10.gif [15]: /embed/inline-graphic-11.gif [16]: F2/embed/inline-graphic-12.gif [17]: /embed/inline-graphic-13.gif [18]: /embed/inline-graphic-14.gif [19]: /embed/inline-graphic-15.gif [20]: /embed/graphic-11.gif [21]: /embed/inline-graphic-16.gif [22]: /embed/inline-graphic-17.gif [23]: /embed/inline-graphic-18.gif [24]: /embed/inline-graphic-19.gif [25]: /embed/inline-graphic-20.gif [26]: /embed/inline-graphic-21.gif [27]: /embed/inline-graphic-22.gif [28]: /embed/inline-graphic-23.gif [29]: /embed/inline-graphic-24.gif [30]: F7/embed/inline-graphic-25.gif [31]: /embed/inline-graphic-26.gif [32]: F8/embed/inline-graphic-27.gif [33]: F8/embed/inline-graphic-28.gif [34]: /embed/inline-graphic-29.gif [35]: /embed/inline-graphic-30.gif