Suppression of Groups Intermingling as Appealing Option For Flattening and Delaying the Epidemiologic Curve While Allowing Economic and Social Life at Bearable Level During COVID-19 Pandemic ============================================================================================================================================================================================== * Ioan Bâldea ## Abstract In this work, we simulate the COVID-19 pandemic dynamics in a population modeled as a network of groups wherein infection can propagate both via intra-group and via inter-group interactions. Our results emphasize the importance of diminishing the inter-group infections in the effort of substantial flattening and delaying of the epi(demiologic) curve with concomitant mitigation of disastrous economy and social consequences. To exemplify with a limiting case, splitting a population into m (say, 5 or 10) noninteracting groups while keeping intra-group interaction unchanged yields a stretched epidemiologic curve having the maximum number of daily infections reduced and postponed in time by the same factor m (5 or 10). More generally, our study suggests a practical approach to fight against SARS-CoV-2 virus spread based on population splitting into groups and minimizing intermingling between them. This strategy can be pursued by large-scale infrastructure reorganization of activity at different levels in big logistic units (e.g., large productive networks, factories, enterprises, warehouses, schools, (seasonal) harvest work). Importantly, unlike total lockdwon strategy, the proposed approach prevents economic ruin and keeps social life at a more bearable level than distancing everyone from anyone. ## 1 Introduction Neither natural immunization nor vaccination or pharmacologic intervention can currently help in fighting the COVID-19 global pandemic.1 Most frequently, to prevent disastrous sanitary consequences, governments across the world responded by imposing, e.g., social (=physical) distancing, wearing face masks, lockdown regulations, and rigid sanitary moves aiming at reducing the infection rate (“flattening the epidemiologic curve”). Still, politics cannot push strict restrictions indefinitely, and “how much is too much?” is a question of time which unavoidably arises sooner or later. Fighting COVID-19 should not ruin economy. 2 This is certainly what a draconian lockdown across the world over months would do. Parenthetically, catastrophic economic consequences inherently make healthcare system itself also collapsing. ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/10/2020.05.06.20093310/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093310/F1) Figure 1: Logistic growth is characterized by a cumulative number of infected cases *n*(*t*) following an S-shaped (sigmoid) curve (panel a) exhibiting an acceleration stage (close to a J-shaped exponential growth), which switches to a deceleration stage beyond the half-time *T*50 and attains *p*(= 95, 99, see equation (12)) percent of the plateau value *N*. Saturation occurs because the effective infection rate ![Graphic][1] is time dependent and gradually decreases to zero (panel b). The epi(demiologic) curve *ṅ*(*t*) = *dn*/*dt* (daily number of new cases) exhibits a peak located at *t* = *T*50 whose shape is controlled by the infection rate *κ* (panel c). Time on the *x*-axis is expressed in units of the characteristic time *Tc* = *κ*−1. In this vein, mathematical modeling may make a notable contribution in providing politics with reasonable suggestions to slowing down epidemic propagation and reducing medical burden while mitigating economy and social crisis. A series of mathematical COVID-19 simulations have recently appeared.3–10 Most of them are based on deterministic continuoustime epidemiologic models, which consider age-independent epidemiologic classes of, e.g., susceptible (S), exposed (E), infected (I), and recovered (R) individuals,11–16 whose numbers S(t), E(t), I(t), R(t) evolve in time (t) according to a system of (deterministic) ordinary differential equations. Open access sources already available17–19 enable one to easily perform various numerical simulations by means of such models. Unfortunately the various SIR-inspired flavors need (too many) input parameters difficult to validate,20 and this would rather mask than enlighten the main idea which the present work aims at conveying. Therefore, to better emphasize this idea, instead of a SIR-based approach (which would pose no special problem), in this paper we prefer to adopt the simpler logistic growth framework. The logistic model is particularly appealing in view of its simplicity and versatility demonstrated in approaching a broad variety of real systems with very different nature,21–28 including population dynamics of epidemic states.29–32 Prior to this study, results based on the logistic model were presented for COVID-19 time evolution in China and USA.33 Fighting against the spread of SARS-CoV-2 virus while allowing economic and social activity to continue to a reasonable extent represents a major challenge for the present era. Extended lockdown does not represent an acceptable response to this challenge. From this perspective, we believe that the results reported below obtained by extending the conventional logistic model may provide useful suggestions on how to sidestep the ongoing difficulty of living under pandemic conditions. While the implementation of the presently proposed strategy via population splitting into smaller groups and reducing intermingling certainly requires considerable effort and fantasy in infrastructure reorganization, it offers the perspective of flattening and delaying the epidemiologic curve by obviating wrecking of economy and maintaining social life to a level more bearable than total lockdown. It can also be pursued as a complementary approach to the massive COVID-19 testing proposed to fostering economy recovery.34 ## 2 Methods The results reported below were obtained by means of the logistic model21–28 extended (see equation (16)) to allow treatment of infection propagation in a network consisting of groups in interaction. To fix the ideas and to make the paper self-contained, in Section 3.1 we will first review the main aspects related to the logistic model applied to an isolated group using a terminology adapted to the specific subject under consideration. The extension of the logistic model to groups in interaction will be presented in Section 3.2. ## 3 Results and Discussion ### 3.1 Logistic Growth in an Isolated Group Uninhibited infected population n growths in time t according to the Malthus law35 ![Formula][2] The intrinsic population-independent rate K entering equation (1) is expressed in terms of the probability *β* of infection per encounter with an infected individual multiplied by the number ![Graphic][3] of encounters per unit time (day) ![Formula][4] This yields an unlimited exponential time growth (*n* = *n* (*t* = 0)) ![Formula][5] depicted by the dark green J-shaped curve of Figure 1a. In a real situation, the exponential growth will gradually slow down and eventually level off. Infections become more and more unlikely because, in a given environment, the increase in the number of infected diminishes the number of individuals that can be infected. Rephrasing, the effective growth rate decreases with increasing population density:![Graphic][6]. By assuming a linear decrease of *f* with *n*, one arrives at the logistic model21–26 ![Formula][7] ![Formula][8] Plotted as a function of time (Figure 1a), equation (5) yields an exponential J-shaped curve only at early times which switches to an S-shaped (sigmoid) curve as the population increases and saturates to the maximum (plateau) value N, which defines the so-called carrying capacity of a given environment. The parenthesis entering the right hand side of equation (5), which acts as Darwin’s “struggle for existence” and suppresses the exponential growth, is similar to the Pauli blocking factor extensively discussed in electron transport theory.36–40 The continuous time representation underlying the differential equation (5) allows to express the cumulative number of cases *n*(*t*) in closed analytical form ![Formula][9] ![Formula][10] By using the half-time *T*50 ![Formula][11] which defines the crossover point where the population attains *p* = 50 percent of its maximum (*N*), the above results can be recast as follows ![Formula][12] ![Formula][13] At *t* = *T*50, there is a substantial infection slowing with respect to the exponential growth. There, the instantaneous infection rate ![Graphic][14] is reduced by 50% as compared to that of the uninhibited growth *κ* (equation (4) and Figure 1b). Saturation occurs within a few characteristic times *Tc* = *κ*−1 beyond the half-time *T*50 (Figure 1a). The moment (“day”) *Tp* when the number *Np* of infected amounts to *p* percent (e.g., *p* = 95 or 99, cf. Figure 1) of the maximum value *N* ![Formula][15] can easily be deduced from equation (9) ![Formula][16] The quantity *κ* is important because it quantifies the daily new infected cases expressed by the time derivative *dn/dt* (equation (13) and the red curve in Figure 1c) ![Formula][17] The height *H* and full width at half maximum FWHM of this so-called epidemiologic curve *dn*/*dt* = *f* (*t*) are expressed by ![Formula][18] ![Formula][19] Large values of *κ*, amounting to a sharp and high peak in Figure 1c, may cause healthcare systems collapses. Diminishing the value of *κ* is of paramount practical importance for a twofold reason: 1. It renders the peak broader and smaller (Figure 1c). Keeping the daily number of infected cases at a manageable level (“flattening of the curve”) is essential for not overwhelming the healthcare system beyond its capacity to treat the sick. 2. Small values of κ yield large values of *T*50 (cf. equation (10)). This means postponement of infection explosion and hence gaining time for a better sanitary and logistic preparation to tackle an upcoming problem: preventing shortage of intensive care unit beds and gaining time for securing and/or producing critical emergency equipment (e.g., masks, ventilators, artificial lungs, personal protective equipment, extracorporeal membrane oxygenation machines, ventilators or other devices) needed to reduce mortality rate. ### 3.2 Modeling Infection in Interacting Groups The above considerations referred to a closed population in which members are neither added nor lost from the group. Neither “imported” nor “exported” infections were included. Let us now focus on a network consisting of groups of individuals {*nj*} = {*n*1; *n*2,…, *nm*} wherein infections can proliferate both by infections within the same group (intra-group infection rates *κj* ≡ *κjj*) and because individuals of one group j can infect or can be infected by individuals of other groups *p* ≠ *j* (inter-group/intermingling infection rates *κpj* and *κjp*, respectively). Figure 2 schematically depicts the case of two groups. ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/10/2020.05.06.20093310/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093310/F2) Figure 2: Schematic representation of a network consisting of two groups 1 and 2 wherein epidemic can spread (panel b) both through intra-group infections (infection rates *κ*1 and *κ*2) and through inter-group infections (infection rates *κ*12 and *κ*21). Limiting cases wherein the network is split into two groups that are completely separated among themselves *(κ*12 = *κ*21 = 0, panel a) or perfectly intermingled (*κ*12 = *κ*12 = *κ*21 = *κ*, panel c). Generalizing the idea underlying equation (5), we will consider below the following extended logistic model ![Formula][20] Although solving equation (16) in general poses no special numerical problem (some examples are presented in Figure 4 and Figure 5 of Section 3.4), to better emphasize the strategy we aim at conveying, let us first focus on the case of identical groups ![Formula][21] ![Formula][22] Equation (18) yields *j*-independent populations ![Formula][23] *nT* (*t*) being the total time dependent population. Upon term-by-term addition (![Graphic][24]) of equation (16) we immediately get ![Formula][25] and hence ![Formula][26] ![Formula][27] ![Formula][28] ![Formula][29] ![Formula][30] By comparing the above formulas with equations (6)-(10) valid for a single group one can conclude that the quantity ![Formula][31] plays the role of a total infection rate. Importantly, both the half-time *T*50 and the maximum number of daily infections *H* deduced from equations (23) and (24) ![Formula][32] ![Formula][33] are controlled by κ*T*. ### 3.3 Analysis of Two Limiting Cases of Practical Importance The results of Section 3.2 allow us to compare how infection propagates in a network (“larger group”) of individuals split into several (*m*) smaller (sub)groups (chosen identical for simplicity) which do not interact with each other against the case of a fictitious splitting, wherein (from the point of view of infection) interactions between members of a given group and between members belonging to different groups are identical (perfect intermingling). Noteworthy, whether the groups are completely separated of each other (label *s*) or perfectly intermingled (label *i*), the total initial population *n**T* = *mn* is taken to be the same in both cases (cf. equation (17)) ![Formula][34] The former case, corresponding to separated groups (label *i*), is characterized by a vanishing inter-group infection rate (*κ*′ ≡ 0). Applied to this case, equations (26), (27) and (28) yield ![Formula][35] At the opposite extreme of perfectly intermingled groups, inter-group interactions are as strong as intra-group interactions (*κ*′ = *κ*). Based on equations (26), (27) and (28) we then get ![Formula][36] The above results show that the epidemiologic curve can be substantially flattened if a larger group is split into several smaller groups separated from each other. By starting from the same number of infections (equation (29)), splitting into groups separated from one another yields a reduction of the infection rate and of the maximum daily cases by a factor *m* and an increase of full width at half maximum by the same factor ![Formula][37] ![Formula][38] ![Formula][39] Equally pleasantly, splitting leads in addition to a time postponement of the infection peak by the same factor *m* ![Formula][40] The results presented in Figure 3 depict these findings for the particular case of a larger group split into two smaller groups (*m* = 2). They also schematically visualize how and why group splitting can relieve the healthcare system. To avoid misunderstandings, one should note that the value *m* = 2 in Figure 3 was chosen just for more clarity. Splitting into more than two noninteracting groups (i.e., making *m* as large as possible) is highly desirable. ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/10/2020.05.06.20093310/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093310/F3) Figure 3: Results demonstrating how preventing intermingling between several groups flattens the epidemic curve with the concomitant time delay of the maximum of daily infections. (a) Cumulative number of cases *n*1,2(*t*) and (b) epidemiologic curve *ṅ*1,2(*t*) = *dn*1,2*(t)* for a group consisting of two subgroups that are either completely separated of each other (green lines) or perfectly intermingled (red curves). These results also schematically depict how flattening and delaying the epidemiologic curve by a factor *m* = 2 by splitting into *m* = 2 groups can prevent overloading the healthcare system capacity. Time on the *x*-axis is expressed in units of the characteristic time *Tc* = κ−1. Before ending this part, we want to emphasize that, however important, flattening and delaying the epidemiologic curve by a factor m achieved by group splitting (cf. Figure 3 and equations (32), (33), and (34)) is not the whole issue. Extremely importantly, the presently proposed group splitting approach does not assume any intra-group (like social distancing and wearing masks) restrictions: economy and social life within individual groups separated of each other can continue. ### 3.4 Additional Results As anticipated in Section 3.2, the logistic model extended as expressed by equation (16) can be used to quantify the impact of mutual infections in networks of interacting groups more general than the particular situations examined in Section 3.3. To briefly illustrate this fact, two examples are presented in Figures 4 and 5. ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/10/2020.05.06.20093310/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093310/F4) Figure 4: (a) Cumulative number of cases *n*1,2 (*t*) and (b) epidemiologic curves *ṅ*1,2(*t*) = *dn*1,2(*t*)/*dt* for two completely separated and perfectly intermingled (red and blue curves, respectively) groups whose populations differ by a factor of four. Notice that, in spite of the equal infection rates, infection of the smaller group 2 (dashed lines) is stronger affected by the inter-group interaction than in the larger group 1 (solid lines). Time on the x-axis is expressed in units of the characteristic time *Tc* = *k*−1. ![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/10/2020.05.06.20093310/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093310/F5) Figure 5: Results showing the impact of the inter-group interaction *K* on the cumulative number of cases *n*1,2(*t*) and epidemiologic curves *ṅ*1,2(*t*) = *dn*1,2(*t*)/*dt* in case of two groups merely differing by the initial number of infected cases: *n*01/*N* = 0.001 and *n*02/*N* = 0.1. They show how infections in the initially less infected group 1 are rapidly triggered by infections in the initially more infected group 2. Time on the *x*-axis is expressed in units of the characteristic time *Tc* = *κ*−1. Figure 4 depicts the case two groups whose populations differ by a factor of four. As visible there, in spite of the equal infection rates, infection of the smaller group 2 (dashed lines) is stronger enhanced by the inter-group interaction than in the larger group 1 (solid lines). Figure 5 presents the case two groups merely differing from each other by the different numbers of initially infected individuals *(n*01*/N* = 0.001 versus *n*02/*N* = 0.1). Comparison of the various panels of Figure 5 reveals that inter-group infection yields an infection rapidly “exported” from the initially more infected group to that which was initially less infected. Intermingling (Figure 5g) quickly wipes out any difference between a initially weaker (or non)infected group and an initially strongly infected group. ## 4 Conclusion In closing, the results presented above indicate that splitting of a large group into smaller groups and reducing intermingling appears to be an appealing strategy for substantially reducing the spread of the SARS-CoV-2 virus while still allowing social life to a more bearable level than distancing everyone from anyone and economy go on, albeit slower. One can expect that (A) it is easier to impose and maintain longer-term regulations on splitting a given population into (say, 5 – 10) groups weakly interacting among themselves than (B) enforcing severe containment measures to all individuals in order to diminish the probability *β* of infection per encounter and the number ![Graphic][41] of daily encounters yielding a reduction by the same factor (5–10) of the infection rate k (cf. equation (2)). The big difference is that with option (A) social life and economy within individual groups go on without intra-group restrictions, while option (B) means attaining the same epidemiologic curve with both social life and economy paralyzed. Devising and implementing an adequate restructuring of large logistic units (large productive networks, factories, enterprises, warehouses, etc) allowing, if/when necessary, society to rapidly switch back and forth between separated groups and intermingled groups can certainly be challenging but may be a long-term strategic goal worth to be pursued when faced with these COVID-19 pandemic times or other similar difficulties that cannot be ruled out in the future. ## Data Availability The manuscript includes all information necessary to generate and reproduce the data reported in the present work. * Received May 6, 2020. * Revision received May 6, 2020. * Accepted May 10, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. 1.Cucinotta, D.; Vanelli, M. Acta Bio Medica Atenei Parmensis 2020, 91, 157–160. 2. 2.1. Baldwin, R.; 2. di Mauro, B. W. Mitigating the COVID Economic Crisis: Act Fast and Do Whatever It Takes; Baldwin, R.; di Mauro, B. W., Eds.; A CEPR Press VoxEU.org eBook, 2020. 3. 3.Chen, T.; Rui, J.; Wang, Q.; Zhao, Z.; Cui, J.-A.; Yin, L. *bioRxiv* 2020, DOI 10.1101/2020.01.19.911669. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4xOS45MTE2Njl2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzEwLzIwMjAuMDUuMDYuMjAwOTMzMTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. 4.Su, L.; Hong, N.; He, J.; Ma, Y.; Wang, X.; Xu, Q.; Jiang, H.; Han, L.; Chang, F.; Zhu, W.; Long, Y. *medRxiv* 2020, DOI 10.1101/2020.03.06.20032177. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4wNi4yMDAzMjE3N3YzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTAvMjAyMC4wNS4wNi4yMDA5MzMxMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. 5.Q, L.; X, G.; et al., W. P. N Engl J Med. 2020, 382, 1199–1207. 6. 6.Roosa, K.; Lee, Y.; Luo, R. e. a. J. Clin. Med. 2020, 9, 596. 7. 7.Shen, M.; Peng, Z.; Xiao, Y.; Zhang, L. *bioRxiv* 2020, DOI 10.1101/2020.01.23.916726. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4yMy45MTY3MjZ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzEwLzIwMjAuMDUuMDYuMjAwOTMzMTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 8. 8.Hui, D.; I Azhar, E.; Madani, T.; Ntoumi, F.; Kock, R.; Dar, O.; Ippolito, G.; Mchugh, T.; Memish, Z.; Drosten, C.; Zumla, A.; Petersen, E. International Journal of Infectious Diseases 2020, 91, 264–266. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.01.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31953166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F10%2F2020.05.06.20093310.atom) 9. 9.Yang, P.; Qi, J.; Zhang, S.; Wang, X.; Bi, G.; Yang, Y.; Sheng, B. *medRxiv* 2020, DOI 10.1101/2020.04.01.20043794. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4wMS4yMDA0Mzc5NHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTAvMjAyMC4wNS4wNi4yMDA5MzMxMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 10. 10.Aggarwal, N. *medRxiv* 2020, doi:10.1101/2020.04.17.20069245. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4xNy4yMDA2OTI0NXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTAvMjAyMC4wNS4wNi4yMDA5MzMxMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. 11.Kermack, W. O.; McKendrick, A. G. Proc. Roy. Soc. 1927, 115A, 700–721, reprinted in Bull. Math. Biol. 53, 3355 (1991). [https://doi.org/10.1007/BF02464423](https://doi.org/10.1007/BF02464423). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1927.0118&link_type=DOI) 12. 12.Kermack, W. O.; McKendrick, A. G. Proc. Roy. Soc. 1932, 138A, 55–83, reprinted in Bull. Math. Biol. 53, 5787 (1991). [https://doi.org/10.1007/BF02464424](https://doi.org/10.1007/BF02464424). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1932.0171&link_type=DOI) 13. 13.Kermack, W. O.; McKendrick, A. G. Proc. Roy. Soc. 1933, 141A, 94–122, reprinted in Bull. Math. Biol. 53, 89118 (1991). [https://doi.org/10.1007/BF02464425](https://doi.org/10.1007/BF02464425). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1933.0106&link_type=DOI) 14. 14.Bailey, N. T. et al. The Mathematical Theory of Infectious Diseases and Its Applications; Charles Griffin & Company Ltd, 5a Crendon Street, High Wycombe, Bucks HP13 6LE., 1975. 15. 15.Hethcote, H. W. Frontiers in Mathematical Biology, Berlin, Heidelberg, 1994; pp 504–515. 16. 16.Hethcote, H. W. SIAM Review 2000, 42, 599–653. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1137/s0036144500371907&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00016567&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F10%2F2020.05.06.20093310.atom) 17. 17.Morris, S. E. shinySIR: Interactive Plotting for Mathematical Models of Infectious Disease Spread; 2020, R package version 0.1.1. 18. 18.Epidemiological Models for Influenza and COVID-19, [https://community.wolfram.com/groups/-/m/t/1896178](https://community.wolfram.com/groups/-/m/t/1896178). 19. 19.Goh, G. Epidemic Calculator, 2020, [http://gabgoh.github.io/COVID/index.html](http://gabgoh.github.io/COVID/index.html). 20. 20.Fong, S. J.; Li, G.; Dey, N.; Crespo, R. G.; Herrera-Viedma, E. International Journal of Interactive Multimedia and Artificial Intelligence 2020, 6, 132–140. 21. 21.Quetelet, L. A. J. Du Système Social et des Lois qui le Régissent; Guillaumin, 1848. 22. 22.Ostwald, W. Journal für Praktische Chemie 1883, 27, 1–39. 23. 23.McKendrick, A. G.; Pai, M. K. Proceedings of the Royal Society of Edinburgh 1912, 31, 649–653. 24. 24.Lloyd, P. Population Studies 1967, 21, 99–108. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11630390&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F10%2F2020.05.06.20093310.atom) 25. 25.Cramer, J. Studies in History and Philosophy of Science Part C: Studies in History and Philosophy of Biological and Biomedical Sciences 2004, 35, 613—626. 26. 26.Vandermeer, J. Nature Education Knowledge 2010, 3, 15. 27. 27.Bâldea, I. Phys. Chem. Chem. Phys. 2017, 19, 30842—30851. 28. 28.Baldea, I. Phys. Chem. Chem. Phys. 2018, 20, 23492–23499. 29. 29.Mansfield, E.; Hensley, C. Journal of the Royal Statistical Society: Series B (Methodological) 1960, 22, 332–337. 30. 30.Waggoner, P. E.; Aylor, D. E. Annual Review of Phytopathology 2000, 38, 71–94, PMID: 11701837. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.phyto.38.1.71&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11701837&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F10%2F2020.05.06.20093310.atom) 31. 31.Koopman, J. Annual Review of Public Health 2004, 25, 303–326, PMID: 15015922. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.publhealth.25.102802.124353&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15015922&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F10%2F2020.05.06.20093310.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221130600015&link_type=ISI) 32. 32.Bangert, M.; Molyneux, D. H.; Lindsay, S. W.; Fitzpatrick, C.; Engels, D. Infect Dis. Poverty 2017, 6, 73. 33. 33.Hermanowicz, S. W. *medRxiv* 2020, DOI 10.1101/2020.03.31.20049486. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4zMS4yMDA0OTQ4NnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTAvMjAyMC4wNS4wNi4yMDA5MzMxMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 34. 34.Baldwin, R. COVID-19 testing for testing times: Fostering economic recovery and preparing for the second wave, 2020, [https://voxeu.org/article/testing-testing-times](https://voxeu.org/article/testing-testing-times). 35. 35.Malthus, T. R. An Essay on the Principle of Population, as It Affects the Future Improvement of Society with Remarks on the Speculations of Mr. Godwin, M. Condorcet, and Other Writers; Printed for J. Johnson, in St. Paul’s Church-Yard, London, 1798. 36. 36.Datta, S. Phys. Rev. B 1992, 45, 1347–1362. 37. 37.Sols, F. Ann. Phys. (NY) 1992, 214, 386—438. 38. 38.Bönig, L.; Schönhammer, K. Phys. Rev. B 1993, 47, 9203–9207. 39. 39.Datta, S. Electronic Transport in Mesoscopic Systems; Cambridge Univ. Press: Cambridge, 1997. 40. 40.Wagner, M. Phys. Rev. Lett. 2000, 85, 174–177. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10991187&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F10%2F2020.05.06.20093310.atom) [1]: F1/embed/inline-graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/inline-graphic-2.gif [4]: /embed/graphic-3.gif [5]: /embed/graphic-4.gif [6]: /embed/inline-graphic-3.gif [7]: /embed/graphic-5.gif [8]: /embed/graphic-6.gif [9]: /embed/graphic-7.gif [10]: /embed/graphic-8.gif [11]: /embed/graphic-9.gif [12]: /embed/graphic-10.gif [13]: /embed/graphic-11.gif [14]: /embed/inline-graphic-4.gif [15]: /embed/graphic-12.gif [16]: /embed/graphic-13.gif [17]: /embed/graphic-14.gif [18]: /embed/graphic-15.gif [19]: /embed/graphic-16.gif [20]: /embed/graphic-18.gif [21]: /embed/graphic-19.gif [22]: /embed/graphic-20.gif [23]: /embed/graphic-21.gif [24]: /embed/inline-graphic-5.gif [25]: /embed/graphic-22.gif [26]: /embed/graphic-23.gif [27]: /embed/graphic-24.gif [28]: /embed/graphic-25.gif [29]: /embed/graphic-26.gif [30]: /embed/graphic-27.gif [31]: /embed/graphic-28.gif [32]: /embed/graphic-29.gif [33]: /embed/graphic-30.gif [34]: /embed/graphic-31.gif [35]: /embed/graphic-32.gif [36]: /embed/graphic-33.gif [37]: /embed/graphic-34.gif [38]: /embed/graphic-35.gif [39]: /embed/graphic-36.gif [40]: /embed/graphic-37.gif [41]: /embed/inline-graphic-6.gif