Effect of specific non-pharmaceutical intervention policies on SARS-CoV-2 transmission in the counties of the United States =========================================================================================================================== * Bingyi Yang * Angkana T. Huang * Bernardo Garcia-Carreras * William E. Hart * Andrea Staid * Matt D.T. Hitchings * Elizabeth C. Lee * Chanelle J. Howe * Kyra H. Grantz * Amy Wesolowksi * Joseph Chadi Lemaitre * Susan Rattigan * Carlos Moreno * Brooke Borgert * Celeste Dale * Nicole Quigley * Andrew Cummings * Alizée McLorg * Kaelene LoMonaco * Sarah Schlossberg * Drew Barron-Kraus * Harrison Shrock * UFCOVID Interventions Team * Justin Lessler * Carl D. Laird * Derek A.T. Cummings ## Abstract Non-pharmaceutical interventions (NPIs) remain the only widely available tool for controlling the ongoing SARS-CoV-2 pandemic. We estimated weekly values of the effective basic reproductive number (*R**eff*) using a mechanistic metapopulation model and associated these with county-level characteristics and NPIs in the United States (US). Interventions that included school and leisure activities closure and nursing home visiting bans were all associated with an Reff below 1 when combined with either stay at home orders (median *R**eff* 0.97, 95% confidence interval (CI) 0.58-1.39)* or face masks (median *R**eff* 0.97, 95% CI 0.58-1.39)*. While direct causal effects of interventions remain unclear, our results suggest that relaxation of some NPIs will need to be counterbalanced by continuation and/or implementation of others. ## Introduction Months after the emergence of SARS-CoV-2 and its subsequent pandemic spread, widespread transmission continues. The United States (US) has been particularly affected, although the burden of disease has been geographically heterogeneous. While Northeastern states like New York, New Jersey, and Massachusetts were hardest hit when the virus first arrived in the US from March to April, 2020, states like Texas, Florida, Arizona, and California, which had previously avoided substantial disease burden, experienced rapidly growing cases between June and August (Fig. 1A) (*1*). Many factors likely contribute to spatial and temporal heterogeneity in COVID-19 incidence, including socio-demographic characteristics, the frequency of importation of infections, and the local use and timing of non-pharmaceutical interventions (NPIs) (*2*–*5*). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/03/2020.10.29.20221036/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/11/03/2020.10.29.20221036/F1) Figure 1. Confirmed COVID-19 cases and reproduction numbers (Reff) in the US. **A)** Daily confirmed COVID-19 cases and deaths. **B**) Distribution of Reff by weeks since the county’s first reported case. Grey horizontal line indicates the threshold of Reff = 1 (same as in C). Medians (points), interquartiles (dark vertical lines) and 95% percentiles (light vertical lines) are shown (same as in C). **C**) Temporal distribution of Reff stratified by county population size. **D**) Map of county Reff for representative weeks. Weeks were selected with a 3-week interval from the last week when Reff were available. Grey indicates no data available. While there is a continuing need to control SARS-CoV-2 transmission, prolonged social and economic strains push the limits of population compliance. Hence, it is critical to identify targeted and effective strategies for disease control. Variation in NPI application and timing across US counties provides an opportunity to compare the effectiveness of interventions in reducing transmission, but these effects are difficult to estimate directly due to spatial and temporal variation in transmission and an imperfect and ever-changing surveillance system. We hypothesized that by using multiple data streams and focusing on estimates of transmissibility rather than raw counts of cases or deaths, some of these challenges may be overcome. Here, we estimated the transmissibility of SARS-CoV-2 in 3,036 (out of 3,142, 97%) US counties using a mechanistic meta-population model that incorporates spatial coupling of transmission between counties to estimate weekly effective basic reproductive numbers (i.e., Reff, the reproductive number adjusted for changes due to factors other than population susceptibility, such as social distancing) from confirmed cases and deaths from January 21 to July 5, 2020. We associated these Reff estimates with NPIs and county level demographics, while accounting for, temporal variation, autocorrelation and uncertainty in our estimates. We aimed to determine which NPIs have been most effective so far to inform future implementations. ## Results A total of 2,846,249 COVID-19 cases and 128,391 deaths were reported in the US as of July 5, 2020 (Fig. 1A). Cases first appeared in coastal counties with large populations in late-January and were reported in 68% of US counties by March 31. Weekly Reff fit to confirmed cases (Fig. 1B-C) suggests that counties with larger population sizes (>90,000) experienced earlier and more efficient transmission (i.e. greater Reff) (median Reff 2.6, interquartile range (IQR) 1.5 to 6.5 in the weeks before April) (Fig. 1C-D; Fig. S1). Later in the epidemic Reff dropped in these large counties (median Reff of 0.8, IQR 0.7 to 1.0 in the week ending May 2), and though Reff was similar in small counties, some had appreciably higher transmission (median Reff of 0.9, IQR 0.6 to 1.8) (Fig 1C-D). We compiled detailed data on the timing of state level NPIs policies and county level data on school closure. Enactment of NPIs started in early March and peaked the week ending April 11, by which time 100% of counties had closed public schools; and, of states, 98% had closed leisure activities (restaurants, gyms and movie theaters), 88% had stay at home orders, 63% had suspended medical services, 59% had banned nursing home visits and 29% had closed daycares (Fig. 2A; Fig. S2). In most counties, the majority of interventions were implemented before a county had its first case; with school closure coming earliest (median 1.4 weeks before first case IQR 2.6 to 0.6 weeks) (Fig. 2B). The one exception is face masking orders which were initiated on average 5.7 weeks (IQR 4.1 to 7.0 weeks) after the first case (Fig. 2B). Many locations started to gradually lift control measures in mid-April, particularly medical service suspensions (remained in only one state as of June 13 (2%)), stay at home orders (12%) and leisure activities closures (45%) (Fig. 2A, S2). At time of analysis no county had lifted school closure. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/03/2020.10.29.20221036/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/11/03/2020.10.29.20221036/F2) Figure 2. Temporal distributions of non-pharmaceutical interventions. **A**) Time series of the proportion of states and counties where interventions were implemented. **B)**Distribution of time differences between intervention times and the occurrence of the county’s first confirmed case. Colored, solid lines indicate the median difference times of each intervention. In March, the average Reff in each state was consistently high across the country (median 4.8, IQR 4.1 to 5.3), but transmission had reduced sharply by May (median 1.2, IQR 0.9 to 1.5). However, these results mask substantial variation within states at all points in the epidemic (mean intra-state coefficient of variation (CV) 1.39, IQR 1.28-1.44 in March, and 1.16, IQR 0.78-1.66 in May). Using generalized estimating equations (GEEs) to account for autocorrelation, we estimated the county-level effects of intervention policies on log-transformed Reff adjusting for log-population size, proportion of individuals in poverty, median household income and other county level covariates (Tables S1-3; see Methods for description of all models considered). Transmission was consistently higher in counties with greater population sizes (21% increase in Reff per 1,000 increase in population size, 95% CI 13 to 29%) and those with a higher proportion of people without college educations (5% per 10% increase, 95% CI 3 to 6%), while transmission was consistently lower in counties with higher proportions of white individuals (2.6% per 10% decrease in the proportion, 95% CI 2.0 to 3.2%) and those with a lower median age (0.6% per 1 year decrease in median age, 95% CI 0.3 to 0.9%). We found that transmission had the strongest association with school closure (37% reduction in Reff, 95% CI 33 to 40%), followed by daycare closure (31%, 95% CI 26 to 35%) and banning nursing home visits (26%, 95% CI 23 to 29%) (Fig. 3A; main model in Tables S1-2). Stay at home orders were associated with a 15% reduction in Reff (95% CI 13 to 17%) while face-mask orders were associated with a 18% reduction (95% CI 16% to 20%). Log-linear regression models that included a lag-term to account for autocorrelation (termed the lag model, an alternative approach to GEE) yielded similar regression coefficients (Fig. 3A and Table S4). To ensure the contribution of interventions in reducing Reff when prioritizing parsimony, we fitted LASSO regression on the lag model model across a range of penalties for model complexity and found that school closure was the intervention most consistently associated with reductions in Reff (Fig. 3B). Including information on NPIs improved model performance across multiple metrics (Tables S2, S3) (see Methods). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/03/2020.10.29.20221036/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/11/03/2020.10.29.20221036/F3) Figure 3. Effects of non-pharmaceutical interventions (NPIs) and county-level characteristics on transmission. **A**) Effects of NPIs and county-level characteristics estimated from the main model. Models were fitted with both GEEs and lag model regression. The order in y-axis (same for C) is according to the importance of covariates in explaining the variances shown in B. **B**) The importance of covariates in explaining the variances. Main models that were formulated for lag model regression were fitted to LASSO with increasing parsimony. **C**) Changes in the estimated effects when each covariate is dropped in the main model fitted with lag model regression. Implementation of interventions was highly temporally correlated both between and within counties (Figs. 2, S2 and S3), presenting challenges to estimating independent associations. For example, school closures occurred at the same time in much of the country, with closure only reported in three unique weeks making separation of the effect of school closure from the effect of week of year impossible (including week of year as a categorical covariate eliminates the observed effect Fig. S4). However, whether causation, coincidence or confounding, reduction in Reff occurred for some reason; hence, we performed several analyses to better disentangle the observed associations with NPIs. Univariate and multivariate models for each NPI showed similar relative associations (Fig. S5) as did models that included an effect for time since the first case in a county (Fig. S6). To assess the impact of collinearities, we reran the main model holding out each NPI one at a time. The largest changes in estimated coefficients were seen when dropping school closure, which substantially impacted the estimated association of Reff with remaining coefficients (e.g. the coefficient of banning nursing home visits indicated a 6.4% larger reduction in models without school closure) (Fig. 3C). To further ensure the observed associations between NPIs and reductions in Reff were not merely the result of spurious associations between changes in Reff and changes in NPI, we permuted data on NPIs across counties (Fig S7A,B) and within counties (Fig. S7C,D), respectively and found all associations to be substantially dampened or eliminated. People changed their behavior in response to the pandemic, whether due to policy or personal choice. An important behavior that impacts SARS-CoV-2 transmission is travel to and social contact in different settings. Data on workplace presence relative to pre-pandemic periods were available on Google users (*6*) and, unlike other measures of movement in these data, were available for the nearly all (98%) of the county-weeks in our analysis. Across the US, large changes in workplace presentation occurred in March (Figure S8). We explored potential mediation/confounding by workplace presence in each county. We found that, at least in part, the relationship with a number of NPIs including school closure (32% [95% CI 28%, 34%] of total effect explained), leisure activities (65% [95% CI 59%, 71%] and stay at home orders (100% [95% CI 94%, 109%]) were mediated and/or confounded by workplace presence (Fig. S9, S10). No single intervention was implemented alone for a sustained period of time in the period of our study, and many combinations of interventions never appear (e.g., there are no cases where stay at home orders are in effect but school closure is not). Hence, we only observe the impact of suites of interventions that may have complex interactions. We fit a GEE model to each of the unique suites of NPI utilized by counties (including no intervention) as categorical variables (Fig 4). We also fit a boosted regression tree model that can account for complex interactions between interventions (see Methods). Estimates from these models were consistent with those calculated from our main model. Interventions that included school and leisure activity closure and nursing home visiting bans were all associated with an Reff below one when combined with either stay at home orders (median *R**eff* 0.97, 95% CI 0.58-1.39) or face masks (median *R**eff* 0.97, 95% CI 0.58-1.39) (Fig. 4). Inclusion of more interventions further reduced Reff, with a minimum median *R**eff* of 0.65 (0.39, 1.08) when all interventions are in place (Fig. 4, Table S7). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/03/2020.10.29.20221036/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/11/03/2020.10.29.20221036/F4) Figure 4. Prediction of reproduction numbers (Reff) associated with different non-pharmaceutical interventions (NPIs) suites. **A**) Non-pharmaceutical interventions (NPIs) suites: no intervention (left column), all interventions (right column) and the ten most frequently used NPI suites in counties in the US (middle columns; ordered by number of NPIs in each suite). Numbers below the columns are percentages of county-weeks in the study in which the suites were active. **B**) Reff with different NPI suites. Median (cross), interquartile range (boxes) and 95% quantile range (vertical lines) of Reff were distributions of predictions from the main model fitted with GEE (red, Fig. 3A), GEE with NPIs fitted as suites (brown), and boosted decision tree (purple). Horizontal dashed line indicates the threshold of one. Multiple sensitivity analyses found broadly similar estimates of associations of NPI with Reff, these included: subsetting our analysis to exclude highly uncertain Reff estimates from early in the pandemic (Fig. S11); allowing interventions to have an ongoing impact after they are lifted (Fig. S12); using alternative reconstruction methods (Fig. S13); and using estimates derived from reported deaths (Figs. S14). In the last analysis, we did observe reduced effects of most NPIs and a larger effect of school closure; though results were otherwise qualitatively the same. ## Discussion A strength of our approach is the detailed data on NPIs that we compiled as well as the robust estimation of transmissibility using a mechanistic meta-population model. We found that six of seven NPIs examined were associated with reductions in Reff, with school closure being associated with the greatest reduction. However, we can only speculate about the mechanistic pathways by which any of these policies may have caused reductions in transmission. School closures, for example, may have had substantial impacts on the social interactions of non-school-aged individuals as parents and workplaces adapt to accommodate changes in children’s schedules, as suggested by our mobility analysis. It may also be a leading indicator of community attitudes about transmission. Our estimates indicate that four interventions were necessary, on average, to reduce Reff below one, and that even with seven interventions, reproductive numbers remained above 0.65 on average (Fig 4). However, we note that our estimates do not include the effect of immunity due to substantial infection in many areas of the US, which will cause additional reductions in R. Overall, even when the greatest number of intervention policies were in place (only observed in 1.2% of county-weeks; Fig. 4A), we never saw reductions as large as those seen in Asia and parts of Europe, where reproductive rates fell as low as 0.44 (*7*). Whether this is the result of poor compliance, structural factors, or states easing restrictions before they could have their full effect remains unknown. However, even though reductions were not as large as those seen in other areas, stark and immediate reductions in reproduction numbers across the US coincided with the use of NPI. We are not the first to attempt to estimate the impact of NPIs on transmission of SARS-CoV-2. Previous analyses agree with ours in several important dimensions, including the clear association of business closures, stay at home orders, and masking with significant reductions in transmission (*8, 9*). Beyond the impact of NPIs, other studies identified population density as a risk factor for transmission (*10*) while we found that county population size had a stronger association (Table S6). There remains debate on the impact of school closure. Consistent with our results, Auger et al. (*2*) found significant associations between school closure and reduced incidence. However, others, e.g., Hsiang et al. 2020 (*8*), have found no association between school closure and the rate of growth in cases. These discrepancies may in part be explained by differences in how investigators treat the temporal association of interventions with changes in observed epidemic patterns (i.e., our estimates are lagged based on modeled epidemic processes and reporting delays), but the fact is that we never see an instance in the data where stay at home orders and many other effective interventions are implemented without school closure, so all other effects are measured in the presence of school closure. How this is handled in analysis can lead to substantial differences in estimated effect. The impact of school closures need not necessarily stem directly from reductions in direct infection by school aged individuals. The associated reductions in Reff we observed suggest substantial confounding or indirect effects of school closure. The association between school closure and time spent in the workplace provides one possible mechanism of indirect effects, however others exist. Further, even if transmission is rare schools bridge many communities, and can play an important role in facilitating epidemic spread by connecting these subpopulations. Regardless of why, the data and past experience show important associations between school closure and transmission which should not be dismissed when setting policy. Aside from the correlation in timing of interventions, other factors may also challenge our inferences. Compliance with policies and lags between implementation and actions by individuals could obscure the impact of policies on transmission. School closure is unique in this regard among the NPI we considered as schools stayed closed with effectively 100% compliance over the period of our analysis. There may be confounders that were unmeasured or not included in our model. Though we were able to obtain county level data on school closures, we were limited to state level data on other interventions. Local changes that affected large populations within a state may lead to misclassification bias. Large spatial and temporal variation in the accuracy of surveillance for confirmed cases or deaths could induce spurious changes in Reff that do not reflect true transmission (*11*), however using both cases and deaths for our inference helps mitigate this possibility. Many counties reported limited numbers of cases and/or deaths and thus infection dynamics could not be reconstructed. We assumed a stable distribution of delays between infection and the time of confirmation or death, though this could have varied over the course of the outbreak (*12*). Not all behavioral change is captured by concrete policies; however, our focus was on the impact of policies enacted by governments rather than actions taken by individuals in the absence of such policies. Despite these limitations our analysis provides critical insight into how individual interventions, or at least commonly used suites of interventions, may affect the spread of SARS-CoV2. These estimates are critical as governments attempt to figure out how to respond to resurgent cases and look for responses that successfully control spread while allowing as much of normal economic and social life to continue as possible. The fact that multiple NPIs were needed to reduce Reff below one, suggests that relaxation of some NPIs might need to be counterbalanced by continuation and/or implementation of others. Our point estimates of the relative contribution of each intervention provide some guidance in making these difficult decisions. ## Materials and Methods Data and code for analysis are made available at [https://github.com/datcummings/SARS-CoV-2.NPI.analysis](https://github.com/datcummings/SARS-CoV-2.NPI.analysis). ### Data on COVID-19 cases and deaths Laboratory-confirmed COVID-19 cases and deaths in the counties of the United States were retrieved from USAFacts ([https://usafacts.org/visualizations/coronavirus-covid-19-spread-map/](https://usafacts.org/visualizations/coronavirus-covid-19-spread-map/)) on July 5, 2020. The data from USA Facts was used for all counties except for New Mexico, where we observed timing offsets in weekend reporting for a few, select dates and counties. Data from the NY Times did not have these issues and thus were used for New Mexico (1). ### Data on county-level demographic characters We obtained data on county-level demographic characteristics (i.e. population size, population median age, land area, median annual income, number in poverty, number reported as white in race, and number with education below college) from the 2014-2018 American Community Survey, which were retrieved through *tidycensus* package version 0.9.9.2 (*13*) in R. Population density was calculated by dividing the population size over land area. Proportions of population in poverty, white in race, and education below college were derived by dividing the numbers with the total number of people surveyed for those characteristics. #### Data on non-pharmaceutical interventions (NPIs) We obtained dates of policy announcements of closures of public schools (K-12) of each county by consulting the government websites of school district, county and state and local news sources. We used the earliest documented date as the county’s school closure date when there were multiple dates available for districts within that county. When closures were announced during or at the end of a planned school break, the date when schools were last in session was reported. We found school closure dates for 94.3% (2,963 out of 3,142) of counties. We limited our analysis to the 2,963 counties for which school closure information was available (from the 3,036 for which we had estimates of Rt). We obtained dates of implementation and termination of the other NPIs at state-level from COVID-19 US State Policy Database (CUSP) on July 6, 2020 (*14*). Eleven interventions that were directly used to reduce transmissions were extracted from the dataset (Table S5). Interventions were grouped if they were semantically similar and clustered in time (e.g. face mask mandated in the public and face mask mandated in businesses; Figure S11). When multiple dates of interventions and terminations were available for the grouped NPIs, we used the earliest documented date for implementations and the latest date for termination. An intervention is considered as implemented if the date of the corresponding Reff was between its implementation and termination date. ### Estimating reproduction number We fit a mechanistic transmission model to confirmed cases of COVID-19 in each of the counties of the US. Separately, we also fit these same models to deaths due to COVID-19 in each county. We used spatially coupled SEIR models to represent the transmission of SARS-CoV-2 with separate metapopulations representing the population in each US county. We estimated the incidence of infection in each county based on confirmed count data (and separately deaths due to COVID-19) using the back-projection method of Becker et al (*15*). For our deconvolution procedure, numbers of cases and deaths were upscaled by sampling from a negative binomial (with probability 1/10 and 1/200 respectively). Second, using the delay distributions for cases and deaths, we applied the function ‘backprojNP’ from package ‘surveillance’ in R to back-project incidence. Finally, to account for right-truncation of incidence values, we followed Abbott et al (*16*) and sampled the estimated incidence to account for infections that had occurred but had not yet been reported or confirmed from a negative binomial distribution, where the probability was given by the respective cumulative delay distributions. The most recent 14 days were then removed. We constructed 100 stochastic realizations of this algorithm for each county. Using the resultant time series of daily infections in each county, we fit our mechanistic model with state variables S (susceptible to SARS-CoV-2), E (exposed, infected but not yet infectious), I (infectious), and R (recovered or deceased and no longer infectious) to each stochastic realization using non-linear optimization tools (IPOPT, Pyomo)(*17, 18*). Models were developed using epi_inference, a software package to be released. Our model included three infectious components to allow the distribution of durations of infectiousness to approximate a gamma distribution. Fitted parameters were piecewise constant transmission coefficients in each county for each week. Estimates used two weeks of incident infections to derive each piecewise constant transmission coefficient in each week (with estimates assigned to the first week of the two weeks used to estimate infections). A mobility matrix derived from US Census commuting data from pre-pandemic time periods was used to specify the fraction of commuters in each county and the fraction of time that those commuters spent in counties other than the ones they reside. ### Natural History of SARS-CoV-2 Time delays from infection to confirmation among those cases that are confirmed was assumed to be log-normally distributed with mean of days (log(8)=2.07 days). This assumption was derived from (*19*) and assumed confirmation comes on average 1.5 days after attendance at a medical facility. We assumed a log-standard deviation of 0.3 of this distribution (*20*). Time delays from infection to death among those that die was assumed to be log-normally distributed with log-mean 2.84 and log-standard deviation of 0.72 and assuming no competing risks (*21, 22*). ### Estimating effects of NPIs on transmission In general, models that included different sets of covariates (i.e. autoregression of Reff, additional temporal marker for county-time, county-level demographic characteristics and NPIs; see details in Table S1 and S2) were fitted with generalized estimating equations (GEEs; *geeglm* of package “geepack”) and autoregressive lag model regressions, separately. We included AR(1) in GEEs and included log Reff in previous week in lag model regression to account for autoregression of Reff. We weighted log-Reff with its inverse coefficient of variation across the above-mentioned 100 stochastic reconstructions. NPI were included as a time-varying covariates with the status of the NPI defined for each week of the analysis as either 1 if in use in that county in a particular week or 0 if not. Table S1 describes the components of the **main** model and alternatives, **base**, the model with the minimum number of covariates that we considered, **time**, which included the covariates included in **base**, but adding time since the first case in each county as a categorical variable and **time and interventions**, which adds NPI to the **time** model. In our main results, we interpreted the estimated effects of NPIs from the main model that was fitted with GEEs (Fig. 3A). Autoregression of Reff, county-level characteristics and NPIs were included in the main model (Table S1). To examine the reduction of different combinations of NPIs on Reff (Fig. 4), we calculated Reff by combining effects for individual NPIs that were estimated from the main model. To further examine whether the effects of NPI suites were robust to interactions between NPIs, we 1) fitted another model with GEEs by including autoregression of Reff, county-level characteristics and NPIs suites (as categorical variables); and 2) fitted a model with XGBoost of individual NPIs (details described below). We then compared the predictions on the Reff for NPIs suites from the GEE suite model and XGBoost model to those calculated from our main model (Fig. 4B). ### Sensitivity analyses of effects of NPIs on transmission #### LASSO To examine the relative importance of NPIs in explaining variance of Reff with increasing parsimony of the model, we performed LASSO (*glmnet* in “glmnet” package) with the main model that was fitted with our lag model regression. The estimated effects of NPIs from lag model regression are highly correlated with those estimated from GEEs (Fig. 3A). Effect sizes of NPIs were presented as the complexity penalizing hyperparameter (A) was increased from 0.0 to 0.1 at intervals of 0.005 (Fig. 3B). Coefficients estimated when A is 0 (i.e. no penalty) are equivalent to the estimates from lag model regressions shown in Fig. 3A. #### Collinearity of covariates We assessed the potential impact of colinearity of NPIs on our estimated associations of NPI with Reff through two approaches: dropping one NPI at a time (“dropped model” hereafter) prior to fitting the model, and by fitting a single-intervention model. For the dropped model, we dropped individual covariates of NPIs and county characteristics one at a time and reran the main model that was fitted to lag model regression. We calculated the ratio of relative changes in Reff that was estimated from the main analysis to that estimated from the dropped models (Fig. 3C). For the single-intervention model, we added one NPI into the base model at a time and compared the estimated effect of that NPI with the original effect size that was estimated from the main model. Single-intervention models were also fitted to GEEs and lag model regression (Fig. S5). #### Impact of temporal declining trend of R*eff**on estimating effect of NPIs* In competing NPIs against temporal markers, we further included a fixed-effect associated with county-time (i.e. time since a county had its first case) into the main model as a weekly discretized coefficient that was shared across all counties (Table S1 and S2). Weeks since the first case in a county were computed for each data point and were included in the models as categorical variables, i.e. from three weeks before to thirteen weeks after the county saw its first case. In order to examine whether the effect of NPIs estimated from our main models was not just capturing the declining trends of Reff over time, we permuted NPIs spatially (i.e. permuted the NPIs suites between counties) and temporally (i.e. permuted the time series of NPIs within counties) and refit the main model. We performed the permutations tests 100 times. #### Impact of uncertainty from early transmission phase Reff estimates at the beginning of an outbreak are often challenged by large uncertainties. In order to examine the effects of these uncertainties on our estimated effects of NPIs, we fitted the main model to a subset of our data set, which only includes Reff since two weeks after the county saw its first case. ### Effects of NPIs relaxations on transmission In order to look at the effects of relaxing NPIs on the transmission, we fitted the main model by further splitting the effect of NPIs into “during the implementation” (intervention on) and “after the implementation was lifted” (intervention off). Relaxations were available for non-essential business, stay at home and medical service suspension. The rest of the covariates were the same as in the main model. ### Using deaths as alternative data to estimate Reff Reff that were estimated from the deconvoluted COVID-19 deaths were derived to assess the robustness of our statistical inference to using these data compared to confirmed cases. We therefore used the above-mentioned methods to estimate weekly Reff from confirmed deaths and refit the main models with GEEs and lag model regression (Fig. S9). Reff were estimated for 1,840 out of 3,142 counties (58.6%). ### Out-of-sample prediction Models were fitted with each of the fifty states and District of Columbia held out as test sets. Prediction performances were measured using root mean squared error (RMSE), mean absolute scaled error (MASE), and coefficient of determination (R2). Fitting procedures for the lag models were as described above. #### Comparative model We employed XGBoost (*23*), a decision tree boosting package, to explore whether more predictive power can be gained through complex model structures. Optimal values of three main hyperparameters (fraction of covariates included in each boosting iteration, fraction of training data included per those iterations, and maximum tree depth) were determined through grid search; ranges (and grid intervals) were 0.3 to 0.9 (0.1), 0.3 to 0.9 (0.1), and 3 to 9 (1), respectively. Performances were evaluated under 10-fold cross-validation. Learning rate was conservatively set to 0.2 and the maximum number of iterations was capped at 200 with early stopping if RMSE does not improve after two iterations to avoid overfitting. The respective optimal values were 0.9, 0.3, and 6. We further optimized the maximum iterations cap when test sets were held out by spatial units; range of 50 to 200 at intervals of 25. Results reported were from the optimal iteration of 75 (which minimized RMSE). ### Effects of NPIs mediated by workplace attendance Google Community Mobility Reports ([https://www.google.com/covid19/mobility](https://www.google.com/covid19/mobility)) were downloaded on September 16, 2020. Mean of daily percentage change of work commutes relative to the baseline of each county were computed for each weekly interval where we have Reff estimates. We chose to focus on workplace presentation among candidate datasets available through Google as this was the only dataset that had less than 50% of county-weeks missing. Mediation analyses were conducted for each NPI separately. For each NPI, we fitted a full model to log transformed weekly Reff and adjusted for county-level characteristics, the workplace attendance and the examined NPI. We then fitted a mediation model which regresses the workplace attendance on the examined NPI. The mediation analyses were conducted using R package “*mediation*” (*24*). ## Supporting information Supplemental Material [[supplements/221036_file02.docx]](pending:yes) ## Data Availability Data and code for analysis are made available at https://github.com/datcummings/SARS-CoV-2.NPI.analysis. [https://github.com/datcummings/SARS-CoV-2.NPI.analysis](https://github.com/datcummings/SARS-CoV-2.NPI.analysis) ## Disclaimers Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc. for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the USDOE or the United States Government. Contributions from Sandia National Laboratories were partially funded by the Laboratory Directed Research and Development (LDRD) program, and the DOE Office of Science through the National Virtual Biotechnology Laboratory, a consortium of DOE national laboratories focused on response to COVID-19, with funding provided by the Coronavirus CARES Act. Material has been reviewed by the Walter Reed Army Institute of Research. There is no objection to its presentation and/or publication. The opinions or assertions contained herein are the private views of the author (AH), and are not to be construed as official, or as reflecting true views of the Department of the Army or the Department of Defense. ## Footnotes * **UFCOVID Interventions Team** Stephanie Khoury, Meenal Indra, Hung-Leong Yau, Ben Cummings, Peter Giannas, Martha-Grace McLean, Ken Hubbard, Camazia Saunders, Caroline Weldon, Caroline Phillips, David Rosenbaum, Dianelys Tabla * * note these numbers round to identical estimates but are different when considering more significant digits * Received October 29, 2020. * Revision received October 29, 2020. * Accepted November 3, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.nytimes, nytimes/covid-19-data, (available at [https://github.com/nytimes/covid-19-data](https://github.com/nytimes/covid-19-data)). 2. 2. K. A. Auger, S. S. Shah, T. Richardson, D. Hartley, M. Hall, A. Warniment, K. Timmons, D. Bosse, S. A. Ferris, P. W. Brady, A. C. Schondelmeyer, J. E. Thomson, Association Between Statewide School Closure and COVID-19 Incidence and Mortality in the US. JAMA (2020),, doi:10.1001/jama.2020.14348. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.14348&link_type=DOI) 3. 3. I. F. Miller, A. D. Becker, B. T. Grenfell, C. J. E. Metcalf, Disease and healthcare burden of COVID-19 in the United States. Nat. Med. (2020), doi:10.1038/s41591-020-0952-y. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0952-y&link_type=DOI) 4. 4. H. J. T. Unwin, S. Mishra, V. C. Bradley, A. Gandy, T. A. Mellan, H. Coupland, J. Ish- Horowicz, M. A. C. Vollmer, C. Whittaker, S. L. Filippi, X. Xi, M. Monod, O. Ratmann, M. Hutchinson, F. Valka, H. Zhu, I. Hawryluk, P. Milton, K. E. C. Ainslie, M. Baguelin, A. Boonyasiri, N. F. Brazeau, L. Cattarino, Z. M. Cucunubá, G. Cuomo-Dannenburg, I. Dorigatti, O. D. Eales, J. W. Eaton, S. L. van Elsland, R. G. FitzJohn, K. A. M. Gaythorpe, W. Green, W. Hinsley, B. Jeffrey, E. Knock, D. J. Laydon, J. Lees, G. Nedjati-Gilani, P. Nouvellet, L. C. Okell, K. V. Parag, I. Siveroni, H. A. Thompson, P. Walker, C. E. Walters, O. J. Watson, L. K. Whittles, A. Ghani, N. M. Ferguson, S. Riley, C. A. Donnelly, S. Bhatt, S. Flaxman, State-level tracking of COVID-19 in the United States,, doi:10.1101/2020.07.13.20152355. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNy4xMy4yMDE1MjM1NXYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTEvMDMvMjAyMC4xMC4yOS4yMDIyMTAzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. 5.State-level variation of initial COVID-19 dynamics in the United States: The role of local government interventions,, doi:10.37473/dac/10.1101/2020.04.14.20065318. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.37473/dac/10.1101/2020.04.14.20065318&link_type=DOI) 6. 6.Google COVID-19 Community Mobility Reports, (available at Google LLC “Google COVID- 19 Community Mobility Reports”. [https://www.google.com/covid19/mobility/](https://www.google.com/covid19/mobility/)). 7. 7. S. Flaxman, S. Mishra, A. Gandy, H. J. T. Unwin, T. A. Mellan, H. Coupland, C. Whittaker, H. Zhu, T. Berah, J. W. Eaton, M. Monod, Imperial College COVID-19 Response Team, A. C. Ghani, C. A. Donnelly, S. M. Riley, M. A. C. Vollmer, N. M. Ferguson, L. C. Okell, S. Bhatt, Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature (2020), doi:10.1038/s41586-020-2405-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2405-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32512579&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F03%2F2020.10.29.20221036.atom) 8. 8. S. Hsiang, D. Allen, S. Annan-Phan, K. Bell, I. Bolliger, T. Chong, H. Druckenmiller, L. Y. Huang, A. Hultgren, E. Krasovich, P. Lau, J. Lee, E. Rolf, J. Tseng, T. Wu, The effect of large-scale anti-contagion policies on the COVID-19 pandemic. Nature (2020), doi:10.1038/s41586-020-2404-8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2404-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32512578&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F03%2F2020.10.29.20221036.atom) 9. 9. J. Zhang, M. Litvinova, Y. Liang, Y. Wang, W. Wang, S. Zhao, Q. Wu, S. Merler, C. Viboud, W. Vespignani, M. Ajelli, H. Yu, Changes in contact patterns shape the dynamics of the COVID-19 outbreak in China. Science. 368, 1481–1486 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjgvNjQ5OC8xNDgxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTEvMDMvMjAyMC4xMC4yOS4yMDIyMTAzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 10. 10. H. M. Korevaar, A. D. Becker, I. F. Miller, B. T. Grenfell, C. J. E. Metcalf, M. J. Mina, Quantifying the impact of US state non-pharmaceutical interventions on COVID-19 transmission. Epidemiology (2020). 11. 11. V. E. Pitzer, M. Chitwood, J. Havumaki, N. A. Menzies, S. Perniciaro, J. L. Warren, D. M. Weinberger, T. Cohen, The impact of changes in diagnostic testing practices on estimates of COVID-19 transmission in the United States. medRxiv (2020), doi:10.1101/2020.04.20.20073338. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNC4yMC4yMDA3MzMzOHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTEvMDMvMjAyMC4xMC4yOS4yMDIyMTAzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 12. 12. S. T. Ali, L. Wang, E. H. Y. Lau, X.-K. Xu, Z. Du, Y. Wu, G. M. Leung, B. J. Cowling, Serial interval of SARS-CoV-2 was shortened over time by nonpharmaceutical interventions. Science (2020), doi:10.1126/science.abc9004. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwNy8xMTA2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTEvMDMvMjAyMC4xMC4yOS4yMDIyMTAzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 13. 13. K. Walker, tidycensus: Load US Census Boundary and Attribute Data as “tidyverse” and “sf”-Ready Data Frames (2020), (available at [https://CRAN.R-project.org/package=tidycensus](https://CRAN.R-project.org/package=tidycensus)). 14. 14. J. Raifman, K. Nocka, D. Jones, J. Bor, S. Lipson, J. Jay, P. Chan, COVID-19 US state policy database, (available at [www.tinyurl.com/statepolicies](http://www.tinyurl.com/statepolicies)). 15. 15. N. G. Becker, L. F. Watson, J. B. Carlin, A method of non-parametric back-projection and its application to AIDS data. Stat. Med. 10, 1527–1542 (1991). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.4780101005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1947509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F03%2F2020.10.29.20221036.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GH68700004&link_type=ISI) 16. 16. S. Abbott, J. Hellewell, R. N. Thompson, K. Sherratt, H. P. Gibbs, N. I. Bosse, J. D. Munday, S. Meakin, E. L. Doughty, J. Y. Chun, Y.-W. D. Chan, F. Finger, P. Campbell, A. Endo, C. A. B. Pearson, A. Gimma, T. Russell, CMMID COVID modelling group, S. Flasche, A. J. Kucharski, R. M. Eggo, S. Funk, Estimating the time-varying reproduction number of SARS-CoV-2 using national and subnational case counts. Wellcome Open Res. 5, 112 (2020). 17. 17. A. Wächter, L. T. Biegler, On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming. Mathematical Programming. 106 (2006), pp. 25–57. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10107-004-0559-y&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234024700002&link_type=ISI) 18. 18. W. E. Hart, C. D. Laird, J.-P. Watson, D. L. Woodruff, G. A. Hackebeil, B. L. Nicholson, J. D. Siirola, Pyomo — Optimization Modeling in Python. Springer Optimization and Its Applications (2017),, doi:10.1007/978-3-319-58821-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-3-319-58821-6&link_type=DOI) 19. 19. S. Sanche, Y. T. Lin, C. Xu, E. Romero-Severson, N. Hengartner, R. Ke, The Novel Coronavirus, 2019-nCoV, is Highly Contagious and More Infectious Than Initially Estimated,, doi:10.1101/2020.02.07.20021154. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMi4wNy4yMDAyMTE1NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTEvMDMvMjAyMC4xMC4yOS4yMDIyMTAzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 20. 20. Q. Bi, Y. Wu, S. Mei, C. Ye, X. Zou, Z. Zhang, X. Liu, L. Wei, S. A. Truelove, T. Zhang, W. Gao, C. Cheng, X. Tang, X. Wu, Y. Wu, B. Sun, S. Huang, Y. Sun, J. Zhang, T. Ma, J. Lessler, T. Feng, Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study. The Lancet Infectious Diseases. 20 (2020), pp. 911–919. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1473-3099(20)30287-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F03%2F2020.10.29.20221036.atom) 21. 21. A. Anzai, T. Kobayashi, N. M. Linton, R. Kinoshita, K. Hayashi, A. Suzuki, Y. Yang, S.-M. Jung, T. Miyama, A. R. Akhmetzhanov, H. Nishiura, Assessing the Impact of Reduced Travel on Exportation Dynamics of Novel Coronavirus Infection (COVID-19). J. Clin. Med. Res. 9 (2020), doi:10.3390/jcm9020601. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm9020601&link_type=DOI) 22. 22. R. Verity, L. C. Okell, I. Dorigatti, P. Winskill, C. Whittaker, N. Imai, G. Cuomo-Dannenburg, H. Thompson, P. G. T. Walker, H. Fu, A. Dighe, J. T. Griffin, M. Baguelin, S. Bhatia, A. Boonyasiri, A. Cori, Z. Cucunubá, R. FitzJohn, K. Gaythorpe, W. Green, A. Hamlet, W. Hinsley, D. Laydon, G. Nedjati-Gilani, S. Riley, S. van Elsland, E. Volz, H. Wang, Y. Wang, W. Xi, C. A. Donnelly, A. C. Ghani, N. M. Ferguson, Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect. Dis. 20, 669–677 (2020). [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%2F11%2F03%2F2020.10.29.20221036.atom) 23. 23. T. Chen, C. Guestrin, in Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (Association for Computing Machinery, New York, NY, USA, 2016), KDD ‘16, pp. 785–794. 24. 24. D. Tingley, T. Yamamoto, K. Hirose, L. Keele, K. Imai, mediation: R Package for Causal Mediation Analysis. Journal of Statistical Software, Articles. 59, 1–38 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v059.i10&link_type=DOI)