key: cord-0908632-j1gi1cxu
authors: Xu, Jie; Hussain, Sabiha; Lu, Guanzhu; Zheng, Kai; Wei, Shi; Bao, Wei; Zhang, Lanjing
title: Associations of Stay-at-Home Order and Face-Masking Recommendation with Trends in Daily New Cases and Deaths of Laboratory-Confirmed COVID-19 in the United States
date: 2020-07-08
journal: Explor Res Hypothesis Med
DOI: 10.14218/erhm.2020.00045
sha: 3479fa33e4899f9780947ab382977b3886ff6d5e
doc_id: 908632
cord_uid: j1gi1cxu

BACKGROUND AND OBJECTIVES: Public health interventions have reduced coronavirus disease 2019 (COVID-19) transmission in several countries, but their impacts on COVID-19 epidemics in the USA are unclear. We examined associations of stay-at-home order (SAHO) and face-masking recommendation with COVID-19 epidemics in the USA. METHODS: In this quasi-experimental interrupted time-series study, we modeled temporal trends in daily new cases and deaths of laboratory-confirmed COVID-19 cases, and COVID-19 time-varying reproduction numbers in the USA between March 1 and April 20, 2020. In addition, we conducted simulation analyses. RESULTS: The number of residents under SAHO increased since March 19 and plateaued at 290,829,980 (88.6% of the U.S. population) on April 7. Trends in COVID-19 time-varying reproduction numbers peaked on March 23, further reduced on April 3, and fell below/around 1.0 on April 13. Early-implementation and early-lift of SAHO would reduce and increase COVID-19 epidemics, respectively. Multivariable piecewise log-linear regression revealed the states’ neighboring relationship with New York was linked to COVID-19 daily new cases and deaths. There were two turning points in daily new-case trend, being March 28 (slope-changes = −0.09) and April 3 (slope-changes = −0.09), which appeared to be associated with implementation of SAHO on March 28 (affecting 48.5% of the US population in 22 states and District of Columbia), and face-masking recommendation on April 3, respectively. There were also two turning points in daily new-death trend, being April 9 (slope-changes = −0.06) and April 19 (slope-changes = −0.90). CONCLUSIONS: We identified two turning points of COVID-19 daily new cases or deaths in the USA, which seem to be linked to implementation of SAHO and the Center for Disease Control’s face-masking recommendation.

As of July 3, 2020, the coronavirus disease 2019 (COVID-19) has been diagnosed in more than 2.7 million people and led to 128,024 deaths in the USA. 1-3 Several blind spots underlying these staggering numbers have been revealed and discussed. 4 In response, 40 states and the District of Columbia (DC) implemented a stay-athome order (SAHO). 5, 6 The Centers for Disease Control and Prevention (CDC) also recommended face-masking on April 3, 2020. 7 It had been demonstrated previously that public health interventions were associated with reduction of the systemic acute respiratory syndrome-coronavirus (SARS-CoV-2) transmission in several countries, [8] [9] [10] [11] but the associations of SAHO and face-masking recommendation with COVID-19 epidemics in the USA are unclear.

Therefore, we examined these associations using observed population data and then performed simulations for outcomes under the scenarios if early-implementation and lift of SAHO occurred. Without the data on these associations, the roles of SAHO and face-masking in controlling COVID-19 pandemics will remain unclear, and the public health officials and state governments will be poorly informed about the best timing, best modalities (or their combination), and best length to implement these interventions.

Quasi-experimental interrupted time series study is the comparative analysis of outcome variable(s) in a cohort with longitudinal (time series) "panel" data sets before and after an intervention. This approach is widely used in situations when randomization is not possible or extremely difficult. 12, 13 Due to the technical difficulty in randomization and proper control of policy implementation for COVID-19, this study design appeared to be the best fit to address the impact of COVID-19 public health interventions. Therefore, we conducted a quasi-experimental interrupted time series study to compare the changes in COVID-19 epidemics before and after SAHO and face-masking recommendation. We extracted the national and statelevel daily new cases and deaths of COVID-19 from the COVID- 19 Tracking Project, which has tracked COVID-19 data since February 28, 2020. 2 Only those cases and deaths that occurred from March 1 to April 20, 2020 in the 50 states and DC were analyzed.

New cases and deaths were defined as laboratory-confirmed positive cases or deaths which were reported by a state's public health authority for better data consistency and quality. 2 Each of these state authorities reported its data in different format, while most, if not all, of them followed the reporting guidelines of the CDC. It is noteworthy that on April 14, 2020, the CDC updated its definition of positive cases and included probable-positive cases. 3 The impact of this change to the data released on April 20, 2020 would be minimal due to the use of confirmatory laboratory-test results and the short time-interval (many states have not adopted the criteria yet). Therefore, the case/death numbers reported here might be smaller than those reported by others after April 14.

