Retrospective Analysis of Interventions to Epidemics using Dynamic Simulation of Population Behavior ==================================================================================================== * Jenna Osborn * Shayna Berman * Sara Bender-Bier * Gavin D’Souza * Matthew Myers ## Abstract Retrospective analyses of interventions to epidemics, in which the effectiveness of strategies implemented are compared to hypothetical alternatives, are valuable for performing the cost-benefit calculations necessary to optimize infection countermeasures. SIR (susceptible-infected-removed) models are useful in this regard but are limited by the challenge of deciding how and when to update the numerous parameters as the epidemic changes in response to population behaviors. We present a method that uses a “dynamic spread function” to systematically capture the continuous variation in the population behavior, and the gradual change in infection dynamics, resulting from interventions. No parameter updates are made by the user. We use the tool to quantify the reduction in infection rate realizable from the population of New York City adopting different facemask strategies during COVID-19. Assuming a baseline facemask of 67% filtration efficiency, calculations show that increasing the efficiency to 75% could reduce the roughly 5000 new infections per day occurring at the peak of the epidemic to 3000. Mitigation strategies that may not be varied as part of the retrospective analysis, such as social distancing, are automatically captured as part of the calibration of the dynamic spread function. Key words * COVID-19 * SIR model * Infection-spread model * personal protective equipment * facemask ## 1. Introduction Retrospective analyses of strategies used to contain epidemics such as COVID-19 are valuable for countering successive waves of the infection, selecting countermeasures for future epidemics, and educating the population regarding the efficacy of implementing certain behavioral modifications. In particular, public-health agencies responsible for recommending types of personal protective equipment (PPE) to stockpile in anticipation of a future pandemic can benefit greatly from the cost-benefit information yielded by retrospective analyses. Mathematical models, including those of the SIR type, can be helpful in providing a quantitative framework for the analyses. SIR models have been applied during the COVID-19 pandemic1,2,3,4, primarily in a predictive capacity. To our knowledge, no studies have attempted to re-create an actual scenario (e.g. CNYC) with different interventions, though in Cooper et al.4 the influence of different mask types on the reproduction number was computed using a generic infection scenario. A formidable challenge in applying SIR models is prescribing the values of the numerous parameters, and updating them to simulate evolving infection dynamics, as population behaviors change in response to interventions. Typically, the behaviors change in a continuous manner, as in the gradual adoption of face masks by an affected population. Such gradual changes are difficult to capture by periodically adjusting the parameters characterizing the behaviors. The challenge is further accentuated by the high sensitivity of the predictions to some of the parameter values2. Often, parameter choices are based upon best guesses, or closeness of fit (sometimes visual) of computed profiles with published curves4. In this paper we introduce a modification of traditional SIR models that incorporates a “dynamic spread” function that captures changes in population behavior in a continuous manner. There is no need to adjust parameters as interventions are implemented during the course of the infection. The spread function satisfies a differential equation with variable coefficients obtained from published infection curves for the epidemic under study. The computed spread function reproduces the infection dynamics resulting from a given intervention strategy. Subsequently, the spread function can be systematically modified to analyze the effect of alternate intervention strategies. We illustrate the process using the COVID-19 crisis in New York City (CNYC). The reduction in infection rate realizable in CNYC with an increased level of mask usage, and with deployment of masks offering higher levels of filtration, is estimated. ## 2. Methods We illustrate the technique using a 3-equation SIR model. It assumes that the infection dynamics are dominated by one transmission mode (e.g. airborne particulates), and that the parameter values are appropriate for all particle sizes contributing to that mode (though interpretation of the resulting equations as an average for a broad particle distribution is possible5). More complicated SIR models can be useful, particularly if it is desired to model the details of the infection dynamic, e.g. symptomatic and asymptomatic individuals1. Our intention is to use the simplest model that can capture the known infection dynamics in a general sense, with the hope that the model can be understood and used by non-experts such as policy makers. Additionally, as noted by Siegenfeld et al.6, simpler models can prove more useful than complex ones, in part because accurate data is often not available to inform complicated formulations. Finally, we expect that some of the technique we present can be extended to more complex models. ### 2.1 Overview of Strategy The dynamic-spread function is part of a systematic procedure for calibrating SIR models. The 5 steps in the procedure are listed below and implemented subsequently. 1. Use the rate of change (measured by the number of new infections per day), *dS/dt*, of the susceptible population, as the variable for calibrating the model, rather than *S*. The derivative profile, which we call *T(t)*, was felt to be the quantity known most accurately. Its determination does not require the number of recovered patients to be tracked. 2. Normalize variables and identify critical dimensionless parameters. Formulating the model in terms of dimensionless quantities reduces the number of property values that must be determined, which can be numerous in SIR models7. 3. Allow the dimensionless parameter *δ*, which is essentially the product of the infection transmission rate and the virus production rate, to vary with time, and account for its time dependence in the governing differential equation for *T(t)*. We denote *δ(t)* the “dynamic-spread” function, as it contains the elements that both vary with time and govern the rate of spread of the infection. 4. Derive the governing equation for *δ(t)*. Provide the required coefficient functions using published *T(t)* profiles for a given infection scenario. 5. Designate as the time origin for the dynamic analysis the point of the first intervention into the epidemic. For CNYC, we identify this as day 17 from the first reported infection, when shelter-in-place was instituted. Prior to that point, it is assumed that *δ* is constant in time, and a standard SIR model applies. The parameters for the standard SIR model can be estimated from the published growth rate and reproduction number. The resulting values serve as initial conditions for the dynamic analysis. ### 2.2 Development of Governing Equations The SIR equations under these conditions listed in sections 2 and 2.1 are as follows 7,5. The primes denote dimensional quantities and will be dropped following nondimensionalization. ![Formula][1] ![Formula][2] ![Formula][3] Here *S’* is the number of susceptible individuals in the total population *N, I’* the number of infected individuals, *D’* the total number of droplets, ![Graphic][4]the transmission rate, *μI* the infection recovery rate, *κ* the droplet production rate, and *ν*-1 the droplet removal rate. We next introduce the maximum number of newly reported infections (roughly 5000 per day for CNYC) *α*, and a time scale Δ, which we take to be the turn-around time (time required to reach dT’/dt’ = 0, about 37 days for CNYC.) The function *T’* is scaled by *α, I’* by *α*Δ, and *D’* by *κα*Δ*2* Additionally, we differentiate Eq. 1a (after nondimensionalization), allowing the spread parameter ![Formula][5] to vary with time. Also, for huge populations as in CNYC, we ignore terms proportional to *α*Δ*/N*, including *S’/N – 1*. Under these assumptions, the SIR equations become: ![Formula][6] ![Formula][7] ![Formula][8] where *γ = Δμ**I* is the dimensionless infection recovery rate and ![Graphic][9]is the dimensionless droplet removal rate. Inserting (2d) in (2b) and using (2a) yields the 2-equation system: ![Formula][10] ![Formula][11] As noted above, we take the time origin to be the time of first intervention. To determine the dynamic spread function, we reformulate Eq. (3a) as an equation for *δ(t)*, assuming *T(t)* and *I(t)* to be known. We then obtain the *T(t)* profile from the published number of new infections per day in the locale of interest (e.g. New York City). We label this published profile *T**p**(t)* and the resulting (from Eq. (3b)) infection profile *I**p**(t)*, and we insert them into the equation for *δ(t)*. The resulting equation for the dynamic spread function is: ![Formula][12] Because the governing equation for *δ(t)* is informed by the published *T**p**(t)* profile, solving Eq’s 3 using this dynamic spread function will reproduce (within numerical tolerances) the published *T**p**(t)* curve. The utility of *δ(t)* derives from modifying it to model alternative intervention strategies and solving Eqs (1) to determine the modified infection rate. Modifications to account for protective strategies were performed in the following manner. ### 2.3 Accounting for Protective Equipment We build upon a previously developed SIR model5,8 that systematically accounts for the presence of protective equipment. Differentiating Eq. (1d) with respect to time yields ![Formula][13] Apportioning a fraction ϵ*κ* (e.g. 1/2) of the change in *δ* to changes in droplet production, we set ![Formula][14] Since from (1d) ![Formula][15] ![Formula][16] which can be integrated to ![Formula][17] In Myers et al5, it was shown that the production rate in the presence of protective equipment can be written as ![Formula][18] Here FE is the filtration efficiency (e.g. the FE for an N95 respirator is 95%) of the mask for the dominant droplet size, and *fi* is the fraction of the infected population wearing the covering at any given time. To perform a retrospective analysis in which a barrier material of different capturing efficiency is investigated, the new FE value would be used in Eq. (10) which, with Eq. (9), would be used to create a new dynamic-spread function. The modified spread function would then be used (in Eqs (3)) to estimate the change in infection rate. To analyze scenarios where different fractions of the infected population deploy a given mask type, Eqs. (9) and (10) can be combined to give ![Formula][19] And ![Formula][20] Eq. (12) can be used to estimate the fraction of infected individuals deploying a mask of a given filtration efficiency for a baseline case (known *δ(t)*). The effect of different fractions of the infected population deploying the mask can be quantified by prescribing *fi(t)*, solving for *δ(t)*, and using this modified spread function in Equations (3). ### 2.4 Solution Technique The initial conditions for Equations (3) are obtained by simulating the dynamics of the infection prior to any intervention, days 1 – 17 for CNYC. In that case, the derivative of the spread function is zero and Equations (3a,b) revert to a traditional SIR model. Seeking solutions that have an exponential time dependence of the form *exp*(Mt) result in the algebraic equation ![Formula][21] The subscript “0” on *δ* implies that the value applies to the initial period of the infection, before intervention occurs. The other parameters do not vary during the course of the epidemic and are not subscripted. An exponentially growing solution will occur when ![Formula][22] The ratio *δ****/(γλ)* is the reproduction number *R*** 5 for the standard SIR model. The growth rate *M* can be obtained from infection rates published during the beginning of the infection. Estimates of the reproduction numbers *R*** for the early stages of epidemics are also published. In the simulations, a range of recovery times *μI* (dimensionless recovery times *γ*) ranging from 2 days to 10 days were considered. For any given value of *γ, λ* and *δ*** were obtained from Eqs (13) and (14) and used in the solution of the dynamic equations (3). The initial value for *T(t)* was obtained from the published profile *T**p**(t)* (published number of new infections per day9), evaluated at day 17. *I**p**(t)* was derived from *T**p**(t)* using Eq. (3b), rather than using a published infection profile, so that it was not necessary to ascertain how well recoveries were tabulated in the published infection curves. The *I**p**(t)* function evaluated at day 17 was used to provide the initial condition for I(t). Equations (3) were solved using a Runge-Kutta method (Matlab *ode45*, Mathworks Inc.). ## 3. Results We performed a retrospective analysis of CNYC during days 17 - 37. This interval was chosen because day 17 is the day of the first intervention (shelter in place), and day 37 is the time of maximum new infections per day, based upon a 7-day average9. Initial reproduction numbers between 2 and 6 were considered, along with recovery times between 2 days and 10 days. The fraction ϵ*κ* was 1/2. We analyzed scenarios where the infected population deployed different types of masks. For baseline, it was assumed that the FE was 67%. This value is representative of homemade masks10, though the filtration capability of homemade masks spans a wide range. Higher-efficiency masks with FE’s of 75%, 80%, and 90% were considered for the hypothetical scenarios. Uncertainty in the calculated results was obtained by performing simulations for different combinations of reproduction numbers and recovery rates and computing the mean and standard deviation for the ensemble of parameter combinations. Figure 1a shows the dynamic spread function as a function of time. A sharp decrease is seen initially, owing to the shelter-in-place restriction. For larger FE, a sharper decrease in the spread function is observed. A large decrease in new infections (Fig. 1b) accompanies a small reduction in spread function value. Increasing FE from 67% to 75%, for example, reduces the spread function value about 10% at day 37, and maximum number of new infections (at day 37) decreases by about 40%. The turn-around time is decreased from about 37 days to 32 days. ![Figure 1a.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/27/2020.11.09.20228684/F1.medium.gif) [Figure 1a.](http://medrxiv.org/content/early/2021/01/27/2020.11.09.20228684/F1) Figure 1a. Dynamic spread function for CNYC with infected population deploying masks with different filtration efficiencies. Shadowed regions denote values within a standard deviation of the mean, for an ensemble of simulations using different reproduction numbers and recovery rates. ![Figure 1b.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/27/2020.11.09.20228684/F2.medium.gif) [Figure 1b.](http://medrxiv.org/content/early/2021/01/27/2020.11.09.20228684/F2) Figure 1b. Number of new infections per day for CNYC, with infected population deploying masks with different filtration efficiencies. Shadowed regions denote values within a standard deviation of the mean, for an ensemble of simulations using different reproduction numbers and recovery rates. For the same increase of FE from 67% to 75%, the number of infected individuals (Fig. 1c) at day 37 is reduced by about 30%. The uncertainty is considerably larger for the spread function (Fig 1a) and the infected population (Fig. 1c) than the number of new infections per day (Figure 1b), as the spread function and infected population are much more strongly influenced by the recovery time than the number of new infections. The recovery time spanned a factor of 5 over all the simulations performed. ![Figure 1c.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/27/2020.11.09.20228684/F3.medium.gif) [Figure 1c.](http://medrxiv.org/content/early/2021/01/27/2020.11.09.20228684/F3) Figure 1c. Infected population as a function of time for CNYC, when infected population deploys masks of different filtration efficiencies. Shadowed regions denote values within a standard deviation of the mean, for an ensemble of simulations using different reproduction numbers and recovery rates. An additional set of simulations was performed in which the fraction of the CNYC population wearing a mask with a filtration efficiency of 0.67 was increased, by 2%, 5%, and 10%. The baseline fraction deploying the mask was determined from Eq. (12), with *δ(t)* (baseline curve in Fig. 1a) derived from Eq. (4). The fraction deploying the mask increased from 0% at day 17 to roughly 75% on day 37. Increasing the fraction by 10% (i.e. multiplying the *fi* value at each time by 1.1) reduced the number of new infections per day from about 5100 to 3600 (Fig. 2). The turn-around time is reduced from approximately 37 days to 32 days. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/27/2020.11.09.20228684/F4.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/01/27/2020.11.09.20228684/F4) Figure 2. Number of new infections per day for scenarios where various fractions of the infected population in CNYC deploy masks with a 67% filtration efficiency. ## 4. Discussion While days 17 to 37 were featured in our simulations, the dynamic-spread-function technique can be applied to any time interval where reliable numbers of new infections are available. The standard SIR model is used prior to the time when either the production rate or the transmission rate is altered by an intervention strategy. At that point the dynamic simulations commence, with the SIR results serving as initial conditions. The model is not intended to be a predictive tool, in the sense of forecasting the future course of an ongoing epidemic. The purpose of the model is to compare different intervention strategies for scenarios where the baseline infection profile (number of new infections per day) is provided. Also required are the initial reproduction number and an estimate of the recovery rate. Though the model is not a forecasting tool, it can be useful for designing future countermeasures, particularly if elements of the anticipated scenario are similar to those of the scenario used to compute the spread function *δ(t)*. These elements include, most importantly, population behaviors such as face mask adoption (affecting both *κ* and ![Graphic][23]in Eq. 5) and social distancing (affecting ![Graphic][24]), but also environmental factors such as the pathogen inactivation rate. We refer the reader to Stilianakis and Drossinos7 for the dependence of infection dynamics on the numerous properties of the pathogen, the population, and the environment. Here we emphasize that the spread function implicitly captures the influence of all these factors, even though the functional dependence of the parameters is not introduced. Only when considering an alternative scenario that varies one of the factors does the does the explicit parametric dependence matter. In this paper, the functional dependence of the production rate *κ* upon the mask filtration efficiency is introduced via Eq. (10) in order to evaluate different PPE strategies. Calibration of the model is based upon number of new infections, rather than the size of the infected population, to eliminate the uncertainty associated with how well recoveries are tracked in the calibration dataset. To compare with published data sets where recoveries are not accounted for, then the integral of *T(t)* (converted to dimensional form) up to the time point of interest would be the proper metric for representing the size of the infected population. Since the dynamic spread function is the product of the transmission rate times the droplet production rate, it quantifies the ability of the infection to spread. The ability of the infection to spread decreased rapidly from day 17 (Fig. 1a). Because adoption of masks in CNYC likely occurred on a continuous basis over the weeks succeeding day 17, the curves in Fig. 1a are smooth and monotonically decreasing. For longer periods of time, where mask usage may eventually decrease, the spread function need not be monotonic. The effects of different protective-equipment strategies in CNYC were investigated without having to update the SIR parameters during the epidemic. The continuous adoption of masks would be difficult to simulate by updating coefficients at various times in standard SIR models. With the dynamic-spread approach, the gradual adoption of masks is captured in a natural manner. As noted above, the dynamic-spread approach allowed social distancing to be captured without being specifically modeled. The pattern of social distancing in CNYC was retained for all the simulated scenarios involving different facemasks. This commonality is largely responsible for the similar shapes of the curves in Fig. 1a. The only decision made relative to social distancing was that the factors contained in the transmission rate ![Graphic][25], which includes social distancing5,7, were responsible for roughly half (ϵ*κ* = 1/2) of the reduction in the spread function (shown in Fig. 1a). Other fractions would result in different reductions in infection in Figs 1 and 2, with higher values of ϵ*κ* resulting in larger reductions in infection rate, and vice versa. The simulated scenarios addressed only changes in masks worn by the infected population. No change in protection for the susceptible population was assumed. The susceptible population deployed masks, but the type was not varied between scenarios. As shown in Myers et al.5, the effect of mask deployment by the susceptible population can be simulated by modifying the transmission term ![Graphic][26] in Eq. (3), using equations analogous to (9) and (10). For the conditions of the simulations, a slight increase in facemask efficiency resulted in a larger benefit than a commensurate increase in compliance. At day 37, for example, a fractional increase in compliance of 0.1 resulted in a reduction in new infections of about 1500 per day (Fig. 2), while a fractional increase in FE of 0.1 reduced the number of new infections by about 1800 (interpolating Fig. 1b). For a higher baseline FE than 67%, increasing the compliance rate would produce a larger decrease in new infections. This comparison between filter efficiency and population compliance illustrates the utility of the model for determining how resources devoted to countermeasures can be optimally spent. In this case, the model can help inform the choice between 1) producing and distributing barriers of higher FE, and 2) educating and incentivizing the population to deploy barriers more readily available. A noteworthy conclusion emerging from the simulations is that considerable benefit can be obtained from higher FE masks without requiring N95 levels of efficiency (Fig 1). It is important to emphasize that for the benefits to be realized, the filtration efficiencies for the barrier material must be attainable for the particle-size range of the dominant transmission mode for the given scenario. One way of assuring this is for the barrier to provide the given FE across the spectrum of particle sizes. Otherwise, knowledge of the material filtration efficiency for the intended application (e.g. reducing airborne particulates generated by coughing or sneezing by infected persons indoors) is required in order to generate useful estimates. The complex issues of dominant transmission mode for COVID-19, and the filtration efficiency of different masks designs for the different modes, will be addressed in future applications of the model. ## Data Availability No human, animal, or in-vitro data is included in the manuscript, which describes a risk-assessment model. Algorithms used to derive the mathematical results have been referenced. ## Data Availability Data used to inform the model was obtained from the Johns Hopkins Coronavirus Resource Center (2020), [https://coronavirus.jhu.edu/map.html](https://coronavirus.jhu.edu/map.html). ## Author Contributions JO and GD developed the computer code and performed the simulations. SB helped derive the governing equations and generated analytical solutions used to verify computations. SBB researched COVID-19 databases and performed operations to optimize data use by the model. MM designed the study. All authors contributed to the writing and editing of the manuscript. ## Competing Interest Statement The authors have no competing interests to report. ## Footnotes * Jenna.osborn{at}fda.hhs.gov * Sberma13{at}terpmail.umd.edu * Sara.benderbier{at}gmail.com * Gavin.dsouza{at}fda.hhs.gov * Matthew.myers{at}fda.hhs.gov * Received November 9, 2020. * Revision received January 26, 2021. * Accepted January 27, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license ## References 1. 1.Stutt, R. O. J. H., Retkute, R., Bradley, B., Gilligan, C. A. & Colvin, J. A modelling framework to assess the likely effectiveness of facemasks in combination with ‘lock-down’ in managing the COVID-19 pandemic. Proc. R. Soc.A 476: 20200376. doi:10.1098/rspa.2020.0376 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.2020.0376&link_type=DOI) 2. 2.Giordano, G., Blanchini, F., Bruno, R., Colaneri, P., Di Filippo, A., Di Matteo, A. & Colaneri, M. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nature Medicine Let. doi:10.1038/s41591-020-0883-7 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0883-7&link_type=DOI) 3. 3.Bertozzi, A. L., Franco, E., Mohlerd, G., Shorte, M. B. & Sledge, D. The challenges of modeling and forecasting the spread of COVID-19. Proc. Nat. Acad. Sci 117, 16732–16738 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzI5LzE2NzMyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMjcvMjAyMC4xMS4wOS4yMDIyODY4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 4. 4.Cooper, I., Mondal, A., & Antonopoulos, C. A SIR model assumption for the spread of COVID-19 in different communities. Chaos, Solitons, and Fractals, 2020. doi:10.1016/j.chaos.2020.110057 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chaos.2020.110057&link_type=DOI) 5. 5.Myers, M., Yan, J., Hariharan, P., Guha, S. A Mathematical Model for Assessing the Effectiveness of Protective Devices in Reducing Risk of Infection by Inhalable Droplets. Mathematical Medicine and Biology, October 2016, [https://doi.org/10.1093/imammb/dqw018](https://doi.org/10.1093/imammb/dqw018). 6. 6.Siegenfeld, A. F., Taleb, N. N. & Bar-Yam, Y. What models can and cannot tell us about COVID-19. Proc. Nat. Acad. Sci. U.S.A. 117, 16092–16095 (2020). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzI4LzE2MDkyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMjcvMjAyMC4xMS4wOS4yMDIyODY4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 7. 7.Stilianakis, N. I. & Drossinos, Y. Dynamics of infectious disease transmission by inhalable respiratory droplets. J. Royal Society Interface 7, 1355–1366 (2010). 8. 8.Yan, J., Hariharan, P., Guha, S., Myers, M. Modeling the effectiveness of respiratory protective devices in reducing influenza outbreak. Risk Analysis, September 2018. doi:10.1111/risa.13181. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/risa.13181&link_type=DOI) 9. 9.Johns Hopkins Coronavirus Resource Center (2020) [https://coronavirus.jhu.edu/map.html](https://coronavirus.jhu.edu/map.html). 10. 10.Howard, J., Huang, A., Lik, Z., et al., Face Masks Against COVID-19: An Evidence Review. Proc. Nat. Acad. Sci. U.S.A. (2020). DOI: 10.20944/preprints202004.0203.v2 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20944/preprints202004.0203.v2&link_type=DOI) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/inline-graphic-1.gif [5]: /embed/graphic-4.gif [6]: /embed/graphic-5.gif [7]: /embed/graphic-6.gif [8]: /embed/graphic-7.gif [9]: /embed/inline-graphic-2.gif [10]: /embed/graphic-8.gif [11]: /embed/graphic-9.gif [12]: /embed/graphic-10.gif [13]: /embed/graphic-11.gif [14]: /embed/graphic-12.gif [15]: /embed/graphic-13.gif [16]: /embed/graphic-14.gif [17]: /embed/graphic-15.gif [18]: /embed/graphic-16.gif [19]: /embed/graphic-17.gif [20]: /embed/graphic-18.gif [21]: /embed/graphic-19.gif [22]: /embed/graphic-20.gif [23]: /embed/inline-graphic-3.gif [24]: /embed/inline-graphic-4.gif [25]: /embed/inline-graphic-5.gif [26]: /embed/inline-graphic-6.gif