Reproduction ratio and growth rates: measures for an unfolding pandemic ======================================================================= * Maíra Aguiar * Joseba Bidaurrazaga Van-Dierdonck * Nico Stollenwerk ## Abstract The initial exponential growth rate of an epidemic is an important measure that follows directly from data at hand, commonly used to infer the basic reproduction number. As the growth rates ***λ*(*t***) of tested positive COVID-19 cases have crossed the threshold in many countries, with negative numbers as surrogate for disease transmission deceleration, lockdowns lifting are linked to the behavior of the momentary reproduction numbers ***r*(*t***), often called *R*. Important to note that this concept alone can be easily misinterpreted as it is bound to many internal assumptions of the underlying model and significantly affected by the assumed recovery period. Here we present our experience, as part of the Basque Country Modeling Task Force (BMTF), in monitoring the development of the COVID-19 epidemic, by considering not only the behaviour of *r*(*t*) estimated for the new tested positive cases - significantly affected by the increased testing capacities, but also the the momentary growth rates for hospitalizations, ICU admissions, deceased and recovered cases, in assisting the Basque Health Managers and the Basque Government during the lockdown lifting measures. Two different data sets, collected and then refined during the COVID-19 responses, are used as an exercise to estimate the momentary growth rates and reproduction numbers over time in the Basque Country, and the implications of using those concepts to make decisions about easing lockdown and relaxing social distancing measures are discussed. These results are potentially helpful for task forces around the globe which are now struggling to provide real scientific advice for health managers and governments while the lockdown measures are relaxed. ## 1 Introduction As the COVID-19 pandemic is unfolding, research on mathematical modeling became imperative and very influential, not only in understanding the epidemiology of COVID-19 but also in helping the national health systems to cope with the high demands of hospitalizations, for example, providing projections and predictions based on the available data. Used as a public health guiding tools to evaluate the impact of intervention measures, governments have already taken important decisions based on modeling results [1, 2, 3]. While for diseases which are long established, such as measles [4], e.g., modeling results are easier to be interpreted, given the availability of long term and well established data collections, and public health interventions, counting with an effective vaccine, are able to be be implemented in time to avoid large outbreaks. For COVID-19 the situation is completely different. We are now dealing with a single disease outbreak in a pandemic scenario and modeling projections, for instance, need to be adjusted for the new scientific information and new data that are generated every day under unprecedentedly fast changes of circumstances. Governments around the globe are relying on quick measurements updates and long established concepts such as the reproduction numbers which are bound to to many internal assumptions of the underlying model, smoothing and approximations [5] and significantly affected by the assumed recovery period. The initial exponential growth rate of an epidemic is an important measure that follows directly from data at hand, commonly used to infer the basic reproduction number *R*, which is the number of secondary cases generated from a primary infected case during its infectiousness before recovering in a completely susceptible population [6]. Both concepts can be extended to larger compartmental models and into the phase when effects of the control measures become visible and parameters slowly change, leading to the so called momentary momentary growth rates *λ*(*t*), and momentary reproduction ratios *r*(*t*). In the beginning of COVID-19 epidemics, the process of collecting data were often not yet well organized or pre-organized in the way that we could immediately use to feed models and extract accurate measurements for the momentary growth rates and the momentary reproduction numbers. To mitigate and suppress COVID-19 transmission, draconian intervention measures were rapidly implemented, crippling our economies as lockdowns were implemented. As research to develop an effective vaccine is ongoing, epidemiologists and public health workers are the frontline of this pandemic, focusing on the well known public health surveillance strategies of testing, isolation and contact tracing of infected COVID-19 individuals. Up to date, more than 4 million cases were confirmed with about 300 thousand deaths, and these numbers are still increasing [7]. After several weeks of social distancing restrictions, lockdowns start now to be lifted and modeling task forces around the globe are struggling to apply the concept of *r*(*t*), often called *R*, to decide whether social distancing relaxation decisions are taken in the right period of time, i.e, when the outbreak is assumed to be controlled, with negative growth rates and a momentary reproduction numbers below 1. Although the absolute value of *r*(*t*) can vary, countries we should rather look at the threshold behavior as it is independent, when we use the growth rates primarily, of those modeling uncertainties and clearly indicates if the outbreak is under control or not when estimations are below or above 1. Complementary measures of growth rates for different variables such as hospitalization, intensive care unites (ICU) admissions and deceased, where data is also collected, should be evaluated when political decisions are taken. In this report we present the growth rates and reproduction numbers for the COVID-19 epidemic in the Basque Country, an autonomous community in northern Spain with 2.2 million inhabitants. For the reproduction number calculation we use a refined stochastic SHARUCD-type model - an extension of the well known simple SIR model that is frequently used to model different disease outbreaks [8, 9, 10], developed within a multidisciplinary task force (so-called Basque Modelling Task Force, BMTF) created to assist the Basque Health managers and the Basque Government during the COVID-19 responses. The model is calibrated using the empirical data provided by the Basque Health Department and the Basque Health Service (Osakidetza), continually collected with specific inclusion and exclusion criteria. Able to describe well the incidences of disease for different variables of tested positive individuals, this framework is now used to monitor disease transmission, including estimations of the momentary growth rates and reproduction numbers, while the country lockdown is gradually lifted [11]. Using two different available data sets for the Basque Country, collected from March 4 to May 9, 2020, the data was revised as variable definition for positive cases was changed in respect to the diagnostic test used, we present results obtained for the momentary growth rates and reproduction ratios during the ongoing COVID-19 epidemic in the Basque Country and discuss the implications of using those concepts during an unfolding pandemic. ## 2 Materials and Methods For the Basque Country we use the cumulative data for the following variables defined as: i) total tested positive patients (*Icum*) in yellow which are recorded in categories for ii) hospital admissions (*CH*), in red, iii) intensive care units admissions (*CU*) in purple, iv) recovered (*CR*) in green and v) deceased (*D*) in black. At the beginning of the outbreak, only patients with severe symptoms admitted to a hospital were tested using the PCR (polymerase chain reaction) method. As testing capacities increased, including also antibody tests used mainly as screening tool in nursing homes, less severe symptomatic cases started to be tested, contributing to enhance the number of confirmed positive cases in the population. This data collection (named “Data set A”) includes, for each category or variable, patients tested with both PCR and rapid antibody tests. “Data set A” has now being revised to include patiences, in all categories, who were tested positive with PCR method only (named “Data set B”). Using the data for the all positive cases, the momentary growth rates (*λ*) and the momentary reproduction numbers (*r*(*t*)) are calculated for both data sets, A and B, and results are compared. ### 2.1 The underlying mathematical model and empirical data We consider primarily SHARUCD model versions as stochastic processes in order to compare with the available data which are often noisy and to include population fluctuations, since at times we have relatively low numbers of infected in the various classes. The stochastic version can be formulated through the master equation in application to epidemiology in a generic form using densities of all variables *x*1:= *S/N, x*2*:= H/N, x*3*:= A/N, x*4*:= R/N, x*5*:= U/N, x*6*:= CH/N, x*7*:= CA/N, x*8*:= CU/N* and *x*9:= *D/N* and *x*10:= *CR/N* hence state vector ![Graphic][1], giving the dynamics for the probabilities ![Graphic][2] as ![Formula][3] with *n* =10 different transitions ![Graphic][4], as described by the mechanisms above, and small deviation from state ![Graphic][5] as ![Graphic][6]. For the basic SHARUCD model we have explicitly the following transitions ![Graphic][7] and its shifting vectors ![Graphic][8] given by ![Formula][9] With these ![Graphic][10] and ![Graphic][11] specified we also can express the mean field ODE system The deterministic version of the refined model is given by a differential equation system for all classes, including the recording classes of cumulative cases *CH*, *CA*, *CR* and *CU* by ![Formula][12] Model parameters and initial conditions are shown in Table 1, where *β* is the infection rate and *ϕ* is the ratio describing the asymptomatic/mild infections contribution to the force of infection. *γ* is the recovery rate, *μ* is the disease induced death rate and *ν* is the ratio of hospitalized going to the ICU. *η* is the proportion of susceptible being infected, develop sever symptoms and being hospitalized whereas 1 − *η* is the proportion of susceptible becoming infected and developing mild disease or asymptomatic. *ξ* is the ratio of detected, via testing, mild/asymptomatic infect individuals. *ϱ* is the import rate needed to describe the introductory phase of the epidemics and for the present study, we assume *ϱ* to be much smaller than the other additive terms of the force of infection, given the strong observational insecurities on the data collected at the beginning of the outbreak. View this table: [Table 1:](http://medrxiv.org/content/early/2020/05/21/2020.05.18.20105528/T1) Table 1: Model parameters and initial condition values. For completeness of the system and to be able to describe the initial introductory phase of the epidemic, an import term *ϱ* should be also included into the force of infection. For the present study, we assume *ϱ* to be much smaller than the other additive terms of the force of infection, given the strong observational insecurities on the data collected at the beginning of the outbreak, when *ϱ* would matter most. ### 2.2 Growth rate After an introductory phase, the epidemic entered into an exponential growth phase, which started in the Basque Country around the March 10, 2020 and due to the effects of the imposed control measures has left to a slower growth around March 27, 2020 [11]. This exponential growth phase is typical for any outbreak with disease spreading in a completely susceptible population, as observed already in the SIR-system, from the dynamics of the infected ![Graphic][13] when *S*(*t*) ≈ *N*, such that a linear differential equation ![Graphic][14] with an exponential growth factor A is obtained. This growth factor then can be measured again from disease data via ![Graphic][15] giving a straight line in a semi-logarithmic plot of the data. For larger compartmental models we obtain similarly an exponential growth factor. For the present SHARUCD model we have the active disease classes *H* and *A* with the dynamics given by ![Formula][16] now including disease induced transition to death via the mortality rate *μ* and transition to ICU admission with admission rate *ν*. For an epidemic in its initial phase, i.e. *S*(*t*) ≈ *N*, we now have constant matrices ![Graphic][17] for entries into the disease classes and ![Graphic][18] for exits from the disease classes, where we had infection rate *β* and recovery rate *γ* in the SIR case. With ![Graphic][19] we now have with *J* = *B − G* the dynamics ![Graphic][20] and its solution ![Formula][21] with matrix exponential including the eigenvalue matrix Λ and the transformation matrix *T* from the eigenvectors of matrix *J*. The eigenvalues of the matrix *J* are given by ![Formula][22] with the parameter dependent trace *tr* = (*η* + *ϕ*(1 − *η*))·*β* − (2*γ* + *μ* + *ν*) and determinant *det = γ*(*γ* + *μ* + *ν*) − ((*γ* + *μ* + *ν*)*ϕ*(1 − *η*) + *γη*)·*β* and the dominating growth factor is given by the largest eigenvalue *λ*1. After an initial introductory phase the exponential growth with *λ*1 dominates the dynamics of *H*(*t*) and *A*(*t*), and from there also all the other variables, because the remaining equations are all inhomogeneous linear differential equations with the inhomogeneities given by the solutions *H*(*t*) and *A*(*t*), and we have ![Formula][23] with constants *KH* and *KA* depending on parameters and initial conditions (from ![Graphic][24] etc.). In the limiting case of a simple SIR-type model (with *ϕ* ≈ 1 and *μ*, *ν* ≪ *γ*) we obtain *tr* = (*η* + *ϕ*(1 − *η*))·*β* − (2 *γ* + *μ* + *ν*) ≈ *β* − 2*γ* and *det* = *γ*(*γ* + *μ* + *ν*) − ((*γ* + *μ* + *ν*)*ϕ*(1 − *η*) + *γη*)·*β* ≈ *γ*2 − *γβ* and hence *λ*1 ≈ *β* − *γ* and *λ*2 ≈ −*γ* The concept of the growth rate can be extended into the phase when effects of the control measures become visible and parameters slowly change, such that for short times the above analysis holds as for constant parameters. The momentary growth rates are analyzed below. ### 2.3 Reproduction ratio Another measure of the spreading of the disease in its initial phase is the basic reproduction number (*R*), the number of secondary cases *Is* from a primary case *Ip* during its infectiveness before recovering in a completely susceptible population. In its simplest version for SIR models a primary case, *Ip*(*t*) = 1, recovers via ![Graphic][25] hence ![Graphic][26]. The number of secondary cases from the primary case is given by ![Graphic][27] with *Is*(*t*) = 0, a simple inhomogeneous linear differential equation in case of a entirely susceptible population *S*(*t*) = *N*. The solution is ![Graphic][28] and gives the total number of secondary cases from a primary case as the long time limit as ![Graphic][29], hence the basic reproduction number is simply ![Graphic][30]. So we have the relation between ![Graphic][31] and the growth rate *λ* here as ![Graphic][32]. Generalized, the reproduction ratio is then given by ![Graphic][33] as the ratio of secondary cases produced by primary cases during their infectiousness. This concept can be also generalized for larger compartmental models, with the notions of matrices *B* and *G* as introduced above. For any primary cases *Hp* or *Ap* we have with ![Graphic][34] the decay dynamics ![Graphic][35] with solution ![Graphic][36], using again the matrix exponential. For secondary cases *Hs* and *As* we have the dynamics of ![Graphic][37] given by ![Graphic][38] with solution analogously to the SIR case as ![Graphic][39] with *F* = *BG−*1 the next generation matrix, since ![Graphic][40] or from generation ![Graphic][41] to generation ![Graphic][42] the discrete iteration ![Graphic][43]. For the present case we have the next generation matrix given as ![Formula][44] with its dominant eigenvalue for the basic SHARUCD-model ![Formula][45] and the other one being zero. In the limiting case of a simple SIR-type model (with *ϕ* ≈ 1 and *μ*, *ν* ≪ *γ*) we obtain again *r*1 = β/*γ* as can be easily seen. This concept of the reproduction ratio can be extended into the phase when effects of the control measures become visible and parameters slowly change. The momentary reproduction ratios (*r*) can be analyzed, as frequently done for the COVID-19 epidemics, but often called “basic reproduction number”. While the momentary growth rate follows directly from the time continuous data at hand, the momentary reproduction ratio depends on the notion of a generation time *γ−*1. we use smoothing of the differences of logarithmic positive cases with *τ* = 7 *d* and Δ*t* = 5 *d* To obtain the momentary growth rates from data directly we use ![Graphic][46] at first applied to the cumulative tested positive cases *Icum*(*t*) obtaining, via a smoothing window, the new cases after time *τ* as ![Formula][47] and hence, the growth rate ![Formula][48] ## 3 Results and discussion From the growth rate, the reproduction ratio is calculated with the recovery period *γ−*1 obtained from our underlying model and recent literature about SARS-CoV-2 interaction with human hosts [12, 13, 14, 15, 16, 17]. Assuming the recovery period to be of 10 days, we use 7 days smoothing of the differences of logarithmic positive cases to include all possible fluctuations during the data collection process such as “weekend effects”, for example, when we often observe a consistent low number of cases reported that are then adjusted shortly after. For this first exercise, both data sets show negative growth rates from April 1st, 2020, confirming a decrease in disease transmission. Nevertheless, when looking at the long term results for the “Data set A” (see Fig. 2 a-b)), an increase of the growth rate over time is estimated, with values crossing the threshold and becoming positive from April 23 to May 1, 2020, whereas “Data set B” (see Fig. 2 c-d), measures were kept constantly negative, without any signal of increasing disease and hence, the growth rate transmission. The momentary reproduction numbers follow the same trends for both data sets respectively, depending on the data set used. The observed signal from “Data set A” would significantly impact decisions on lockdown lifting, as the national plan for lifting the restrictions imposed during the state of alarm, called “Plan for the Transition towards a new normality”, was announced on April 28, 2020 [18]. Taking place over 4 phases, with a gradual de-escalation to “a new normality”, the plan is dependent on the on-going progress of COVID-19 epidemic’s control across the different regions of Spain. However, results obtained by “Data set B” would, alternatively, support the already started lockdown lifting with its “Phase Zero” initiated on May 4, 2020. When assuming a short recovery period of *γ* = 4 days, see Fig. 4, similar results are observed between the different data sets, only with variation on the absolute values. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/21/2020.05.18.20105528/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/05/21/2020.05.18.20105528/F1) Figure 1: Ensemble of stochastic realizations of the SHARUCD-type model. a) Cumulative tested positive cases *Icum*(*t*). From April 6, 2020, we note an increase of reported positive cases as the testing capacities were increasing. In b) cumulative hospitalized cases *CH*(*t*), c) cumulative ICU admission *CU*(*t*), d) cumulative deceased cases *D*(*t*). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/21/2020.05.18.20105528/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/05/21/2020.05.18.20105528/F2) Figure 2: Momentary growth rates estimation from the data on positive tested infected cases in a) PCR + rapid tests and c) PCR alone. The momentary reproduction ratios from the same data respectively are shown in c) and d), for *γ*−1 = 10. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/21/2020.05.18.20105528/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/05/21/2020.05.18.20105528/F3) Figure 3: Momentary growth rates estimation from the data on positive tested infected cases in a) PCR + rapid tests and c) PCR alone. The momentary reproduction ratios from the same data respectively are shown in c) and d) for *γ*−1 = 4. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/21/2020.05.18.20105528/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/05/21/2020.05.18.20105528/F4) Figure 4: Using data on PCR + rapid tests we plot the momentary growth rates estimation from the data on positive tested infected cases (yellow), hospitalizations (red) and ICU admission (purple) are plotted in a) and recovered (green) and deceased cases (black) in b). Using data on PCR tests only we plot the momentary growth rates estimation from the data on positive tested infected cases (yellow), hospitalizations (red) and ICU admission (purple) are plotted in c) and recovered (green) and deceased cases (black) in d). The momentary growth rates for the various variables are also calculated to verify and support the interpretation of the estimated *r*(*t*) threshold behaviour since for any assumed recovery period *γ*−1, results obtained for the various variables are the same, changing only when considering the different data sets. Figure 4 shows the behavior of three variables that are synchronized in the Basque Country, *Icum*, *CH*, and *CU*. They also cross the threshold to a negative growth rate on April 1st, 2020, confirming the observed r(t) trend obtained by looking at data on *Icum* alone. Recovered and deceased cases, shown in Fig. 4 b), follow 1-2 weeks later, due to the delay between onset of symptoms, hospitalization and eventually death, reaching negative growth rate on April 7 and April 11, 2020, respectively. Besides the observed deviation of *I*, for the “Data set A”, the other variables are kept below the threshold, constantly negative until May 9, 2020, supporting the political decision of starting lifting the lockdown measures rather sooner than later in time. So which measure should be considered to guide political decisions? Here, the answer is simple. When the available data is consistently collected and defined, the momentary growth rates for different variables, *I*, *H*, *U*, *R* and *D*, measured directly from the data at hand, should also be considered as complementary investigation. As the concept of *R* used alone can be easily misinterpreted, specially now when testing capacity is increasing and consequently the number of new notified cases, the BMTF now monitors the development of the COVID-19 epidemic in the Basque Country by considering not only the behaviour of the momentary growth rates *λ*(*t*) and momentary reproduction numbers *r*(*t*) for the positive cases *Icum*(*t*), but also the *λ*(*t*) for hospitalizations (*CH*), ICU admissions (*CU*), deceased (*D*) and recovered cases (*CR*), assisting the Basque Health Managers and the Basque Government with results that are obtained by the model framework, based on available data and evidence as scientific advise. Without interfering in any political decision, we now use “Data set B”, with a clearer definition of tested positive cases *Icum* and all other variables that follow, *CH*, *CU*, *CR* and *D*, and recovery period of *γ* = 4, shown in Figure 3 d) and Figures 4 c-d). At the moment, the reproduction ratio *r* is estimated to be below the threshold behavior of *r* =1, but still close to 1, meaning that although the number of new cases reported in the Basque Country are decelerating, the outbreak is still in its linear phase and careful monitoring of the development of the dynamics of the new cases from all categories, based on new information and data, to support the upcoming political decisions that will change the current life situation of millions of people is required. Using the available data for the Basque Country, a small community with short path for data collection and validation, we developed a modeling framework able to predict the course of the epidemic, from introduction to control measure response, potentially helpful for task forces around the globe which are now struggling to provide real scientific advice for health managers and governments while the lockdown measures are relaxed.. ## Data Availability Epidemiological data used in this study are provided by the Basque Health Department and the Basque Health Service (Osakidetza), continually collected with specific inclusion and exclusion criteria. [https://www.euskadi.eus/boletin-de-datos-sobre-la-evolucion-del-coronavirus/web01-a2korona/es/](https://www.euskadi.eus/boletin-de-datos-sobre-la-evolucion-del-coronavirus/web01-a2korona/es/) ## Acknowledgments Maíra Aguiar has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skodowska-Curie grant agreement No 792494. Data were provided by Basque Health Department. We thank the huge efforts of the whole COVID-19 BMTF, specially to Eduardo Millán for collecting and preparing the data sets used in this study. We thank Adolfo Morais Ezquerro, Vice Minister of Universities and Research of the Basque Goverment for the fruitful discussions. * Received May 18, 2020. * Revision received May 18, 2020. * Accepted May 21, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. [1]. R. M. Anderson, H. Heesterbeek, D. Klinkenberg, T. D. Hollingsworth, How will country-based mitigation measures influence the course of the COVID-19 epidemic? The Lancet. 395, 931–934 (2020). 2. [2]. S. Flaxman et al., “Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand” (Imperial College COVID-19 Response Team published on 16 March, 2020. Retrieved from [https://www.imperia/.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-9-impact-of-npis-on-covid-19/](https://www.imperia/.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-9-impact-of-npis-on-covid-19/) 3. [3]. 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 14 Apr 2020:eabb5793 4. [4].World Health Organization. Measles in the European region” (WHO Situation report 2, published in August 2019). Retrieved from [http://www.euro.who.int/\_\_data/assets/pdf\_fi/e/0010/410779/Measles-Sitrep-2-Aug-2019\_.pdf?ua=1](http://www.euro.who.int/\__data/assets/pdf_fi/e/0010/410779/Measles-Sitrep-2-Aug-2019_.pdf?ua=1) 5. [5]. P.L. Delamater, E.J. Street, T.F. Leslie, Y. Yang, K.H. Jacobsen, Complexity of the Basic Reproduction Number (*R*). Emerg. Infect. Dis. 25(1), 1–4 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2501.171901&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30560777&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.18.20105528.atom) 6. [6]. J. A. P. Heesterbeek, K. Dietz, The concept of *R* in epidemic theory. Stat. Neerl. 50, 89–110 (1996). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1467-9574.1996.tb01482.x&link_type=DOI) 7. [7].World Health Organization. “Situation report - 113 Coronavirus disease 2019 (COVID-19)” (WHO report published on 12 May, 2020. Retrieved from [https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200512-covid-19-sitrep-113.pdf?sfvrsn=feac3b6dJ2](https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200512-covid-19-sitrep-113.pdf?sfvrsn=feac3b6dJ2)) 8. [8]. M. Aguiar, W. B. Kooi, F. Rocha, P. Ghaffari, N. Stollenwerk, How much complexity is needed to describe the fluctuations observed in dengue hemorrhagic fever incidence data? Ecological Complexity, 16, 31–40 (2013). 9. [9]. M. Aguiar, S. Ballesteros, B. W. Kooi, N. Stollenwerk, The role of seasonality and import in a minimalistic multi-strain dengue model capturing differences between primary and secondary infections: complex dynamics and its implications for data analysis. Journal of Theoretical Biology, 289, 181–196 (2011). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21907213&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.18.20105528.atom) 10. [10]. N. Stollenwerk, S. van Noort, J. Martins, M. Aguiar, F. Hilker, A. Pinto, G. Gomes, A spatially stochastic epidemic model with partial immunization shows in mean field approximation the reinfection threshold. Journal of Biological Dynamics, 4, 634–649 (2010). 11. [11]. M. Aguiar et al., [https://medrxiv.org/cgi/content/short/2020.05.10.20086504v1](https://medrxiv.org/cgi/content/short/2020.05.10.20086504v1) (2020). 12. [12]. Stephen A. Lauer et al. (2020). The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Annals of Internal Medicine doi: [https://doi.org/10.7326/M20-0504](https://doi.org/10.7326/M20-0504) 13. [13]. Yang Liu et al. (2020). Viral dynamics in mild and severe cases of COVID-19. The Lancet Infectious Diseases. doi:[https://doi.org/10.1016/S1473-3099(20)30232-2](https://doi.org/10.1016/S1473-3099(20)30232-2) 14. [14]. J. A. Tetro. (2020). Is COVID-19 receiving ADE from other coronaviruses? Microbes Infect. 22, 72–73. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.micinf.2020.02.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32092539&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.18.20105528.atom) 15. [15]. Y. Fu, Y. Cheng, Y. Wu. (2020). Understanding SARS-CoV-2-mediated inflammatory responses: From mechanisms to potential therapeutic tools. Virol. Sin. doi:10.1007/s12250-020-00207-4pmid:32125642 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12250-020-00207-4pmid:32125642&link_type=DOI) 16. [16]. Francois-Xavier Lescure et al. (2020). Clinical and virological data of the first cases of COVID-19 in Europe: a case series. The Lancet Infectious Diseases doi:[https://doi.org/10.1016/S1473-3099(20)30200-0](https://doi.org/10.1016/S1473-3099(20)30200-0) 17. [17]. Ruiyun Li et al. (2020). Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). (2020). Science, 368(6490), 489-493 DOI: 10.1126/science.abb3221 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80ODkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8yMS8yMDIwLjA1LjE4LjIwMTA1NTI4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 18. [18]. Gobierno de España, April 28th, 2020. Plan para la Transicin hacia una nueva normalidad. Retrieved from [https://www.mscbs.gob.es/profesionales/saludPublica/ccayes/alertasActual/nCov-China/planDesescalada.htm](https://www.mscbs.gob.es/profesionales/saludPublica/ccayes/alertasActual/nCov-China/planDesescalada.htm) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/graphic-1.gif [4]: /embed/inline-graphic-3.gif [5]: /embed/inline-graphic-4.gif [6]: /embed/inline-graphic-5.gif [7]: /embed/inline-graphic-6.gif [8]: /embed/inline-graphic-7.gif [9]: /embed/graphic-2.gif [10]: /embed/inline-graphic-8.gif [11]: /embed/inline-graphic-9.gif [12]: /embed/graphic-3.gif [13]: /embed/inline-graphic-10.gif [14]: /embed/inline-graphic-11.gif [15]: /embed/inline-graphic-12.gif [16]: /embed/graphic-5.gif [17]: /embed/inline-graphic-13.gif [18]: /embed/inline-graphic-14.gif [19]: /embed/inline-graphic-15.gif [20]: /embed/inline-graphic-16.gif [21]: /embed/graphic-6.gif [22]: /embed/graphic-7.gif [23]: /embed/graphic-8.gif [24]: /embed/inline-graphic-17.gif [25]: /embed/inline-graphic-18.gif [26]: /embed/inline-graphic-19.gif [27]: /embed/inline-graphic-20.gif [28]: /embed/inline-graphic-21.gif [29]: /embed/inline-graphic-22.gif [30]: /embed/inline-graphic-23.gif [31]: /embed/inline-graphic-24.gif [32]: /embed/inline-graphic-25.gif [33]: /embed/inline-graphic-26.gif [34]: /embed/inline-graphic-27.gif [35]: /embed/inline-graphic-28.gif [36]: /embed/inline-graphic-29.gif [37]: /embed/inline-graphic-30.gif [38]: /embed/inline-graphic-31.gif [39]: /embed/inline-graphic-32.gif [40]: /embed/inline-graphic-33.gif [41]: /embed/inline-graphic-34.gif [42]: /embed/inline-graphic-35.gif [43]: /embed/inline-graphic-36.gif [44]: /embed/graphic-9.gif [45]: /embed/graphic-10.gif [46]: /embed/inline-graphic-37.gif [47]: /embed/graphic-11.gif [48]: /embed/graphic-12.gif