For quality control of the released and curated data, the COV-ID-19 Tracking Project employed a 4-tier score system, which included four simple components, namely, reporting positive test results reliably, reporting negatives sometimes, reporting negatives reliably and reporting all commercial tests. Based on the sum of these scores, each state corresponded to a letter grade A, B, C or D, with A for the best quality. All states scored A or B.

This study was exempt (Category 4) from an Institutional Review Board review owing to the use of data from a publicly available, de-identified database.

Several population-based factors were included in the multivariable piecewise log-linear regression analyses. Specifically, the timing of SAHO and populations of the states were obtained to calculate the number of subjects and the proportion of the USA population under SAHO on a given date. 5 The state populations were extracted from the USA Census (up to July 2019). 14 To adjust for confounders, the proportions of daily positive results in all daily tests, and state-level daily new cases and deaths were obtained from the COVID-19 Tracking Project. 2

The time-varying reproduction number (Rt) was defined as the mean number of secondary cases generated by a typical primary case at the time (t) in a population, and estimated using previously-reported mean serial-intervals 15-17 and the R package (version 3.6.3). 8, 18 The mean serial interval was defined as the average time from the illness onset to successive cases in the chains of transmission. 15-17 Three-day moving averages of the Rt were reported for their better sensitivity than 5-day moving averages, due to the greater amount of data points (data not shown). We estimated the segmental coefficients using piecewise log-linear models and two presumptive turning-points, including the day as the unit for x axis (i.e. the per-day change in the logarithmic values of COVID-19 daily new cases or deaths). Simulation analysis was performed using the prediction function. All statistical analyses were performed using Stata (version 15) and the Joinpoint program (NCI, version 4.7.0.0) with the Poisson variance option. All p-values were twosided, with a cut-off of 0.05 for significance.

Considering the recently reported state heterogeneity of COV-ID-19 epidemics in the European Union, 19 we examined the association between the states' neighboring relationship with New York (epicenter at the time) and COVID-19 epidemics (daily COVID-19 new cases or new deaths). The relationship was coded as 0 for New York, 1 for the six neighbor states (Connecticut, Massachusetts, New Jersey, Pennsylvania, Rhode Island, and Vermont), and 2 for the other states and DC. We compared these epidemics using Wilcoxon rank sum test and Student's t-test. We then conducted sensitivity studies using multivariable piecewise log-linear regression with the states' relationship with New York as a covariate. The COVID-19 trend lines were smoothed in the figure using a locally weighted scatter smoothing (known as LOWESS) algorithm. 20,21

On March 19, 2020, the state of California started a SAHO which affected 39,512,223 residents (12.0% of the U.S. population) (Fig.  1) . The number and proportion of residents under the SAHO subsequently implemented by other 39 states and DC continued to increase until April 7, and plateaued at 290,829,980 and 88.6%, respectively. Approximately 85.2% of the U.S. population were under SAHO on April 3, 2020, while only 3.4% of the U.S. population were added from April 3 to April 25. The number of the states that implemented a state-wide SAHO peaked at 41 on April 25, 2020, and some states started to lift the state-wide SAHO on April 26, 2020 ( Table 1) .

The log-linear models used by Joinpoint program and Stata identified similar turning-points. In the multivariable piecewise Our simulation analysis shows early-implementation of SAHO would be associated with a significant reduction in daily new cases and deaths, while lift of SAHO would be associated with a significant increase in daily new cases and deaths (Fig. 2) . The estimates of Rt based on the three reported mean serial-intervals of COVID-19 all started to decline on March 19, when SAHO was first implemented in the USA, and declined faster after March 23 (Fig. 3) . After a short plateau, Rt continued to decline after April 3 and fell below/around 1.0 on April 13.

