Face Masks, Public Policies and Slowing the Spread of COVID-19: Evidence from Canada ==================================================================================== * Alexander Karaivanov * Shih En Lu * Hitoshi Shigeoka * Cong Chen * Stephanie Pamplona ## Abstract We estimate the impact of mask mandates and other non-pharmaceutical interventions (NPI) on COVID-19 case growth in Canada, including regulations on businesses and gatherings, school closures, travel and self-isolation, and long-term care homes. We partially account for behavioural responses using Google mobility data. Our identification strategy exploits variation in the timing of indoor face mask mandates staggered over two months in the 34 public health districts in Ontario, Canada’s most populous province. We find that mask mandates are associated with a 25 percent or larger weekly reduction in new COVID-19 cases in July and August, relative to the trend in absence of mask mandate. Additional analysis with province-level data provides corroborating evidence. Counterfactual policy simulations suggest that mandating indoor masks nationwide in early July could have reduced the number of new cases in Canada by 25 to 40 percent in mid-August, which corresponds to 700 to 1,100 fewer cases per week. **JEL codes** I18, I12, C23 Keywords * COVID-19 * face mask mandates * non-pharmaceutical interventions * counterfactuals ## 1 Introduction When government policies to stem the spread of COVID-19 were introduced in early 2020, the best available evidence supporting them came from research on previous epidemics, epidemiological modeling and case studies (OECD, 2020). Even when the efficacy of a given precaution in reducing COVID-19 transmission has been established, doubts regarding the usefulness of specific policy measures may persist because of uncertainty regarding adherence to the rules and other behavioural responses. For example, even though there is significant agreement in the medical literature that respiratory transmission of COVID-19 is the dominant vector (Meyerowitz et al., 2020), and many clinical studies have shown that face masks reduce the transmission of COVID-19 and similar diseases (Chu et al. 2020, Prather et al. 2020, Leung et al. 2020, Greenhalgh et al. 2020), a mask mandate may not be effective in practice if it fails to increase the prevalence of mask wearing (compliance) or if it leads to increased contacts due to a false sense of security. The relatively low economic cost of mask mandates relative to other COVID-19 containment measures has generated keen interest worldwide for studying their effectiveness. This interest has been compounded by the substantial variation in official advice regarding mask use, with national health authorities and the World Health Organization giving inconsistent or contradictory recommendations over time, ranging from ‘not recommended’ to ‘mandatory’. Figure B1 in the Appendix plots self-reported mask usage in select countries (Canada, United States, Germany and Australia) in the left panel, and across Canadian provinces in the right panel. The figure shows large differences in mask usage, both across countries and within Canada.1 An added challenge is to disentangle the impact of mask mandates from that of other policies, behavioural responses, or third factors (Chernozhukov et al. 2020, Mitze et al. 2020). Given the absence of large-scale randomised controlled trials on mask mandates (Howard et al., 2020), observational studies like ours are essential for informing health policy and public opinion, by formally analyzing the relationship between the various policy measures and the rate of propagation of COVID-19. We estimate and quantify the impact of mask mandates and other non-pharmaceutical interventions (NPI) on the growth of the number of COVID-19 cases in Canada. The Canadian data have the important advantage of allowing two complementary approaches to address our objective. First, we estimate the effect of mask mandates by exploiting within-province geographic variation in the timing of indoor face mask mandates in the 34 public health regions (PHUs) in Ontario, Canada’s most populous province with roughly 40% (15 mln) of the country’s population (Statistics Canada, 2020). The advantage of this approach is that it uses variation over a relatively small geographic scale (PHU), holding all province-level policies and factors constant. The adoption of indoor face mask mandates in the 34 sub-regions was staggered over approximately two months, creating sufficient intertemporal policy variation. Second, we evaluate the impact of NPIs in Canada as a whole, using the variation in the timing of policies in the country’s ten provinces. By studying inter-provincial variation in the COVID-19 policy measures, we can analyze the impact of not only mask mandates, but also other NPIs, for which there is little or no variation across Ontario’s PHUs, including regulations on businesses and gatherings, schooling, travel and long-term care. In addition, the province-level data include both the ‘closing’ period (March-April) and the gradual ‘re-opening’ period (May-August), providing variation from both the imposition and the relaxation of policies. Our panel-data estimation strategy broadly follows the approach of Chernozhukov, Kasahara and Schrimpf (2020), hereafter CKS (2020), adapted to the Canadian context. We allow for behavioural responses (using Google Community Mobility Reports geo-location data as proxy for behaviour changes and trends), as well as lagged outcome responses to policy and behavioural changes. Our empirical approach also allows past epidemiological outcomes to impact current outcomes, as information variables affecting unmeasured behaviour and policy, or directly, as in the SIR model framework. We find that, two weeks after implementation, mask mandates are associated, on average, with a reduction of 29 to 37 log points in the weekly case growth rate, which can be interpreted as a 25 to 31% weekly reduction in new diagnosed COVID-19 cases in Ontario, relative to the trend in absence of mask mandate. We find corroborating evidence in the province-level analysis, with a 36 to 46% weekly reduction in cases relative to the no-mandate trend, depending on the empirical specification. Furthermore, using additional survey data (Jones et al. 2020), we show that mask mandates increased self-reported mask usage in Canada by 30 percentage points shortly after implementation, suggesting that the policy has a significant impact on behaviour. Jointly, our results suggest that mandating mask wearing in indoor public places can be a powerful policy tool to slow the spread of COVID-19, with little associated economic disruption at least in the short run.2 Counterfactual policy simulations using our empirical estimates suggest that mandating indoor masks Canada-wide in early July could have reduced weekly new cases in the country by 25 to 40% on average by mid-August relative to the actually observed numbers, which corresponds to 700 to 1,100 fewer new cases per week. In addition, we find that the most stringent policy restrictions on businesses and gatherings observed in the data are associated with a weekly decrease of 48 to 57% in new cases (65 to 85 log point reduction in the case growth rate), relative to absence of restrictions. The business/gathering estimates are, however, noisier than our estimates for mask mandates and do not retain statistical significance in all specifications; they appear driven by the smaller provinces and the re-opening period (May to August). School closures and travel restrictions are associated with a large decrease in weekly case growth in the closing period (March-April). Our results on business/gathering regulations and school closures suggest that relaxed restrictions and the associated increase in business or workplace activity and gatherings (including retail, restaurants and bars) or school re-opening can offset in magnitude, in whole or in part, the estimated effect of mask mandates on COVID-19 case growth, both in our sample and subsequently. A further contribution of this research is that we assemble, from official sources only, and make publicly available a complete data set of COVID-19 cases, deaths, tests and policy measures in all 10 Canadian provinces.3 Specifically, we constructed, from official public health orders and announcements, time series for 17 policy indicators regarding face masks, regulations on businesses and gatherings, school closures, travel and self-isolation, and longterm care homes. Our paper relates most closely to two other recent empirical studies on the effects of mask mandates using observational data.4 CKS (2020) and Mitze et al. (2020) study the effect of mask mandates in the United States and Germany, respectively. CKS (2020), whose estimation strategy we broadly follow, exploit U.S. state-level variation in the timing of mask mandates for employees in public-facing businesses, and find that these mandates are associated with 9 to 10 percentage point reduction in the weekly growth rate of cases. This is substantially smaller that our estimates, possibly because the mask mandates that we study are much broader: they apply to all persons rather than just employees, and most apply to all indoor public spaces rather than just businesses. Mitze et al. (2020) use a synthetic control approach and compare the city of Jena and six regions in Germany that adopted a face mask policy in early to mid April 2020, before their respective state mandate. They find that mandatory masks reduce the daily growth rate of cases by about 40%. Our paper has several advantages compared to the above two studies. First, we exploit both regional variation within the same province (like Mitze et al., 2020 but without synthetic controls) and provincial variation in the whole country (like CKS, 2020) and find similar results, which strengthens the validity of our findings. Second, we additionally show that self-reported mask usage has increased after introducing mask mandates. We view this “firststage” result on mask usage as informative, since the effectiveness of any NPI or public policy critically depends on the compliance rate. Moreover, this result mitigates possible concerns that the large estimated mask mandate effect on COVID-19 case growth may be caused by factors other than the mask policy. Third, an important difference between our paper and CKS (2020) is that we evaluate the effect of *universal* (or *community*) mandatory indoor mask wearing for the public rather than the effect of mandatory mask wearing for *employees only*.5 While other factors such as differences in mask wear compliance between Canada and the U.S. may also contribute to the different estimated magnitude of the policy impact, our results suggest that more comprehensive mask policies can be more effective in reducing the case growth rate. ### Other Related Literature Abaluck et al. (2020) study the effectiveness of universal adoption of homemade cloth face masks and conclude that it could yield large benefits, in the $3,000–$6,000 per capita range, by slowing the spread of the virus. The analysis compares countries with pre-existing norms that sick people should wear masks and countries without such norms.6 In the medical literature, Prather et al. (2020) conclude that masks can play an important role in reducing the spread of COVID-19. Howard et al. (2020) survey the medical evidence on mask efficiency and recommend public use of masks in conjunction with existing hygiene, distancing, and contact tracing strategies. Greenhalgh et al. (2020) provide evidence on mask usage during non-COVID epidemics (influenza and SARS) and conclude that even limited protection could prevent some transmission of COVID-19. Leung et al. (2020) study exhaled breath and coughs of children and adults with acute respiratory illness and find that the use of surgical face masks could prevent the transmission of the human coronavirus and influenza virus from symptomatic individuals. Meyerowitz et al. (2020), a recent comprehensive review of the evidence on COVID-19 transmission, conclude that there is strong evidence, from case and cluster reports, that respiratory transmission is dominant, with proximity and ventilation being key determinants of transmission risk, as opposed to direct contact or fomite transmission. Our paper also complements recent work on COVID-19 policies in Canada.7 Mohammed et al. (2020) use survey data to study the effect of changes in mask-wear policy recommendations, from discouraged to mandatory, on the rates of mask adoption and public trust in government institutions. They show that Canadians exhibit high compliance with mask mandates and that trust in public health officials remained consistent over time. Yuksel et al. (2020) use an outcome variable constructed from Apple mobility data along with weather data and lagged COVID-19 cases or deaths as dependent variables to study compliance with social distancing measures.8 ## 2 Data We use three main data sources, respectively for epidemiological variables, NPI and mask mandates, and behavioural responses. The time period is from the start of detected community transmission in Canada in March to mid-August, 2020. We located and accessed the original official sources to collect a complete data set of COVID-19 cases, deaths, tests and policy measures in all ten Canadian provinces.9 In addition, our data include cases and policy measure indicators for each of the 34 public health units (PHUs) in Ontario. A detailed description is provided in the data files shared at the project’s Github repository. Implementation dates of NPIs and other public policies were collected from government websites, announcements, public health orders and staged re-opening plans, collected from the respective official sources. In the national data, the raw policy measures data contain the dates or enactment and relaxation (if applicable) of 17 policy indicators including: mandatory mask wear; closure and re-opening of retail and non-essential businesses, restaurants, recreation facilities, and places of worship; school closures; limits on events and gatherings; international and domestic travel restrictions and self-isolation requirements; restrictions on visits and staff movement in long-term care homes. All policy indicator variables are defined in Table C1 in the Appendix.10 Since many of these indicators are highly correlated with each other, we combine them into five policy aggregates in the empirical analysis (see Table A17 and Section 3.2). In the Ontario PHU data, the implementation dates of mask mandates and the relaxation dates of policies for businesses and gatherings vary across PHUs. Decisions about the former were made at the PHU level, while decisions about the latter were made by the province, which classified PHUs into three groups, with some exceptions. Regarding behavioural responses, we use the Google COVID-19 Community Mobility Reports, which summarize daily cellphone geo-location data for each province as indices calculated relative to the median value for the same day of the week in the five-week baseline period January 3 to February 6, 2020.11 In Ontario, these location data are available for each of the 51 first-level administrative divisions (counties, regional municipalities, single-tier municipalities and districts).12 ## 3 Empirical method We follow the approach of CKS (2020), but modify and adapt it to the Canadian context. The empirical strategy uses the panel structure of the outcome, policy and behavioural proxy variables, and includes lags of outcomes as information, following the causal paths suggested by the epidemiological SIR model (Kermack and McKendrick, 1927). Specifically, we estimate the effect of policy interventions on COVID-19 outcomes while controlling for information and behaviour. In contrast to CKS (2020) and Hsiang et al. (2020), who study variation in NPIs across U.S. states or across countries, our identification strategy exploits policy variation at the sub-provincial level (Ontario’s PHUs) in addition to cross-province variation, and our data captures both the closing down and gradual re-opening stages of the epidemic. ### 3.1 Estimation strategy The main data used in our empirical analysis are summarized below; Section 3.2 describes the variables in detail. Everywhere *i* denotes health region (PHU) for the Ontario analysis and province for the national analysis, and *t* denotes time measured in days. 1. Outcomes, *Y**it* – the growth rate of weekly cases or deaths. 2. Policy/NPIs, *P**it* – for the national analysis, five policy aggregates by province and date; for the Ontario analysis, two policy variables (mask mandates and regulations on businesses and gatherings) by PHU and date. 3. Behavioural responses, *B**it* – proxied by Google mobility data capturing changes in people’s geo-location relative to a baseline period in January-February. 4. Information, *I**it* – lagged outcomes, i.e., level or growth rate of cases (or deaths). We also consider a specification that includes past cases (deaths) and case (deaths) growth at the national level as additional information variables. 5. Controls, *W**it* – province or PHU fixed effects, growth rate of weekly new COVID-19 tests, and a time trend. To assess and disentangle the impact of NPIs and behavioural responses on COVID-19 outcomes, we estimate the following equation: ![Formula][1] where *l* denotes a time lag measured in days. Equation (1) models the relationship between COVID-19 outcomes, *Y**it* and lagged policy measures, *P**it*−*l*, lagged behaviour, *B**it*−*l*, and information (lagged outcomes), *I**it*. For case growth as the outcome, we use a 14-day lag, *l* = 14. For deaths growth as the outcome, we use *l* = 28.13 The choice of these lags is discussed in detail in Appendix D, and alternative lags are explored in Section 4.3. By including lagged outcomes, our approach allows for possible endogeneity of the policy interventions *P**it*, that is, the introduction or relaxation of NPIs based on information about the level or growth rate of cases or deaths. Also, past cases may be correlated with (lagged) government policies or behaviours that may not be fully captured by the included policy and behaviour variables. In Appendix Table A18, we also report estimates of the following equation: ![Formula][2] which models the relationship between policies *P**it*, information, *I**it* (weekly levels or growth of cases or deaths) and behaviour, *B**it*. It is assumed that behaviour reacts to the information without a significant lag. We find strong correlation between policy measures and the Google mobility behavioural proxy measure (see Table A18). Equation (1) captures both the direct effect of policies on outcomes, with the appropriate lag, as well as the potential indirect effect on outcomes from changes in behaviour captured by the geo-location trends proxy, *B**it*−*l*. In Appendix Tables A19 and A20, we also report estimates of equation (1) without including the behavioural proxy, that is, capturing the total effect of policies on outcomes. Since our estimates of the coefficient *α* in equation (1) are not significantly different from zero, the results without controlling for the behavioural proxy *B**it*−*l* are very similar to those from estimating equation (1). ### 3.2 Data variables and descriptive analysis #### Outcomes Our main outcome of interest is the weekly growth rate of new positive COVID-19 cases as defined below.14 We use weekly outcome data to correct for the strong day-of-the-week effect present in COVID-19 data.15 Weekly case growth is a metric that can be helpful in assessing trends in the spread of COVID-19, and it is highlighted in the World Health Organization’s weekly epidemiological updates (World Health Organization, 2020). Specifically, let *C**it* denote the cumulative case count up to day *t* and define Δ*C**it* as the weekly COVID-19 cases reported for the 7-day period ending at day *t*: ![Formula][3] The (log) weekly case growth rate, *Y**it* is then defined as: ![Formula][4] that is, the week-over-week growth in cases in region *i* ending on day *t*.16 The weekly death growth rate is defined analogously, using cumulative deaths data. #### Policy In the Ontario analysis, we exploit regional variation in the timing of indoor mask mandates staggered over two months in the province’s 34 regions (“public health units” or PHUs). Figure 1 displays the gradual introduction of mask mandates in the 34 PHUs. The exact implementation dates of the mask mandates are reported in Table C2. Mandatory indoor masks were introduced first in the Wellington-Dufferin-Guelph PHU on June 12 and last in the Northwestern PHU on August 17.17 ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F1) Figure 1: Ontario - mask mandates over time Notes: There are a total of 34 public health units (PHU) in Ontario. See Table C2 for the exact date of mask mandate implementation in each PHU. In the province-level analysis, we assign numerical values to each of 17 policy indicators listed in Table C1 in Appendix C. The values are on the interval [0,1], with 0 meaning absence or lowest level of restrictions and 1 meaning maximal restrictions. A policy value between 0 and 1 indicates partial restrictions, either in terms of intensity (see Table C1 for more details and definitions) or geographical coverage (in large provinces). The numerical values are assigned at the daily level for each region (PHU or province), while maintaining comparability across regions. Many NPIs were implemented at the same time, both relative to each other and/or across regions, especially during the March closing-down period. This causes many of the policy indicators to be highly correlated with each other (see Appendix Table A4). To avoid multi-collinearity issues, we therefore group, via simple averaging, the 17 policy indicators into 5 policy aggregates: (i) *travel*, which includes international and domestic travel restrictions and self-isolation rules; (ii) *school*, which is an indicator of provincial school closure (including Spring and Summer breaks); (iii) *business/gathering*, which comprises regulations and restrictions on non-essential businesses and retail, personal businesses, restaurants, bars and nightclubs, places of worship, events, gyms and recreation, and limits on gathering; (iv) *long-term care (LTC)*, which includes NPIs governing the operation of long-term care homes (visitor rules and whether staff are required to work on a single site) and (v) *mask* which takes value 1 if an indoor mask mandate has been introduced, 0 if not, or value between 0 and 1 if only part of a province has enacted the policy.18 The five policy aggregates are constructed at the daily level and capture both the closing-down period (an increase in the numerical value from 0 toward 1) and the re-opening period (decrease in the numerical value toward zero). In comparison, the policy indicators compiled by Raifman (2020) for the USA used in CKS (2020) are binary “on (1)”/”off (0)” variables.19 For consistency with the weekly outcome and information variables and the empirical model timing, we construct the policy aggregates ![Graphic][5] used in the regressions (where *j* = 1,…, 5 denotes policy type) by taking a weekly moving average of the raw policy data, from date *t* − 6 to date *t*. Figure 2 plots the values of the 5 policy aggregates over time for each of the 10 provinces. Travel restrictions, school closures and business closures were implemented in a relatively short period in the middle of March. There is some variation in the travel policy aggregate since some Canadian provinces (the Atlantic provinces and Manitoba) implemented interprovincial domestic travel or self-isolation restrictions in addition to the federal regulations regarding international travel. Restrictions on long-term care facilities were introduced more gradually. In the re-opening period (May-August), there is more policy intensity variation, especially in the business and gatherings category, as the different provinces implemented their own re-opening plans and strategies. Mask mandates were introduced in Ontario starting from June in a few smaller PHUs and in early July in some of the most populous PHUs, Toronto, Ottawa and Peel (see Appendix Table C2). In Quebec, indoor masks were mandated province-wide on July 18. Nova Scotia (province-wide) and Alberta’s two main cities Calgary and Edmonton implemented mask mandates on July 31 and August 1, respectively. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F2) Figure 2: Policy aggregates – Canada Notes: The figure plots the numerical values of the 5 policy aggregates (Mask, Business/gathering, School, Travel and Long-term care, LTC) over time, for each of the 10 provinces. The mask policy values for ON reflect the gradual adoption of mandates (see Fig. 1) and the respective PHUs population sizes. There are two empirical challenges specific to the Canadian context and data. The first challenge is the presence of small provinces or sub-regions with very few COVID-19 cases or deaths. In Section 4.3, we perform robustness checks using different ways of handling the observations with very few cases (in particular zero cases). The second data limitation is that there are only 10 provinces in Canada and 34 public health units in Ontario, unlike the 51 U.S. jurisdictions in CKS (2020). To account for the resulting small number of clusters in the estimation, we compute and report wild bootstrap standard errors and p-values, as proposed by Cameron et al. (2008).20 On the positive side, our data has the advantage of a longer time horizon (March to August) and non-binary, more detailed policy variables compared to Raifman et al. (2020). #### Behaviour proxy We follow CKS (2020) and other authors in interpreting the location change indices from the Google Community Mobility reports as proxies for changes in people’s behaviour during the pandemic, keeping in mind that location is only one aspect of behaviour relevant to COVID-19. The general pattern in the data (see Figure B3) shows sharply reduced frequency of recorded geo-locations in shops, workplaces and transit early in the pandemic (March), with a subsequent gradual increase back toward the baseline (except for transit), and a flattening out in July and August. Several of the six location indicators (retail, grocery and pharmacy, workplaces, transit, parks and residential) are highly correlated with each other (see Tables A1 and A2) and/or contain many missing observations for the smaller PHUs and provinces. To address these data limitations and the possible impact of collinearity on the estimation results, we use as proxy for behavioural changes the simple average of three mobility indicators: “retail”, “grocery and pharmacy” and “workplaces”.21 To be consistent with the weekly outcome variables and to mitigate day-of-the-week geo-location variation, we construct the Behaviour proxy *B**it* by taking a weekly moving average of the ![Graphic][6] (retail + grocery and pharmacy + workplaces) data, from date *t* − 6 to date *t*.22 As a result, our empirical analysis uses weekly totals (for cases, tests and deaths) or weekly moving averages (for policies and the behavior proxy) for all variables recorded on daily basis. Tables A3 and A4 display the correlation between the behaviour proxy *B**it* and the five policy aggregates ![Graphic][7]. The behaviour proxy and mask mandate variable are not highly correlated, suggesting that the effect of mask mandates on COVID-19 outcomes is likely not dependent of location behaviour changes. #### Information We use the weekly cases and case growth variables defined above, Δ*C**it* and *Y**it*, to construct the information variables *I**it* used in estimation equation (1). Specifically, we use as information the lagged value of the weekly case growth rate, Δ log(Δ*C**it*−*l*) and the log of past weekly cases, log(Δ*C**it*−*l*). We also use the lagged provincial (in the Ontario analysis) or national (in the Canada analysis) case growth rate and log of weekly cases as additional information variables in some specifications. A two-week information lag, *l* = 14, is used in all baseline analysis. In the supplementary analysis using the death growth rate as the outcome, we use information on past deaths and a four-week lag (see Section 4.6). #### Control variables In all regressions, we control for region fixed effects (PHU or province) and the weekly COVID-19 tests growth rate Δ log(Δ*T**it*), where *T**it* denotes cumulative tests in region *i* until date *t* and Δ*T**it* is defined analogously to Δ*C**it*. We also include a time trend: our baseline results use a cubic polynomial in days, and we also report results with no time trend and with week fixed effects. In robustness checks, we also include news or weather variables as controls (see Section 4.3). #### Time period We use the period May 15 to August 13 for the analysis with Ontario PHU data and the period March 11 to August 13 for the analysis with provincial data. The end date reflects the data availability at the time of empirical analysis and writing. The start date for the Ontario sample (May 15) is approximately two weeks after the last restrictive measures were implemented and four weeks before the first mask mandate was introduced in Ontario. Sensitivity checks with different initial dates (May 1, June 1 and June 15) are reported in Section 4.3, with our main results remaining robust. The initial date for the Canada sample (March 11) was chosen as the first date on which each province reported at least one COVID-19 test (so that cases could be potentially reported). Alternative initial dates are also considered in Section 4.3. ## 4. Results ### 4.1 Mask mandates in Ontario public health regions We start with a simple graphical illustration of the effect of mask mandates on COVID-19 cases growth. Figure 3 displays the average log case growth, *Y**it* = Δ log(Δ*C**it*) in Ontario PHUs with and without mask mandates. On average, the PHUs with a mask mandate two weeks prior have lower case growth than the PHUs without a mask mandate two weeks prior. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F3) Figure 3: Ontario - mask mandates and weekly case growth Notes: The figure plots the average log weekly case growth Δlog(ΔC) in the PHUs with mask mandate (blue) and without (red) mask mandate 14 days prior. Table 1 shows the estimates of equation (1), in which we control for other policies, behaviour, and information, as explained in Section 3.1.23 We report wild bootstrap p-values clustered at the PHU level to account for the small number of clusters.24 The odd-numbered columns in Table 1 use lagged cases and lagged case growth at the PHU level as information; the even-numbered columns also include lagged cases and lagged case growth at the province level as additional information variables.25 View this table: [Table 1:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T1) Table 1: Main Results - Ontario public health regions We present estimates of equation (1) from three specifications that handle time effects differently. Columns (1) and (2) in Table 1 are the most basic specifications, without a time trend. To control for additional possible province-wide factors affecting the spread of COVID-19, e.g. income support policies, adaptation to the pandemic over time, or so-called COVID fatigue, we also estimate (1) including a cubic time trend in days from the beginning of the sample, in columns (3) and (4) of Table 1, and including week fixed effects, in columns (5) and (6). The estimates in Table 1 imply that, controlling for other policies, information, testing, and geo-location mobility behaviour, mandatory indoor face masks are associated with a decrease of 29.1 to 37.6 log points (*p* < 0.05) in the weekly growth rate of new cases Δ*C**it*, two weeks after their implementation. This result can be interpreted as a 25–31% weekly reduction in new cases, relative to the trend without mask mandate.26 Columns (3)-(6) of Table 1 show that our estimates of the mask mandate policy remain robust to the inclusion of cubic time trend or week fixed effects. The magnitude of the mask policy estimates is not very sensitive to whether lagged province-level data are included as additional information. The results in Table 1 indicate that indoor mask mandates can be a powerful preventive measure in the COVID-19 context. Our estimates of the mask mandate impact in Ontario’s PHUs are larger than the 9–10 percentage point reduction in case growth estimated by CKS (2020) for the U.S. A possible explanation is that Ontario’s mask mandates are more comprehensive: we evaluate the effect of universal indoor mask-wearing for the public rather than the effect of mask-wearing for *employees only* in CKS (2020). Differences in the compliance rate may also be a factor; we discuss this potential mechanism in Section 4.4. The results in Table 1 also show a statistically significant negative association between information (log of past cases, log(ΔC)_14) and current weekly case growth (*p* < 0.01 in all specifications), indicating that a higher level of cases two weeks prior is correlated with lower current case growth. While *B**it* allows for certain behavioural responses to information, the negative estimate on past cases, log(ΔC)_14 in Table 1 suggests that our location-based proxy may not capture other important aspects of behaviour, e.g., frequent hand-washing or physical distancing. In fact, our coefficient estimate on the behaviour proxy *B**it* is very close to zero (both in Table 1 and in Section 4.2’s province-level results), unlike in CKS (2020).27 In Appendix Table A18, we find strong contemporaneous correlations between the policy measures, log cases, and the Google mobility behavioural proxy from estimating equation (2). This suggests that the information (lagged cases) and the lagged policy variables included in equation (1) may absorb lagged behavioural responses proxied by *B**it*−*l* or other latent behavioural changes not captured by *B**it*−*l*. ### 4.2 Province-level results We next evaluate the impact of NPIs on COVID-19 case growth in Canada as a whole by exploiting variation in the timing of policy interventions across the 10 provinces. Here, we also examine NPIs for which there is no variation across Ontario’s PHUs (i.e., school closures, travel, and long-term care regulations), in addition to mask mandates. In addition, the provincial data contain variation in the timing and intensity of policy changes in both the closing and re-opening phases, allowing us to study both the imposition and relaxation of restrictions. As in the Ontario analysis, we begin with a graphical illustration of mask mandates and COVID-19 case growth across Canadian provinces, from March 11 to August 13, 2020. Figure 4 plots the average log weekly case growth in the provinces with and without mask mandates two weeks prior. While mask mandates are implemented relatively late in our sample period, average case growth in the provinces with a mask mandate (Ontario and Quebec) diverged from the average case growth in the provinces without a mandate begin roughly two to four weeks after the mandates are imposed.28 ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F4) Figure 4: Canada - mask mandates and weekly case growth Notes: The figure plots the average weekly case growth Δ log(ΔC) in the provinces with mask mandate (blue) vs. without mask mandate (red) 14 days prior. Table 2 displays the estimates of equation (1) for weekly case growth, along with wild bootstrap p-values, clustered at the province level (see Table A9 for other methods of computing the standard errors). The odd-numbered columns use lagged cases and lagged case growth at the provincial level as information while the even-numbered columns additionally include lagged cases and case growth at the national level. View this table: [Table 2:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T2) Table 2: Main Results – Canada As in Table 1, we present in Table 2 estimates from three specifications: no time trend (columns (1)-(2)), with cubic time trend in days (columns (3)-(4)) and with week fixed effects (columns (5)-(6)). The most robust result is the estimated effect of mask mandates: they are associated with a large reduction in weekly case growth of 45 to 62 log points, which corresponds to a 36 to 46% weekly reduction in new cases across the different specifications, relative to the trend in absence of mandate. The estimates are statistically significantly different from zero in all cases, with a p-value of less than 0.001 in columns (1)-(4). It is reassuring that these province-level results regarding mask mandates are consistent with our results for Ontario PHUs in the previous section. Table 2 further shows that restrictions on businesses and gatherings are associated with a reduction in the weekly case growth of 65 to 85 log points or, vice versa, that relaxing business/gathering restrictions is associated with higher case growth. The estimate corresponds to a 48 to 57% weekly decrease in new cases in our sample period, relative to the trend in absence of mandate. The business/gathering estimates are, however, noisier than our estimates for the mask mandate policy and do not retain statistical significance in the specifications with week fixed effects (*p* = 0.15 and 0.14). Tables A8 and A15 in the Appendix further suggest that the results on business and gathering NPIs are driven by the smaller provinces and the re-opening period (May to August). Still, these findings suggest that relaxed restrictions and the associated increase in business and workplace activity or gatherings (including restaurants, bars and retail) can be an important offsetting factor for the estimated impact of mask mandates on COVID-19 case growth, both in our sample and at later dates. We also find that school closures (the School_14 variable in Table 2) are negatively associated with case growth. However, the estimates are statistically significant from zero only in the specifications with cubic time trend (columns (3) and (4)). As seen in Figure 2, provincial school closures occurred in a very short time interval in March, so we lack statistical power to separately identify their effect from that of other NPIs (especially travel-related restrictions). Hence, we interpret this result with caution. Finally, as in Table 1, past cases, log(Δ*C*) 14, are negatively and statistically significantly associated with current weekly case growth in columns (1)-(4). Since the specification with cubic time trend in Tables 1 and 2 allows for non-monotonic aggregate time trends in case growth in a parsimonious way, we choose it as our baseline specification for robustness checks and counterfactual exercises.29 ### 4.3 Robustness #### Policy collinearity A possible concern about our data for the national analysis is that some NPIs (e.g., international travel restrictions or closing of schools) were implemented within a very short time interval.30 Thus, we may lack sufficient regional variation to distinguish and identify the separate effect of each policy.31 Collinearity could also affect the standard errors and the signs of the estimated coefficients. To check robustness with respect to potential collinearity in the policy variables, Tables A7 and A10 report estimates from our baseline specifications with cubic time trend, omitting one policy at a time, for Ontario and Canada respectively. First, the mask mandate estimates are hardly affected by omitting any of the other policies. This is expected since mask mandates were imposed during a period where other NPIs changed little (see Figure 2). Similarly, the effects of business/gathering regulations and school closures in Table A10 are not sensitive to omitting other policies one at a time, which suggests that there is sufficient statistical power and variation to identify them in the national analysis. #### Treatment of zero weekly cases Another possible concern is that the dependent variable *Y**it* = Δ log(Δ*C**it*) is not well defined when the weekly case totals Δ*C**it* or Δ*C**it*−7 is zero. As in CKS (2020), we replace log(0) with −1 in our baseline analysis.32 We now check the robustness of our estimates to alternative treatments of zero weekly cases. For easier comparison, the first two columns in Table A5 repeat the baseline results from columns (3) and (4) in Table 1 for Ontario. Our main results on mask mandates in Ontario PHUs are robust to replacing log(0) with 0 and to adding 1 to all Δ*C**it* observations before taking logs, as shown in columns (3)-(6) of Table A5. Another way to address the issue of PHUs with very few cases is to estimate a weighted least squares regression in which the PHUs are weighted by population. Columns (7) and (8) in Table A5 show that the resulting mask mandate estimate has a slightly smaller magnitude and, due to the reduced effective sample size, weaker statistical significance. Similarly, Table A8 shows that our province-level estimates, in particular for mask mandates, are also robust to the same alternative specifications as above. In columns (9) and (10) of Table A8, we also restrict the sample to only the largest 4 provinces (British Columbia, Ontario, Quebec and Alberta), for which only 0.3% of observations (2 out of 624) are zeroes. Again, the estimated mask mandate effect changes very little. #### Alternative dates Figure B4 shows that our estimates and confidence intervals for the effect of mask mandates in the Ontario baseline regressions do not vary much with the initial date of the sample. Similarly, Figure B5 shows that, in the national analysis, our results about mask mandates and business/gathering restrictions are also robust to alternative initial dates. #### Alternative lags We explore alternative time lags, of either shorter or longer duration, centered around the baseline value of 14 days. Figure B6 (with Ontario data) and Figure B7 (with province-level data) plot the estimates and confidence intervals from the baseline regressions and show that our mask mandate estimates remain fairly consistent for different plausible lag values. #### Omitted variables Our behaviour proxy variable (Google geo-location trends) likely misses some aspects of behaviour that could be relevant for COVID-19 transmission. One factor that may impact behaviour is weather. For example, good weather may cause more people to spend time outside, lowering the chance of viral transmission. Columns (3) and (4) in Table A11 report national estimates with lagged weather variables (daily maximum and minimum temperatures and precipitation for the largest city in each province33) as additional regressors. Our NPI estimates, in particular for mask mandates, are little changed from the baseline results in columns (1) and (2). Another possible concern is that our information variables (lagged cases and lagged case growth) may not fully capture the information based on which people react or adjust their behaviour, possibly affecting the observed weekly case growth. Columns (5) and (6) in Table A11 add a national-level “News” variable to the baseline specification. The variable is defined as the number of daily search results from a news aggregator website (*ProQuest Canadian Newsstream*) for the terms “coronavirus” or “COVID-19” (see Appendix C for more details). In column (6), the lagged news variable approaches the 10% significance level (*p* = 0.103). Our estimates on masks and business/gathering remain very close to those in the baseline. ### 4.4 Mask usage The effectiveness of any NPI or public policy crucially depends on whether and how it affects behaviour. In this section, we use self-reported survey data on mask usage to examine whether mask mandates indeed increase mask use in Canada (“first-stage” analysis). We use data from the YouGov COVID-19 Public Monitor (Jones et al. 2020), which includes multiple waves of public opinion surveys fielded regularly since early April 2020 in many countries.34 Here, we focus on inter-provincial comparison within Canada. Our variable of interest is based on the responses to the question “Thinking about the last 7 days, how often have you worn a face mask outside your home (e.g. when on public transport, going to a supermarket, going to a main road)?” The answer choices are “Always”, “Frequently”, “Sometimes”, “Rarely”, and “Not at all”. We create a binary variable taking value 1 if the response is “Always” and 0 otherwise, as well as another variable taking value of 1 if the respondent answered either “Always” or “Frequently” and 0 otherwise. We begin with a simple illustration of self-reported mask usage in Canada from April to August 2020. Figure B2 in the Appendix plots the average self-reported mask usage (the response “Always”) in the provinces with and without mask mandates.35 The figure clearly shows that self-reported mask usage is higher, by up to 50 percentage points, in the provinces with a mask mandate than in the provinces without mask mandates. Since Figure B2 does not account for compositional changes in the data, we formally estimate equation (2), using self-reported mask usage as the behavioural outcome.36 Figure 5 shows a graphical event study analysis of mask mandates and the change in mask usage. The event study approach is appropriate for the mask usage outcome variable, since the policy impact is expected to be immediate, unlike the other outcomes that we study, for which any impact is expected to occur with a lag and for which we use weekly totals or moving averages. We replace the mask policy variable in equation (2) by the interaction of a variable denoting being in the treatment group (i.e. under a mask mandate) with a series of dummies for each week, ranging from 6 weeks before the mask mandate to 5 weeks after the mask mandate (*T* = –6 to +5, where *T* = 0 is the implementation date of the mask mandate). The reference point is one week before the implementation of the mask mandate (*T* = –1), and we use the same y-axis scale in both panels. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F5) Figure 5: Event study of self-reported mask usage – Canada Notes: The data source is YouGov. The outcome is a binary variable taking value 1 if the respondent respectively answered “Always” (in the left panel) or “Always” or “Frequently” (in the right panel) to “Thinking about the last 7 days, how often have you worn a face mask outside your home?” The figure plots the estimates from a version of equation (2) where the mask policy variable is replaced by the interaction of a variable denoting being in the treatment group (imposed mask mandate) with a series of dummies for each week, ranging from 6 weeks before the mask mandate to 6 weeks after (*T* = –6 to +5, where *T* = 0 is the mandate implementation date). The reference point is 1 week before the implementation (*T* = –1). Wild bootstrap (cgmwildboot) standard errors clustered by province with 5,000 repetitions are used to construct the confidence intervals. Sample weights are used. The left and right panels of Figure 5 present the results from the event study analysis for the “Always” and “Always” or “Frequently” mask usage answers, respectively. We make several observations. First, neither panel shows a pre-trend – the estimates are close to zero before the mask mandates are implemented. This addresses the potential concern that provinces that implemented mask mandates may have had a different trend in mask usage than provinces that did not. Second, the effect of mask mandates on mask usage is immediate: an increase of roughly 20 percentage points as soon as the mask policy is implemented at (*T* = 0). Third, the effect appears persistent rather than transitory, since mask usage after *T* = 0 does not revert to its level before *T* = 0. Table 3 displays the estimates on self-reported mask usage (answer “Always”) in equation (2) along with wild bootstrap p-values clustered at the province level. The odd-numbered columns use lagged cases and lagged case growth at the provincial level as information while the even-numbered columns additionally include lagged cases and case growth at the national level. As in Table 1 and Table 2, we present estimates without time trend, including cubic time trend (in days), and including week fixed effects. Our preferred specification with cubic time trend, column (4) of Table 3, shows that mask mandates are associated with 31.5 percentage point increase in self-reported mask usage (*p* < 0.001), from a base of self-reported mask usage without mask mandate of 29.8%.37,38 View this table: [Table 3:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T3) Table 3: Self-reported mask usage – Canada These “first-stage” results show that mask mandates exhibit broad compliance in Canada and establish a basis for the significant impact of mask mandates on the spread of COVID-19 that we find. That said, given that mask mandates do not change everyone’s behaviour, our estimates in Tables 1 and 2 represent intent-to-treat effects. The full effect of the entire population shifting from not wearing to wearing masks is likely larger.39 There is a heated debate on whether community use of masks may create a false sense of security that reduces adherence to other preventive measures. We investigate this question using the YouGov survey data. As Tables A13 and A14 indicate, we find no evidence that mask mandates in Canada have had an offsetting effect on other preventive measures such as hand washing, using sanitizer, avoiding gatherings, and avoiding touching objects in public during the period we study. On the contrary, mask mandates may slightly increase social distancing in one out of the eight precaution categories (avoiding crowded areas) (*p* < 0.10).40 ### 4.5 Counterfactuals We evaluate several counterfactuals in which the actual mask policy in a province or Canadawide is replaced by an alternative policy. Letting *t* be the implementation date of a counterfactual policy, we set the counterfactual weekly case count, ![Graphic][8], equal to Δ*C**it* for all *t* < *t*. For each date *t* ≥ *t*, using the definition of *Y**it* in (3), we compute ![Graphic][9] and the counterfactual weekly case growth rate, ![Graphic][10], as follows: ![Formula][11] where *Ŷ**it* is the regression-fitted value of weekly case growth; *β**Mask_*14 is the coefficient estimate (−0.376 or −0.613) on the mask mandate variable Mask_14 in baseline specification (column 4) in Table 1 or 2, depending on the counterfactual; Mask*c*_14 is the counterfactual mask policy (e.g., different implementation date, wider geographic coverage or absence of mask mandate); and *β*log(Δ*C*)_14 is the coefficient estimate (−0.209 or −0.227) on lagged cases log(Δ*C*) 14 in Table 1 or 2, column 4. The coefficient *β*log(Δ*C*)_14 adjusts the counterfactual case growth rate for the negative statistically significant association between the weekly case total two weeks prior and time-*t* case growth. This feedback effect may be due to people being more careful when they perceive the risk of infection to be higher. Figures 6, 7 and 8 show results from two counterfactual policy evaluations. The first counterfactual, depicted in the left panels of the figures, assumes that masks are adopted everywhere at the earliest date observed in the data. Specifically, Figure 6 considers the counterfactual of all Ontario PHUs adopting mask mandates on June 12, while Figures 7 and 8 assume that all Canadian provinces adopt a mask mandate on July 7.41 The counterfactual line is set equal to the observed weekly case value and diverges from it after the 14-day lag from implementation (i.e., after June 26 on Figure 6 and July 21 on Figures 7 and 8). ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F6) Figure 6: Counterfactuals - Ontario public health regions Notes: The left panel assumes that mask mandates were adopted in all PHUs on June 12 (date of the first mask mandate in ON). The right panel assumes that mask mandates were not adopted in any PHU. We use the mask estimate (−0.376) from column (4) of Table 1. The counterfactual mean value (the green or red solid lines) and confidence bands (shaded area) are displayed. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F7) Figure 7: Counterfactuals – Canada (Table 1 mask estimate) Notes: The left panel assumes that mask mandates were adopted in all provinces on July 7 (the adoption date in Toronto and Ottawa). The right panel assumes that mask mandates were not adopted in any province. We use the mask estimate (−0.376) from column (4) of Table 1. The counterfactual mean value (the green or red solid lines) and confidence bands (shaded area) are displayed. ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F8) Figure 8: Counterfactuals – Canada (Table 2 mask estimate) Notes: The left panel assumes that mask mandates were adopted in all provinces on July 7 (the adoption date in Toronto and Ottawa). The right panel assumes that mask mandates were not adopted in any province. We use the mask estimate (−0.613) from column (4) of Table 2. The counterfactual mean value (the green or red solid lines) and confidence bands (shaded area) are displayed. Using our mask policy estimate from Table 1, Figure 6 shows that an early face mask mandate implemented in all Ontario PHUs could have led to an average reduction of about 300 cases per week as of August 13, holding all else equal. For Canada as a whole, a nationwide adoption of mask mandates in early July is predicted to reduce total cases per week in the country by 700 to 1,100 cases on average as of August 13, depending on whether we use the more conservative mask estimate (−0.376) from column (4) of Table 1 (see Figure 7) or the larger estimate (−0.613) from column (4) of Table 2 (see Figure 8). In all cases, the indirect feedback effect via *β*log(Δ*C*)_14 (lagged cases as information) starts moderating the decrease in cases two weeks after the initial effect of the counterfactual mask policy. In the right side panels of Figures 6, 7 and 8, we perform the opposite exercise, namely assuming instead that mask mandates were *not* adopted in any Ontario PHU or any Canadian province. Our estimates imply that the counterfactual absence of mask mandates would have led to a large increase in new cases, both in Ontario and Canada-wide, especially when using the larger mask policy coefficient estimate from Table 2 (see Figure 8). Finally, in Figure B11 in the Appendix, we evaluate the counterfactual in which British Columbia and Alberta, the third and fourth largest Canadian provinces by population, adopt province-wide mask mandates on July 15. The results, using the Mask 14 estimate from Table 2, suggest a reduction of about 300 cases per week in each province by mid-August. The counterfactual simulations assume that all other variables, behaviour and policies (except the mask policy and *t*−14 case counts) remain fixed at their values observed in the data. This is a strong assumption, but it may be plausible over the relatively short period that we analyze. In addition, the counterfactuals assume that regions without a mask mandate would react in the same way, on average, as the regions that have imposed a mandate. Therefore, these results should be interpreted with caution and only offer a rough projection of the estimated impact of mask mandates on COVID-19 cases. ### 4.6 Additional analysis #### Closing and re-opening periods We investigate whether the NPIs impact varied in different phases of the pandemic by splitting the full sample period into two sub-periods: “closing” (March 11 to May 14) and “re-opening” (May 15 to August 13). The dividing date of May 15 (referring to the NPIs in place around May 1) was chosen because very few policies were relaxed before May 1 and very few non-mask policies were tightened after May 1 in our sample period (see Figure 2). In Table A15, we report estimates and wild bootstrap standard errors using our preferred specification with cubic time trend, separately for the closing and re-opening periods. We find that the imposition of school closures and travel restrictions early in the closing period is associated with a very large subsequent reduction in weekly case growth, as can be also seen on Figure B8 – the average observed log growth rate of cases Δ log(ΔC) falls from 2.4 (ten-fold growth in weekly cases) to −0.4 (33% decrease in weekly cases) between March 15 and April 5. Long-term care restrictions are also associated with reduced case growth two weeks later during the March to May closing period. We interpret these results with caution, however, since many of the policy measures and restrictions were enacted in a brief time interval during March and there is not much inter-provincial variation (see Figure 2). No mask mandates were present in the closing period. In the re-opening period (May - August), our results in Table A15 are in line with the full-sample results for mask mandates and business/gathering regulations (Table 2), with slightly larger coefficient estimates and less statistically significant p-values, possibly due to the smaller sample. Travel restrictions and school closures are not statistically significant in the re-opening period. This is unsurprising: any relaxation of travel policies was minor and endogenous (only re-opening to safe areas within Canada), and the schools that re-opened (in parts of Quebec and part-time only in British Columbia) did so on voluntary attendance basis with small class sizes. #### Deaths We also examine the weekly death growth as an outcome variable. We only have disaggregated deaths data at the province level. We estimate equation (1) using *Y**it* = Δ log(Δ*D**it*) as the dependent variable, for each province *i* and date *t* in the sample period, where Δ*D**it* denotes weekly deaths from date *t*−6 to date *t*. We use a 28-day lag for the policy, behavior proxy, and information variables to reflect the fact that deaths occur on average about two weeks after case detection; see Appendix D for details and references.42 Table 4 reports the estimates from the same specifications as those for case growth in Table 2. In all specifications, a mask mandate is associated with a large reduction in the observed weekly deaths growth rate four weeks later (more than 90 log points reduction in growth, or equivalently more than 60% weekly reduction in deaths relative to the trend in the absence of mask mandate). These estimates are larger than our results for case growth in Tables 2, but not inconsistent with them considering the size of the standard errors. See also Figure B12, which plots the average weekly death growth in the provinces without a mask mandate four weeks prior and in Ontario, the only province with mask mandate four weeks prior in our sample period. View this table: [Table 4:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T4) Table 4: Canada – deaths growth rate and policies The robustness checks in Table A16, however, show that, unlike for case growth, the mask mandate estimates in Table 4 are not robust to weighing by population or to restricting the sample to the largest 4 provinces. This suggests that the estimated effect is largely driven by observations from the small provinces, which have a disproportionately larger number of zero or small weekly death totals.43 Furthermore, there are only 9 days with observations at least 28 days after the first significant adoption of mask mandates in Ontario, on July 7 (see footnote 41). Due to these data limitations, the relationship between mask mandates and COVID-19 deaths in Table 4 is suggestive at best and caution should be used in reporting or extrapolating from these results. That said, our findings about mask mandates and COVID-19 case growth in Tables 1 and 2 may have implications about future growth in deaths, particularly if the affected demographics become less skewed toward the young in later periods. ## 5 Conclusion The wearing of face masks by the general public remains a contentious policy issue during the COVID-19 pandemic, with large variation in official advice across jurisdictions and over time. “Conspiracy theories” and misinformation surrounding mask wearing abound in social media, sometimes fuelled by the perception that mask mandates constitute restrictions on individual freedoms. We estimate the impact of face mask mandates and other public policy measures on the spread of COVID-19 in Canada. We use both within-province and cross-province variation in the timing of mask mandates and find a robust significantly negative association between mask mandates and subsequent COVID-19 case growth: a 25 to 46% average weekly reduction in new cases in the first few weeks after adoption, depending on the data sample and empirical specification used. These results are supported by our analysis of survey data on compliance with the mask mandates, which shows that the mandates significantly increase the proportion of people reporting always wearing mask in public by around 30 percentage points. While our sample period and data do not allow us to definitively say whether the mask mandate effect persists beyond several weeks after implementation, the presence of a mandate may serve as an important behavioural anchor, especially in (temporary) periods of plateauing or declining case growth. We conclude that mandating mask wearing in indoor public spaces can be a powerful policy tool to keep COVID-19 transmission at a manageable level, especially given the mandates’ relatively low cost to the economy. Mask mandates were introduced in Canada at a time when other policy measures were relaxed, as part of the economy’s re-opening. We find that reduced restrictions on businesses or gatherings are positively associated with subsequent COVID-19 case growth – a factor that can offset and obscure the public health benefits of mask mandates. Past case totals are also found to matter for subsequent COVID-19 outcomes, suggesting that riskier behaviour may follow prior information perceived as favourable. This effect may limit how low the number of new cases can be pushed by mask mandates or other restrictions – short of a lockdown. Importantly, the estimated effect of mask mandates is relative to their absence and not absolute; depending on the situation, a mask mandate may not be sufficient to prevent an increase in new infections on its own and should be considered in conjunction with other policy measures. We have deliberately refrained from studying direct economic impacts of COVID-19, focusing instead on the unique features of the Canadian data for identifying the effect of NPIs, in particular mask mandates, on COVID-19 case growth. Future research jointly considering the epidemiological impact and the economic benefits and costs of the various public policies and restrictions would enrich the ongoing debate and provide further policy guidance. ## Data Availability All data used in this manuscript is based on information available in public domain, and is available via the data availability link and/or links provided in the manuscript. [https://github.com/C19-SFU-Econ/Public-Data](https://github.com/C19-SFU-Econ/Public-Data) ## Appendix A. Additional tables View this table: [Table A1:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T5) Table A1: Ontario – Correlations between the Google mobility indicators View this table: [Table A2:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T6) Table A2: Canada – Correlations between the Google mobility indicators View this table: [Table A3:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T7) Table A3: Ontario – Correlations between policies and location behaviour View this table: [Table A4:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T8) Table A4: Canada – Correlations between policies and location behaviour View this table: [Table A5:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T9) Table A5: Ontario – Robustness (treatment of zero weekly cases) View this table: [Table A6:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T10) Table A6: Ontario – Robustness (standard errors) View this table: [Table A7:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T11) Table A7: Ontario – Robustness (policy collinearity) View this table: [Table A8:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T12) Table A8: Canada - Robustness (treatment of zero weekly cases) View this table: [Table A9:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T13) Table A9: Canada – Robustness (standard errors) View this table: [Table A10:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T14) Table A10: Canada - Robustness (policy collinearity) View this table: [Table A11:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T15) Table A11: Canada – Robustness (news and weather) View this table: [Table A12:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T16) Table A12: Self-reported mask usage (“Always” or “Frequently”) – Canada View this table: [Table A13:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T17) Table A13: Self-reported precautions – Canada View this table: [Table A14:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T18) Table A14: Self-reported precautions – Canada (continued) View this table: [Table A15:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T19) Table A15: Canada – Closing vs. Re-opening sub-periods View this table: [Table A16:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T20) Table A16: Canada – Deaths growth (treatment of zero weekly deaths) View this table: [Table A17:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T21) Table A17: Canada - Correlations between the policy indicators View this table: [Table A18:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T22) Table A18: Canada – Location behaviour and policies View this table: [Table A19:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T23) Table A19: Ontario – Policies and information only View this table: [Table A20:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T24) Table A20: Canada – Policy and information only ## Appendix B. Additional Figures ![Figure B1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F9.medium.gif) [Figure B1:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F9) Figure B1: Self-reported mask usage in selected countries and Canadian provinces Notes: The data source is [YouGov](https://today.yougov.com/covid-19). The figure plots the average self-reported mask usage by week (the fraction of respondents that answered “Always” to the survey question “Thinking about the last 7 days, how often have you worn a face mask outside your home?”). Sample weights are used to compute the national and provincial averages. ![Figure B2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F10.medium.gif) [Figure B2:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F10) Figure B2: Canada - mask mandates and self-reported mask usage Notes: The data source is [YouGov](https://today.yougov.com/covid-19). The figure plots the average self-reported mask usage by week (the fraction of respondents that answered “Always” to the survey question “Thinking about the last 7 days, how often have you worn a face mask outside your home?”) in the provinces with and without mask mandates. Sample weights are used to compute the averages. ![Figure B3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F11.medium.gif) [Figure B3:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F11) Figure B3: Canada - Behaviour proxy, *B**it* Notes: The Behaviour proxy *B**it* is the average of the “retail”, “grocery and pharmacy”, and “workplaces” Google mobility indicators. Province-level 7-day moving averages are plotted. ![Figure B4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F12.medium.gif) [Figure B4:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F12) Figure B4: Ontario - different initial dates Notes: We plot the coefficient estimates on mask policy, with 95% confidence intervals, from equation (1), for different initial dates of the sample. The initial sample date in the baseline specifications reported in Table 1 is May 15. The left panel corresponds to baseline column (3) in Table 1; the right panel corresponds to column (4) in Table 1. ![Figure B5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F13.medium.gif) [Figure B5:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F13) Figure B5: Canada - different initial dates Notes: We plot the coefficient estimates on mask policy, with 95% confidence intervals, in the upper panel and the estimates on business/gathering policy in the lower panel, from equation (1) for different initial dates of the sample. The initial date in our baseline specification (Table 2) is March 11. The left panels correspond to column (3) in Table 2; the right panels correspond to column (4) in Table 2. ![Figure B6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F14.medium.gif) [Figure B6:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F14) Figure B6: Ontario - different lags Notes: We plot the coefficient estimates on mask policy, with 95% confidence intervals, in equation (1) for different lag values. The lag used in the baseline specifications in Table 1 is 14 days. The left panel corresponds to column (3) in Table 1; the right panel corresponds to column (4) in Table 1. ![Figure B7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F15.medium.gif) [Figure B7:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F15) Figure B7: Canada - different lags Notes: We plot the estimates on mask policy in the upper panel and the business/gathering policy in the lower panel, in equation (1) for different lag values. The lag in our baseline specification (Table 2) is 14 days. The left panels correspond to column (3) in Table 2; the right panels correspond to column (4) in Table 2. ![Figure B8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F16.medium.gif) [Figure B8:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F16) Figure B8: Canada – Weekly cases, deaths and tests (growth rate) ![Figure B9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F17.medium.gif) [Figure B9:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F17) Figure B9: Canada – Weekly cases, deaths and tests (level) ![Figure B10:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F18.medium.gif) [Figure B10:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F18) Figure B10: Canada – Daily cases, deaths and tests ![Figure B11:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F19.medium.gif) [Figure B11:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F19) Figure B11: Counterfactuals – Mask mandates in Alberta and British Columbia Notes: The figure assumes mask adoption on July 15 for two provinces that have not yet adopted mask mandates, specifically British Columbia (BC, left panel) and Alberta (AB, right panel). The counterfactual uses the mask mandate estimate −0.613 from Table 2, column (4). The counterfactual mean value (the solid line) and confidence bands (shaded area) are displayed. ![Figure B12:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/16/2020.09.24.20201178/F20.medium.gif) [Figure B12:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/F20) Figure B12: Canada - mask mandates and weekly deaths growth Notes: Average log weekly death growth in provinces with vs. without mask mandates 28 days prior. ## Appendix C. Definitions and data sources View this table: [Table C1:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T25) Table C1: Policy indicators and aggregates View this table: [Table C2:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T26) Table C2: Ontario public health regions and date of mask mandate View this table: [Table C3:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T27) Table C3: Canada COVID-19 official data sources ### Weather We downloaded historical weather data for the largest city in each province from [Weather Canada](https://weather.gc.ca/canada_e.html). The data provide daily information on 11 variables: maximum temperature (C), minimum temperature (C), mean temperature (C), heating degree-days, cooling degree-days, total rain (mm), total snow (cm), total precipitation (mm), snow on the ground (cm), direction of maximum wind gust (tens of degrees), and speed of maximum wind gust (km/h). We only use the temperature and precipitation data in Table A11 as possible factors determining outside vs. inside activity. ### News We collected data from *ProQuest Canadian Newsstream*, a subscription service to all major and small-market daily or weekly Canadian news sources. We recorded the number of search results for each day from Feb 1, 2020 to Aug 20, 2020 by searching the database for the keywords “Coronavirus” or “COVID-19”. We only counted the results with source listed as “newspaper” since other sources, such as blogs or podcasts, tend to duplicate the same original content. View this table: [Table C4:](http://medrxiv.org/content/early/2020/10/16/2020.09.24.20201178/T28) Table C4: YouGov survey questions All data used in the paper are available at [github.com/C19-SFU-Econ/data](http://github.com/C19-SFU-Econ/data). ## Appendix D. Lags Determination As discussed in Section 3.1, we assume a lag of 14 days between a change in policy or behaviour and its hypothesized effect on weekly case growth, and a lag of 28 days between such a change and its effect on weekly death growth. First, we consider the lag between infection and a case being reported. As most identified cases of COVID-19 in Canada are symptomatic, we focus on symptomatic individuals. For most provinces cases are listed according to the date of report to public health. In provinces where the dates instead refer to the public announcement, we shifted them back by one day, as announcements typically contain the cases reported to public health on the previous day. The relevant lag therefore has two components: 1. Incubation period: Most studies suggest an average incubation period of 5-6 days (e.g. 5.2 days in Li et al. (2020), 5.5 days in Lauer et al. (2020), 5.6 days in Linton et al. (2020), 6.4 days in Backer et al. (2020)). 2. Time between symptoms onset and reporting of the case to public health: the Ontario data contain an estimate of the symptom onset date (“episode date”) for each case. For our sample period the average difference between the date of report and the episode date is 4.8 days (median: 4 days) including only values from 1 to 14 days, and 6.3 days (median: 5 days) including only values from 2 to 28 days. We assume that the lags in Ontario and in other provinces are similar, and use a value of 5-6 days between symptom onset and report to public health authorities. Adding these together implies that the typical lag between infection and a positive case being reported to public health is around 11 days. Second, we consider the effect of weekly averaging on the appropriate lag for our analysis. Suppose a policy or behavioural change starts on date *t*, impacting the daily growth in infections between dates *t* − 1 and *t* and in each subsequent day. Then, assuming a lag of 11 days between infection and case reporting, case counts *C* are affected from date *t*+11 onward. Our outcome variable Δ log(ΔC) thus would react to the original policy or behavioural change on date *t* + 11. The change is complete on *t* + 23, when the week from *t* + 17 to *t* + 23 is compared to the week from *t* + 10 to *t* + 16. The midpoint of the change is *t* + 17. Choosing a lag of *l* days implies that the policy/behaviour variable phases in from *t* + *l* to *t* + *l* + 6. To match the midpoint of this phase-in to the midpoint of the change in the outcome variable, we set *l* = 14. The chosen lag matches the lag used by other authors who study COVID-19 policy interventions, e.g., CKS (2020). We explore sensitivity to alternative lags in Section 4.3. With respect to deaths, our data are, in most cases, backdated (revised by the authorities ex-post) to the actual date of death. The medical literature suggests that the mean time from symptom onset to death is around 19 days (20 days in Wu et al. (2020), 17.8 days in Verity et al. (2020), 20.2 days when accounting for right truncation in Linton et al. (2020), 16.1 days in Sanche et al. (2020), etc.), that is, two weeks longer than our estimate of the time from symptom onset to reporting of a positive test result. We accordingly set the lag used in our analysis of the death growth rate (Section 4.6) to 28 days. ## Footnotes * † NBER * * We thank Hiro Kasahara and Kevin Schnepel for excellent comments and suggestions. * Improvements to clarity and exposition throughout the text. Some references were also added. * 1 We show mask usage for the U.S. and Germany because of the related work on mask mandates by Chernozhukov et al. (2020) and Mitze et al. (2020), respectively. We show Australia as an example of a country which did not mandate mask usage (except Melbourne in late July). See Hatzius et al. (2020) for more cross-country comparisons of mask usage. * 2 Hatzius et al. (2020) estimate that a national mask mandate in the USA could replace alternative restrictions costing 5% of GDP. * 3 All data are available for download at [github.com/C19-SFU-Econ/data](http://github.com/C19-SFU-Econ/data). The COVID-19 cases, deaths and tests data that we collected and use incorporate all official ex-post revisions as of mid-August, unlike data from the Government of Canada COVID-19 website or other aggregator websites. * 4 Howard et al. (2020), a comprehensive review of the medical literature, stresses that “no randomized controlled trial (RCT) on the use of masks as source control for SARS-CoV-2 has been published.” It is unlikely that an RCT on masks’ effectiveness against COVID-19 will be feasible or ethical during the pandemic. * 5 Lyu and Wehby (2020) provide suggestive evidence that community mask mandates are more effective than employees-only mandates. * 6 The authors report average daily case growth rate of 18% in countries with no pre-existing mask norms vs. 10% in countries with such norms (S. Korea, Japan, Hong Kong and Taiwan). On a weekly basis, this is a difference of 49 log points (100(log(1.187) − log(1.17))) in case growth. * 7 Several papers study the effect of COVID-19 policies elsewhere or around the world. See, e.g., Hsiang et al. (2020), Dergiades et al. (2020), Abouk and Heydari (2020). * 8 See also Armstrong et al. (2020), who find strong evidence of policy measures on subsequent mobility behaviour in 75 US and Canadian cities. * 9 The provinces differ in the ease of accessibility of their official time series of COVID-19 cases, deaths and test numbers. In some cases, we located and used the hidden json sources feeding the public dashboard charts. In few instances in which data were not available, we used the numbers reported in the daily provincial government announcements. All COVID-19 outcome data sources are referenced and web-linked in Appendix Table C3. * 10 Additional survey data on mask usage is described and used in Section 4.4. * 11 The reports are available for download at [www.google.com/covid19/mobility/](http://www.google.com/covid19/mobility/). * 12 Each of these divisions is either entirely (in most cases) or predominantly located within a single PHU. In cases where a PHU corresponds to multiple divisions, 2016 Census population counts were used as weights to compute the PHU’s mobility index. * 13 Our chosen lag for deaths is one week longer than that used in CKS (2020). The difference is due to additional evidence from the medical literature and the construction of weekly variables (see Appendix D). * 14 We also report results using the growth rate of deaths as supplementary analysis in Section 4.6. * 15 Figures B9 and B10 in the Appendix respectively display the weekly and daily cases, deaths and tests in each Canadian province over time. There are markedly lower numbers reported on weekends or holidays. * 16 To deal with zero weekly values, which mostly occur in the smaller regions, as in CKS (2020), we replace log(0) with −1. We also check the robustness of our results by adding 1 to all Δ*C**it* observations before taking logs, by replacing log(0) with 0, and by using population weighted least squares; see Tables A5 and A8. * 17 There was no PHU-wide mask mandate in Lambton as of August 31, but its main city, Sarnia, enacted a mask mandate on July 31. * 18 We do not use provincial declarations of emergency in our analysis as they are mostly legal tools enabling other restrictions rather than restrictions *per se*. * 19 The daily numerical values for each of the 17 basic policy indicators and the 5 policy aggregates for each province and date are available on the project’s [Github repository](https://github.com/C19-SFU-Econ/data). * 20 Alternative methods of computing the standard errors are explored in Section 4.3. * 21 We drop the “transit”, “parks”, and “residential” location indicators because, respectively, 10.6%, 13.7%, and 2.8% of the observations are missing in the provincial data, and 20.7%, 52.1%, and 11.1% are missing in the Ontario data. The “transit” and “residential” variables are also highly correlated with the three indicators included in our location behaviour proxy *B**it*. Furthermore, the “parks” indicator does not have clear implication for COVID-19 outcomes. * 22 In the Ontario analysis, 1.4% of the *B**it* values were imputed via linear interpolation. * 23 Mask mandates and regulations on businesses and gatherings vary at the PHU level. Long-term care policy changed only province-wide. The school and travel policies do not vary in the sample period and hence are omitted from the analysis with Ontario PHU data. * 24 Table A6 in the Appendix reports alternative standard error specifications: regular clustering at the PHU level (Stata command “cluster”), wild bootstrap standard errors clustered at the PHU level, and wild bootstrap standard errors two-way clustered by PHU and date. Our results are robust to alternative ways of computing standard errors. * 25 In all tables, *Variable*_14 denotes the 14-day lag of *Variable*. * 26 Using equations (1) and (3), a coefficient estimate ![Graphic][12] on Mask 14, as the latter changes from 0 to 1, corresponds to a ![Graphic][13] percent change in the ratio of current-week to past-week cases, Δ*C**it*/Δ*C**it*−7. For example, suppose weekly cases in weeks 1, 2 and 3 are 100, 150 and 225 respectively, showing a 50% weekly case growth. Our estimates suggest that, all else equal, a mask mandate implemented at the beginning of week 1 would reduce week 3 cases to 155–169 (instead of 225 in absence of mandate), with continued reductions in the next several weeks. * 27 We also tried including each location change measure separately and the results were similar (not shown). * 28 Figure 4 assumes a July 7 mask mandate implementation date for Ontario (when its most populous PHU, Toronto, adopted a mask mandate, along with Ottawa), and July 18 for Quebec (province-wide mandate). * 29 Sensitivity analysis with the other specifications is available upon request. * 30 Table A4 shows a correlation of 0.61 between the Travel and School policy aggregates. * 31 Aggregating the 17 basic policy indicators into five groups mitigates this issue. Here, we test whether any remaining collinearity poses a problem. * 32 535 out of 3,094 observations (17%) in Table 1 and 230 out of 1,560 observations (15%) in Table 2 are affected, mostly in the small provinces or PHUs. When both Δ*C**it* and Δ*C**it*−7 equal zero, the weekly case growth rate is *Y**it* = 0. * 33 Vancouver, BC; Calgary, AB; Saskatoon, SK; Winnipeg, MB; Toronto, ON; Montreal, QC; Moncton, NB; Halifax, NS; Charlottetown, PE; and St. John’s, NL. * 34 The YouGov data is available at: [https://yougov.co.uk/covid-19](https://yougov.co.uk/covid-19). * 35 As on Figure 4, we use July 7 as the mask mandate implementation date for Ontario. * 36 Since mask usage is reported only for specific dates within each survey wave, we use the mask policy variable daily values for these same dates instead of the weekly moving average. * 37 Similarly, in Table A12, column (4) shows that “Always” or “Frequent” mask usage increases by 21.5 percentage points. The finding that the increase in mask usage among the “Always” respondents is larger than among the “Always” or “Frequent” respondents is consistent with some people switching from wearing masks “frequently” to “always”. * 38 Hatzius et al. (2020) document that state mask mandates in the US increased mask usage roughly by 25 percentage points in 30 days. The compliance with mask mandates may differ across countries or regions based on social norms, peer effects, political reasons or the consequences of noncompliance (e.g., fines). * 39 If we take the increase of about 30 percentage points in reported mask usage induced by mask mandates at face value, the full effect of mask wearing (treatment-on-the-treated effect) would be roughly triple our estimates. It could be larger still if there is desirability bias in answering the mask usage survey question, so that the actual increase in mask use may be smaller than our estimate. * 40 Consistent with this result, Seres et al. (2020) find that wearing masks increased physical distancing based on a randomized field experiment in German stores. * 41 June 12 is the date of the earliest mask mandate in Ontario. For the national analysis, July 7, the effective date for Toronto and Ottawa, is considered Ontario’s first significant date of mask mandate enactment: PHUs with earlier mandates account for less than 10% of Ontario’s population. * 42 In Table 4, *Variable*_28 denotes the *Variable* lagged by 28 days. * 43 205 out of the 1,470 observations (14%) had log(0) replaced by −1. * Received September 24, 2020. * Revision received October 16, 2020. * Accepted October 16, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. [1].Abaluck, J., J. Chevalier, N. Christakis, H. Forman, E. Kaplan, A. Ko and S. Vermund (2020), “The Case for Universal Cloth Mask Adoption and Policies to Increase Supply of Medical Masks for Health Workers”, COVID Economics (5), CEPR Press. 2. [2].Abouk, R. and B. Heydari (2020), “The Immediate Effect of COVID-19 Policies on Social Distancing Behavior in the United States”, working paper. 3. [3].Armstrong, D., M. Lebo and J. Lucas (2020), “Do COVID-19 Policies Affect Mobility Behaviour? Evidence from 75 Canadian and American Cities”, Canadian Public Policy 46(S2): S127-S144. 4. [4].Backer, J., D. Klinkenberg, J. Wallinga (2020), “Incubation period of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China, 20-28 January 2020”, Euro Surveillance 25(5): pii=2000062. 5. [5].Cameron, A., J. Gelbach and D. Miller (2008), “Bootstrap-based improvements for inference with clustered errors”, Review of Economics and Statistics 90: 414–427 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1162/rest.90.3.414&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258260500003&link_type=ISI) 6. [6].Chernozhukov, V., H. Kasahara and P. Schrimpf (2020), “Causal Impact of Masks, Policies, Behavior on Early COVID-19 Pandemic in the U.S.”, medRxiv pre-print. 7. [7].Chu, D., E. Akl, S. Duda, K. Solo, S. Yaacoub and H. Schünemann (2020), “Physical distancing, face masks, and eye protection to prevent person-to-person transmission of SARS-CoV-2 and COVID-19: a systematic review and meta-analysis”, Lancet 395: 1973–87. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)31142-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F16%2F2020.09.24.20201178.atom) 8. [8].Dergiades, T., C. Milas, E. Mossialos and T. Panagiotidis (2020), “Effectiveness of Government Policies in Response to the COVID-19 Outbreak”, working paper. 9. [9].Greenhalgh, T., M.B. Schmid, T. Czypionka, D. Bassler and L. Gruer (2020), “Face masks for the public during the covid-19 crisis”, BMJ 369: m1435. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjkvYXByMDlfMTEvbTE0MzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMC8xNi8yMDIwLjA5LjI0LjIwMjAxMTc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 10. [10].Hatzius, J., D. Struyven and I. Rosenberg (2020), “Face Masks and GDP” (2020), Goldman Sachs Research, [https://www.goldmansachs.com/insights/pages/face-masks-and-gdp.html](https://www.goldmansachs.com/insights/pages/face-masks-and-gdp.html). 11. [11].Howard, J., A. Huang, Z. Li, Z. Tufekci, V. Zdimal, H.-M. v.d. Westhuizen, A. v. Delft, A. Price, L. Fridman, L-H. Tang, V. Tang, G. L. Watson, C.E. Bax, R. Shaikh, F. Questier, D. Hernandez, L.F. Chu, C.M. Ramirez and A.W. Rimoin (2020), “Face Masks Against COVID-19: An Evidence Review”, pre-print. 12. [12].Hsiang, S., D. Allen, S. Annan-Phan, K. Bell, I. Bolliger, T. Chong, H. Druckenmiller, L. Huang, A. Hultgren, E. Krasovich, P. Lau, J. Lee, E. Rolf, J. Tseng and T. Wu (2020), “The effect of large-scale anti-contagion policies on the COVID-19 pandemic”, Nature 584: 262–267. 13. [13].Greenhalgh, T., M. Schmid, T. Czypionka, D. Bassler and L. Gruer (2020), “Face masks for the public during the covid-19 crisis”, BMJ 369:m1435 doi: 10.1136/bmj.m1435. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjkvYXByMDlfMTEvbTE0MzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMC8xNi8yMDIwLjA5LjI0LjIwMjAxMTc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 14. [14].Jones, S., Imperial College London Big Data Analytical Unit and YouGov Plc. (2020), Imperial College London YouGov Covid Data Hub, v1.0, YouGov Plc, retrieved September 2020. 15. [15].Kermack, W. and A. McKendrick (1927), “A contribution to the mathematical theory of epidemics”, Proceedings of the Royal Society A 115(772): 700–721. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1927.0118&link_type=DOI) 16. [16].Lauer, S., K. Grantz, Q. Bi, F. Jones, Q. Zheng, H. Meredith, A. Azman, N. Reich and J. Lessler (2020), “The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application”, Annals of Internal Medicine 172: 577–582. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-0504&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32150748&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F16%2F2020.09.24.20201178.atom) 17. [17].Leung, N., D. Chu, E. Shiu, K. Chan, J. McDevitt, B. Hau, H. Yen, Y. Li, D. Ip, J. Malik Peiris, W. Seto, G. Leung, D. Milton and B. Cowling (2020), “Respiratory virus shedding in exhaled breath and efficacy of face masks”, Nature Medicine 26: 676–680. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0843-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32371934&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F16%2F2020.09.24.20201178.atom) 18. [18].Li, Q., X. Guan, P. Wu, X. Wang, L. Zhou, Y. Tong, R. Ren, K. Leung, E. Lau, J. Wong, X. Xing, N. Xiang, Y. Wu, C. Li, Q. Chen, D. Li, T. Liu, J. Zhao, M. Liu, W. Tu, C. Chen, L. Jin, R. Yang, Q. Wang, S. Zhou, R. Wang, H. Liu, Y. Luo, Y. Liu, G. Shao, H. Li, Z. Tao, Y. Yang, Z. Deng, B. Lui, Z. Ma, Y. Zhang, G. Shi, T. Lam, J. Wu, G. Gao, B. Cowling, B. Yang, G. Leung and Z. Feng (2020), “Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia”, The New England Journal of Medicine 382(13): 1199–1207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F16%2F2020.09.24.20201178.atom) 19. [19].Linton, N., T. Kobayashi, Y. Yang, K. Hayashi, A. Akhmetzhanov, S. Jung, B. Yuan, R. Kinoshita and H. Nishiura (2020), “Incubation Period and Other Epidemiological Characteristics of 2019 Novel Coronavirus Infections with Right Truncation: A Statistical Analysis of Publicly Available Case Data”, J. of Clinical Medicine 9(2): 538. 20. [20].Lyu, W. and G. Wehby (2020), “Community Use Of Face Masks And COVID-19: Evidence From A Natural Experiment Of State Mandates In The US”, Health Affairs 39(8): 1419–1425. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1377/hlthaff.2020.00818&link_type=DOI) 21. [21].Meyerowitz, E., A. Richterman, R. Gandhi and P. Sax (2020), “Transmission of SARS-CoV-2: A Review of Viral, Host, and Environmental Factors”, Annals of Internal Medicine, Reviews, 17 September 2020. 22. [22].Mitze, T., R. Kosfeld, J. Rode and K. Wälde (2020), “Face Masks Considerably Reduce COVID-19 Cases in Germany”, working paper. 23. [23].Mohammed, A., R. Johnston and C. van der Linden, (2020), “Public Responses to Policy Reversals: The Case of Mask Usage in Canada during COVID-19”, Canadian Public Policy 46(S2): S119–S126. 24. [24].Organization for Economic Co-operation and Development (OECD) (2020), “Flattening the COVID-19 peak: Containment and mitigation policies”, OECD Policy Responses to Coronavirus (COVID-19), [https://www.oecd.org/coronavirus/policy-responses/flattening-the-covid-19-peak-containment-and-mitigation-policies-e96a4226/](https://www.oecd.org/coronavirus/policy-responses/flattening-the-covid-19-peak-containment-and-mitigation-policies-e96a4226/). 25. [25].Prather, K., C. Wang and R. Schooley (2020), “Reducing transmission of SARS-CoV-2” Science 368(6498): 1422–1424. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjgvNjQ5OC8xNDIyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTAvMTYvMjAyMC4wOS4yNC4yMDIwMTE3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 26. [26].Raifman, J., K. Nocka, D. Jones, J. Bor, S. Lipson, J. Jay, P. Chan, M. Brahim, C. Hoffman, C. Corkish, E. Ferrara, E. Long, E. Baroni, F. Contador, H. Simon, M. Simko, R. Scheckman, S. Brewer, S. Kulkarni, F. Heykoop, M. Patel, A. Vidyasagaran Chiao, C. Safon and S. Burkhart (2020), “COVID-19 US state policy database”, [https://tinyurl.com/statepolicies](https://tinyurl.com/statepolicies). 27. [27].Sanche, S., Y.T. Lin, C. Xu, E. Romero-Severson, N. Hengartner and R. Ke (2020), “High Contagiousness and Rapid Spread of Severe Acute Respiratory Syndrome Coro-navirus 2”, Emerging Infectious Diseases 26(7): 1470–1477. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F16%2F2020.09.24.20201178.atom) 28. [28].Seres, G., A. Balleyer, N. Cerutti, J. Friedrichsen, and M. Sauer (2020), “Face Mask Use and Physical Distancing before and after Mandatory Masking: Evidence from Public Waiting Lines”, working paper. 29. [29]. Statistics Canada (2020), “Table 17-10-0009-01 Population estimates, quarterly”, https://doi.org/10.25318/1710000901-eng,retrieved September 2020. 30. [30].Verity, R., L. Okell, I. Dorigatti, P. Winskill, C. Whittaker, N. Imai, G. Cuomo-Dannenburg, H. Thompson, P. Walker, H. Fu, A. Dighe, J. Griffin, M. Baguelin, S. Bhatia, A. Boonyasiri, A. Cori, Z. Cucunubá, R. Fitzjohn, K. Gaythorpe, W. Green, Hamlet, W. Hinsley, D. Laydon, G. Nedjati-Gilani, S. Riley, S. v. Elsland, E. Volz, H. Wang, Y. Wang, X. Xi, C. Donnelly, A. Ghani and N. Ferguson (2020), “Estimates of the severity of coronavirus disease 2019: a model-based analysis”, Lancet Infectious Diseases 20: 669–677. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30243-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32240634&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F16%2F2020.09.24.20201178.atom) 31. [31].World Health Organization (2020), “Weekly Epidemiological Update, Coronavirus disease 2019 (COVID-19), 21 September 2020”, [https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200921-weekly-epi-update-6.pdf](https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200921-weekly-epi-update-6.pdf). 32. [32].Wu, J., K. Leung, M. Bushman, N. Kishore, R. Niehus, P. de Salazar, B. Cowling, M. Lipsitch and G. Leung (2020), “Estimating clinical severity of COVID-19 from the transmission dynamics in Wuhan, China”, Nature Medicine 26: 506–510. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0822-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F16%2F2020.09.24.20201178.atom) 33. [33].Yuksel, M., Y. Aydede and F. Begolli (2020), “Dynamics of Social Mobility during the COVID-19 Pandemic in Canada”, IZA Discussion Paper No. 13376. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/inline-graphic-1.gif [6]: /embed/inline-graphic-2.gif [7]: /embed/inline-graphic-3.gif [8]: /embed/inline-graphic-4.gif [9]: /embed/inline-graphic-5.gif [10]: /embed/inline-graphic-6.gif [11]: /embed/graphic-13.gif [12]: /embed/inline-graphic-7.gif [13]: /embed/inline-graphic-8.gif