We found that New York state and its neighbor states had higher daily new cases and deaths than other states ( 

The population under SAHO in the USA grew from March 19, 2020 to April 6, 2020, and reached 290,829,980 (88.6%) by April 7. The multivariable piecewise log-linear regression models identified two turning points of COVID-19 daily new cases and the time-varying reproduction number, Rt, in the USA, as well as the link between states' relationship with New York and COVID-19 epidemics. The two turning points were associated with implementation of SAHO on March 28 (affecting 48.5% of USA population in 22 states and DC) and the CDC's face-masking recommendation on April 3. Similar turning points of COVID-19 daily new deaths were April 3 and April 15, which represented 10 and 12 days of delay, respectively. Simulation on early-implementation and lift of SAHO also reveals considerable impacts on COVID-19 daily new cases and deaths.

Several recent studies showed that public health interventions, including school closure, cordons sanitaire, traffic restriction, social distancing and others, were linked to reduction of Rt and daily incidence of COVID-19. 8-11 However, these data were mostly from other countries where the socioeconomic systems and the modalities and extents of public health interventions were different from those in the USA. 5,6 Among several modalities of social distancing, including school closure, gathering restrictions and restaurant restriction, 6 the implementing time of SAHO appeared to be the only factor that matched to the identified turning point (March 28) in this study. Moreover, the Berkeley Interpersonal Contact Study also showed reduction in interpersonal contact in the USA between Mach 22 and April 8. 22 Therefore, our and others' data together suggest that the implementation of SAHO, which is the strictest social distancing modality, may be required to reach the turning point/timing in trends of new cases. In addition, the studies in China also support our findings that travel restriction alone would reduce SARS-CoV2 transmission, 8-10 but additional reduction in transmission in the community (e.g., facemasking in the USA) may be required to attenuate or reverse the epidemic trajectory.

The effects of masking on the epidemic of COVID-19 in the USA was simulated using the filtering efficacy of masks on influenza. 23, 24 The influence of public intervention on COVID-19 Rt is one of the most widely used metrics for assessing transmission rate of infectious diseases, 8,25 and linked to the incidence decay with exponential adjustment model (known as IDEA) and Farr's law. 32 However, Rt is difficult to estimate. One of the challenges is the variances in the SARS-CoV2 serial intervals. [15] [16] [17] The reported serial intervals of COVID-19 ranged from 3.96 to 7.5 days, and were first used to rigorously examine the changes of Rt associated with SAHO and face-masking recommendation in the USA. Two similar turning points were identified in the Rt trends estimated using three different serial intervals, and further support the findings discovered using our piecewise log-linear model.

Geographic attributes have been shown to be linked to COV-ID-19 epidemics in several countries. 33-36 The geographic differences of COVID-19 incidence and deaths in states have been reported in a previous study, 37 which, unfortunately, did not conduct quantitative trend analyses. The states' neighboring relationship with the epicenter (New York) in the USA was linked to the COV-ID-19 epidemics in our multivariable models, while the number of states under SAHO, proportion of positive COVID-19 testing and the number and proportion of the populations under SAHO were not linked to them. Additional works are needed to better understand the geographic differences in state-level COVID-19 trends.

Many COVID-19 studies were published as preprints. Some of them were of inferior methodological merit and low reporting quality, as discussed before. 30, 38, 39 This observation is concerning, and calls for more collaborative efforts in the efficient review and timely dissemination of studies on COVID-19.

The major strength of this study is that we used state-level and national data of laboratory-confirmed cases in the analyses. 2 Moreover, a very rigorous study-design was used. We applied multivariable models to exclude confounders and a quasi-experimental interrupted time series study-design to infer the impacts of public health interventions. The quasi-experimental interrupted time series study-design is likely the most rigorous method available for when randomization and controlling are not possible, 12,13,19 as the case for when COVID-19 SAHO and face-masking were implemented in the USA. Further, two piecewise log-linear regression methods were used to rigorously examine trends changes, according to the guidelines on trend analyses of population data and other methodological considerations. 26-29 Finally, the simulation analysis provides comprehensive estimates of trends changes linked to early-implementation of SAHO at various time points and early-lift of SAHO with various extents. Indeed, the recent upward trends in the USA COVID-19 epidemics are consistent with our simulated scenarios of early SAHO lift (even in some states or counties). 2,3,31 Thus, the proper timing and duration of SAHO implementation warrant further research.

Several limitations of this study are noteworthy. First, the COVID-19 positive rates varied among states and by time, suggesting under-testing of the potential patients. The exact COV-ID-19 case numbers are thus not available, although efforts were made to estimate them using the Johns Hopkins' data repository of COVID-19 cases. 40 Given the data inconsistence we noticed, 1 such an estimation was not optimal in our view. We were more confident in the reliability of the laboratory-confirmed case numbers. Inclusion of positive-test rate in our models may alleviate the variances in test rate across time. Indeed, the positive-test rate was not a factor linked to COVID-19 epidemics in our multivariable models. Second, there was a lag in COVID-19 reporting, 41 which may lead to inaccurate estimation of the case numbers. However, the increases in the proportion of the tested population appeared stable in the USA, 2 and was not a factor linked to COVID-19 epidemics in our multivariable models. It suggests that the lag in reporting may not change significantly as the time changes, and will have minimal impact on the trend analyses. Third, the CDC changed the diagnostic criteria of COVID-19 on April 14, 2020. The timing and extent of these changes in each state were largely unknown, while we consider the data released before April 20 acceptable owing to the transition time of different criteria. Therefore, the data released after April 20 would have an unmeasurable bias and would not be suitable for comparing with those before April 20. Thus, we were not able to reliably compare the observed and the simulated COVID-19 trends after April 20. Finally, we did not report the daily incidence. There were no significant changes in the USA population during the study-period. The daily new cases of COVID-19 thus should be proportional to its daily incidence in the USA but are easier to interpret than daily incidence and were used here.

Future research should be focused on how state-level SAHO changes the COVID-19 epidemics. This is particularly interesting since state heterogeneity has been reported in 27 European Union states, 19 and in the USA, states' neighboring relationship with New York was linked to COVID-19 epidemics (as this study showed). As more data and more-reliable data become available, we will keep looking into the association of SAHO and face-masking policies with changes in the COVID-19 epidemic.

There were two turning points of COVID-19 daily new cases or deaths in the USA, which appeared to be associated with the implementation of SAHO and the CDC's face-masking recommendation. Simulation on early-implementation and lift of SAHO reveals considerable impact of these interventions on COVID-19 daily new cases and deaths. These findings may inform decision-making of lifting SAHO and face-masking recommendation.

Xu J. et al: Policy and COVID-19 trends Explor Res Hypothesis Med

All raw data used here are publicly available at their respective websites. No additional data are available.

The work was supported in part by a grant from the National Institutes of Health (R03 HD100708 to WB).

Trends and Prediction in Daily New Cases and Deaths of COVID-19 in the United States: An Internet Search-Interest Based Model

The COVID Tracking Project

CDC. Cases of Coronavirus Disease (COVID-19) in the U.S. Available from

Blind Spots in Fighting the Outbreak of Coronavirus Disease 2019

state and local government response to the COVID-19 pandemic

Pandemic Politics: Timing State-Level Social Distancing Responses to COVID-19

Especially in Areas of Significant Community-Based Transmission

Association of Public Health Interventions With the Epidemiology of the COVID-19 Outbreak in Wuhan, China

Changes in contact patterns shape the dynamics of the COVID-19 outbreak in China

The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak

Impact of climate and public health interventions on the COVID-19 pandemic: a prospective cohort study

Quasi-experimental study designs series-paper 5: a checklist for classifying studies evaluating the effects on health interventions-a taxonomy without labels

Regression based quasi-experimental approach when randomisation is not an option: interrupted time series analysis

Census. State Population Totals and Components of Change

Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia

Serial interval of novel coronavirus (COVID-19) infections

Serial Interval of COVID-19 among Publicly Reported Confirmed Cases

A new framework and software to estimate time-varying reproduction numbers during epidemics

State heterogeneity of human mobility and COVID-19 epidemics in the European Union

Speaking Stata: Smoothing in various directions

Robust locally weighted regression and smoothing scatterplots

Quantifying interpersonal contact in the United States during the spread of COVID-19: first results from the Berkeley Interpersonal Contact Study

To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic

Testing the efficacy of homemade masks: would they protect in an influenza pandemic

Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel Coronavirus

National Center for Health Statistics Guidelines for Analysis of Trends. National Center for Health Statistics. Vital and Health Statistics (series

Statistics and pitfalls of trend analysis in cancer research: a review focused on statistical packages

Trends in Alcoholic Fatty Liver Disease

Trend analysis of diabetic mortality

Trends and prediction in daily incidence of novel coronavirus infection in China, Hubei Province and Wuhan City: an application of Farr's law

A Weekly Surveillance Summary of U.S. COVID-19 Activity

Relatedness of the incidence decay with exponential adjustment (IDEA) model

Epidemiology of 2019 novel coronavirus in Jiangsu Province, China after wartime control measures: A population-level retrospective study

Epidemiology of the 2020 Pandemic of COVID-19 in the State of Texas: The First Month of Community Spread

Geographic risk assessment of COVID-19 transmission using recent data: An observational study

Effect of weather on COVID-19 spread in the US: A prediction model for India in 2020

Geographic Differences in COVID-19 Cases, Deaths, and Incidence -United States

Medical Preprints-A Debate Worth Having

Will research preprints improve healthcare for patients?

Estimating the number of SARS-CoV-2 infections in the United States

lags in coronavirus testing after slow response to outbreak

This work has been deposited at MedRxiv https://doi.org/10.1101 /2020.05.01.20088237.

No disclosures were reported.

Supervision (LZ), concept and design (JX, LZ), drafting of the manuscript (JX), statistical analysis (LZ), acquisition, analysis, or interpretation of data (JX, SH, GL, KZ, SW, WB, LZ), critical revision of the manuscript for important intellectual content (all authors). Dr. Zhang (LZ) had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.