key: cord-0952502-zsp8l3sp authors: Jinjarak, Yothin; Ahmed, Rashad; Nair-Desai, Sameer; Xin, Weining; Aizenman, Joshua title: Accounting for Global COVID-19 Diffusion Patterns, January–April 2020 date: 2020-09-04 journal: Econ Disaster Clim Chang DOI: 10.1007/s41885-020-00071-2 sha: c2c1a20eb04f22cc59a775bdaf2e411be5a5f626 doc_id: 952502 cord_uid: zsp8l3sp Key factors in modeling a pandemic and guiding policy-making include mortality rates associated with infections; the ability of government policies, medical systems, and society to adapt to the changing dynamics of a pandemic; and institutional and demographic characteristics affecting citizens’ perceptions and behavioral responses to stringent policies. This paper traces the cross-country associations between COVID-19 mortality, policy interventions aimed at limiting social contact, and their interactions with institutional and demographic characteristics. We document that, with a lag, more stringent pandemic policies were associated with lower mortality growth rates. The association between stricter pandemic policies and lower future mortality growth is more pronounced in countries with a greater proportion of the elderly population and urban population, greater democratic freedoms, and larger international travel flows. Countries with greater policy stringency in place prior to the first death realized lower peak mortality rates and exhibited lower durations to the first mortality peak. In contrast, countries with higher initial mobility saw higher peak mortality rates in the first phase of the pandemic, and countries with a larger elderly population, a greater share of employees in vulnerable occupations, and a higher level of democracy took longer to reach their peak mortalities. Our results suggest that policy interventions are effective at slowing the geometric pattern of mortality growth, reducing the peak mortality, and shortening the duration to the first peak. We also shed light on the importance of institutional and demographic characteristics in guiding policy-making for future waves of the pandemic. population mortality rates per capita during the first phase of the pandemic, controlling for policy and structural factors subject to data availability and quality. We plan to revisit these issues with better quality and longer-term data in the coming quarters. A fair share of the countries reached a local peak of the COVID-19 population daily new mortality rate curve during the sample period [see Figs. 1, 2 and 3] . 3 Applying various techniques, we study the factors accounting for the empirical shape of the mortality curve from the onset of the pandemic to the local peak, with a focus on the impact of policy intensity interacting with structural variables. Like similar studies, one should use healthy skepticism in reading the results. First, data quality and availability are a major limitation, as each country has its challenges with data collection, aggregation, and reporting. Second, 'better performance' in the first mitigation phase of a pandemic does not guarantee superior future performance, as the dynamics of a new viral pandemic are yet unknown. By design, flattening the pandemic curve shifts some mortality incidence forward. The susceptibility to secondary waves of infection remains a looming threat. 4 Policies adopted in the second quarter of 2020, and the realized pandemic infections, containment, and treatment will explain the future performance of each country. Furthermore, only time and much more medical research will tell the degree to which infected persons that recovered gained immunity for a long enough period to allow smooth convergence to 'herd immunity. ' Our study relies on daily COVID-19 policy and case data reported by Oxford and John Hopkins University, as well as Apple mobility data and various controls. Our baseline estimation study examines OECD and Emerging Market (EM) sub-samples based on data from 1/23/2020-4/28/2020; or the first 97 days of the pandemic. Below, we summarize the main results. First, we investigate the evolution of weekly mortality growth rates over time and across countries. Applying dynamic panel analysis, local projections (Jordà 2005) suggest that administering more stringent pandemic policies were associated with significantly lower mortality growth rates, with a lag of 2 to 4 weeks, during the first pandemic phase (i.e., during the time from the first death through the end of the sample period). Countries with a Stringency Index (SI) 10 units higher than average had, two weeks later, mortality growth rates that were on average 22 percentage points lower (Oxford's SI is normalized between 0 to 100; where 100 = strictest response). The reductions in mortality growth rate are smaller after three and four weeks, roughly 17 and 13 percentage points respectively. While the reduction in growth rates seem quite large, it is important to put these numbers in perspective. Given the geometric nature of disease spread, oftentimes, week-over-week mortality growth rates can be anywhere from +50% to +100% or greater. Taking slow-moving country fundamentals from the period pre-COVID-19 as exogeneous, we find that countries in the 75th percentile in terms of proportion elderly (65 or older) saw a much stronger reduction in mortality growth rates from the same 10 unit rise in SI, compared to countries with relatively low proportions of elderly (25th percentile). Countries in the 75th (25th) percentile saw mortality growth rates of about −25% (−10%) after two weeks. Countries 3 Using a conservative threshold of countries which reached the peak mortality one week prior to the end of the sample period, roughly 39 countries out of 59 in the sample reached their first peak by the end of the first wave. Additional visualizations on the dates of peak mortalities in the sample and which countries did not reach the peak can be found in the, GitHub, https://github.com/snairdesai/COVID-19/tree/master/output_png. Data and codes available at Github Repository: https://github.com/snairdesai/COVID-19. 4 A "better performance" in the first mitigation phase of a pandemic may induce a too fast opening of the economy towards "business as usual," reducing the public vigilance of keeping social distancing and masking. Hence, we chose to avoid making any arguments related to the overall performance in terms of total death toll. with a greater proportion of the elderly are unconditionally more susceptible to the pandemic, but for this same reason, they are likely to benefit more under stringent policies. In countries with a higher proportion of urban population, SI measures had a stronger impact on mortality growth than in countries with a lower proportion of urban population. This heterogeneity is consistent with the fact that urban populations may face greater exposure due to their high population densities. Greater policy stringency is also more strongly associated with lower mortality growth during the first phase of the pandemic in countries with greater inbound travel, greater GNI per capita, higher health expenditures and a greater level of democracy (measured with the EIU Democracy Index). 5 While international travel flows are quite intuitive risk factors for a pandemic like COVID-19, the role of democratic freedom is an ongoing topic of debate (Ang 2020) . Our results are consistent with the view that greater individual rights may be detrimental in this situation, making it more difficult for the government to administer strict quarantines and have citizens abide by them. Next, we turn to cross-country regression results, where the dependent variables include the logged peak new mortality rate (calculated as the new deaths out of the population at the peak of daily new mortality, by country), and the ratio of 'peak new mortality rate'-to-'pandemic duration to first peak' measured in days (a proxy for flatness/steepness of the mortality rate curve). Definitions of these peak-related dependent variables are illustrated in Fig. 4 using the daily new mortality curve of Czechia. Countries with more aggressive policy interventions in place prior to the first death tend to exhibit a lower new mortality rate at the peak ("Early SI" or "proactive stringency policy"). 6 We find that a one unit increase in Early SI (average SI level prior to the first COVID mortality) is associated with peak new mortality rates on average − 12% lower. Quantitatively, early SI was also associated with flatter mortality curves, but the estimate is statistically insignificant. Unsurprisingly, higher early mortality growth rates (growth rate of new mortality rate in the first week following the first death) are significantly associated with higher peak mortalities and steeper mortality curves, while higher levels of the mortality rate early on (cumulative mortality rate within the first week following the first death) were associated with flatter overall mortality curves. 7 Greater early mobility was also significantly associated with higher mortality peaks. Higher population density is, somewhat surprisingly, associated with flatter mortality curves (column [2] ). However, this association may be driven by the fact that many of the high-density countries are in Asia. Overall, the evidence suggests (but does not necessarily assert) that policy stringency directly reduced the peak mortality rates, and that other forces were also at play (e.g., early mobility, initial pandemic conditions, institutional features). Not only do mortality rates during the first pandemic phase differ across countries, there is also considerable variation in how long new deaths continued to climb in terms of days. We term this as the 'pandemic duration to the first peak' (PD). One should be careful when interpreting the effects of covariates on the PD in terms of altering the shapes of mortality curves, as a longer PD could be accompanied by a higher peak mortality rate and thus a steeper curve, or a lower peak mortality rate and thus a flatter curve. Fitting a Kaplan-Meyer curve for 5 The EIU Democracy Index is based on five categories: electoral process and pluralism; civil liberties; the functioning of government; political participation; and political culture. 6 Higher Early SI measures in some countries, for example, often reflect the use of international travel bans or 'shelter-in-place' policies that were announced in anticipation of a pandemic breakout. 7 Including early mortality growth in the cross-country regression also knocks out the significant effect on elderly population. It's possible that countries with higher early mortality growth rates are also countries with higher proportion elderly. These countries would realize more deaths early on given the at-risk population is larger. the PD over all countries in the sample, in number of days, suggests that countries with stronger Early SI measures (above the average) had significantly lower PDs on the way to the first local peak of the mortality/day curve, compared to countries which did not. The probability of peaking after 40 days is close to 0% for above-average Early SI countries (virtually all of these countries experienced their first local peak by day 40). For below-average Early SI countries, it took approximately 80 days, twice the amount of time, until the first peak was reached. To better understand the cross-country variation in PD, we estimate a Cox proportional hazards model. Higher mortality rates early on are associated with shorter pandemic durations to the peak, while countries realizing higher mortality peaks tend to have, Fig. 1 COVID-19 mortality rate curves, by country. LHS (a): Cumulative logged mortality rate. RHS (b): New Mortality Rate. 7-day rolling averages. Series starts from the 5th COVID-assigned death unsurprisingly, longer pandemic durations to the peak. Although the negative association between the strictness of early policy interventions and the pandemic duration to the peak is not significant under the Cox model, more aggressive policy interventions early on could still be associated with a shorter pandemic duration to the peak through reducing mortality peaks as evidenced by the cross-country analysis. Countries that responded faster prior to the first death had shorter pandemic duration to the peak. Additionally, countries with greater proportions of elderly populations, greater proportions of urban population, greater shares of vulnerable employment tend to exhibit shorter pandemic durations to the peak. 8 Moreover, under certain specifications, the level of democratic freedom appears to be a highly significant determinant in pandemic duration to the peak. That is, countries that have a higher level of democracy saw longer pandemic durations to the peak. Countries further away from the equator also tended to experience longer PD. While at this stage we are reporting suggestive statistical associations, more data and research are needed to get fuller identifications of all these factors. We focus on two aspects of COVID-19 mortality over the first three months of the pandemic spanning from January 23rd, 2020 to April 28th, 2020, which we refer to as the first pandemic phase: first, the dynamic of COVID-19 mortality rates out of the total population, examining the weekly growth rate of the new mortality rate per capita; and second, the empirical shape of the mortality rate curve from the onset to the local peak of the first pandemic phase, examining three dependent variables discussed below. To filter out noise in the daily mortality data, we construct a 7-day rolling average of the daily mortality rate per capita and use these series of 8 This result, which seems counterintuitive, is consistent with our findings in the dynamic panel analysis, where countries exhibiting these same risk factors also had more effective stringency policies in terms of reducing mortality growth rates. Economics of Disasters and Climate Change averages in our estimations and peak identifications. For simplicity, the mortality rate mentioned hereafter is referred to as the 7-day rolling average of the mortality rate. Although examining the growth rate of the new mortality rate per capita provides evidence on how mortality rates evolve over time, it is not trivial to conclude which variables can characterize the cross-country difference in the empirical shape of the mortality curve from the onset to the local peak of the first pandemic phase. We argue that three outcome variables related to the local peak matter: first, the new mortality rate at the peak; second, the pandemic duration to the peak; third, the ratio of the new mortality rate at the peak to the pandemic duration to the peak. We illustrate this idea by comparing country cases. Figure 5 shows the daily new mortality rate curves for several countries: the left one comparing that of Hungary and Norway, the middle one comparing to that of Denmark and Norway, and the last one comparing that of Austria, Estonia, and Greece. The left figure shows that although the peak mortality rates of Hungary and Norway are around the same, their durations to the peak are different: Norway reached the peak faster, and thus had a steeper mortality curve. The middle figure shows that although the durations to the peak of Denmark and Norway are around the same (i.e., both around 25 days), their peak mortality rates are different: Denmark's mortality rate climbed to a higher level before going down, and thus had a steeper mortality curve. Hence, a lower peak mortality rate or a longer duration to the peak implies a flatter mortality curve; however, only on the condition that all other things are held equal. Once both the peak mortality rate and the duration to the peak change in the same direction, it may be ambiguous whether the mortality curve is flattened or not, as demonstrated in the right figure of Fig. 5 . Both Estonia and Greece realized a lower peak mortality rate and a shorter duration to the peak than Austria; however, it is obvious that Estonia has a much steeper mortality curve while Greece has a much flatter mortality curve compared with Austria. This implies that the ratio of the peak mortality rate to the duration to the peak also plays a key role in understanding the empirical mortality curve. Hence, it follows that all three outcome variables related to the peak (including the daily new mortality rate at the peak, the duration to the peak, and the ratio of the peak rate to the duration) together characterize the empirical shape of the mortality curve from the onset to the local peak of the first pandemic phase. By knowing how these three outcome variables are impacted by government pandemic policies or country-specific structural variables can we uncover to what extent these factors account for the pattern of mortality's climb to the local peak of the first pandemic phase. We start with a panel study of mortality growth rate dynamics, using the week-over-week growth rate of the new mortality rate per capita, accounting for containment and closure policy interventions (see Oxford's COVID-19 Government Response Tracker). 9 Specifically, our dependent variable y i, t in country i on date t is defined as where MortalityRate i, t is the new mortality rate in country i on date t. A lower growth rate of the new mortality rate implies a flattening of the mortality curve. Our first benchmark estimation uses the method of local projections (Jordà 2005) , examining future (or current) mortality growth rate as a function of current (or past) mortality growth rate and degree of policy stringency. We aim to understand to what degree policy interventions are associated with future mortality growth, and therefore, the evolution of the pandemic. Local projections do not only simplify our problem but also produce robust estimates under misspecification. 10 Specifically, our model is where y i, t + h is the week-on-week growth rate of the new mortality rate in country i at date t + h for h = {14, 21, 28}. SI i, t is the Stringency Index constructed in the Oxford COVID-19 Government Response Tracker, an aggregate measure of the overall stringency of containment and closure policies, 11 at date t and y i, t − 1 is the one-day lagged mortality growth rate. Fixed effects are denoted as a i and δ t , representing the country and time fixed effects, respectively. We choose to examine the response of new mortality growth with respect to SI at horizons no shorter than two weeks, with reference to studies on the incubation and death periods: The incubation period is 6 days on average (McAloon et al. 2020) , and the death period (number of days from symptom onset to death) ranges from 12 to 15 days for different age groups on average according to CDC estimates. 12 The collection of estimates b γ h ð Þ for h = {14, 21, 28} trace out the dynamic impact of stringency policies on mortality growth at the weekly frequency. 13 Additionally, we study the heterogeneity in the association between policy interventions and mortality growth by estimating the model with interaction terms between the Stringency Index and country-specific social and economic variables where x i is the country-specific variable of interest. We consider: the proportion of the elderly population (people aged 65 and over), the proportion of the urban population, proportion of employment in vulnerable sectors, population density, the logarithm of GNI per capita, health expenditure (% of GDP), population-weighted exposure to ambient PM2.5 pollution, the logarithm of tourist arrivals and departures, level of democracy, and country location measured with latitude and longitude. 14 10 Local projection methods have become widely popular for macroeconomic analysis (see Ramey and Zubairy (2018) for a recent application), specifically for estimating impulse response functions. They are shown to be robust to misspecification and very easily allow for non-linearities (which we introduce when studying heterogeneity). The basic idea is that a series of regressions are estimated, where in each regression the dependent variable is shifted one period forward for h periods. This results in a series of h regressions outputs (referred to as the set of local projections) with h regression coefficient sets, which are used to trace out the dynamic response of the dependent variable from a shock to a particular covariate. 11 Stringency Index (SI) is calculated by taking the ordinal value and adding a weighted constant if the policy is general rather than targeted, if applicable, which are then rescaled by their maximum value to create a score between 0 and 100, with a missing value contributing 0. Hence it captures the average intensity of containment and closure policies. More information can be found at Oxford's Government Response Tracker, https://www. bsg.ox.ac.uk/research/research-projects/coronavirus-government-response-tracker. 12 See https://www.cdc.gov/coronavirus/2019-ncov/hcp/planning-scenarios.html#table-2. 13 Reporting week-on-week estimates instead of daily ensures that the results are estimated on non-overlapping periods. 14 The range of the latitude is from −90 to 90 and the range of the longitude is from −180 to 180, indicating that each country is uniquely identified by an interaction of the latitude and longitude. We follow with a cross-country analysis examining the mortality rate at the peak, a key moment in the first quasi-bell curve, which puts hospitals' capacity to their most severe test. The quasi-bell shapes are normalized by the day of the first significant death. As discussed before, we consider three outcome variables related to the peak: first, the logged peak mortality rate, second, the pandemic duration to the first peak (PD), and third, the ratio of the logged peak mortality rate to the PD. We opt to use linear regression analysis to examine the cross-country difference in the logged peak mortality rate and the ratio of the logged peak mortality rate to the PD and use survival analysis to examine the cross-country difference in the PD. Our cross-country peak mortality data is calculated from the sample from January 23rd, 2020 to April 28th, 2020, during which many OECD countries and emerging market economies finished their ride up to the first peak of the first quasi-bell in terms of contagion per capita and fatality per capita. Our cross-sectional linear regression model is where y i is the logged peak mortality rate or the ratio of the logged peak mortality rate to the PD in country i. A higher peak mortality rate implies a larger inflow of patients, stretching hospitals' capacity. Accounting for the pandemic duration to the first peak, a higher ratio of the logged peak mortality rate to the PD implies a steeper mortality curve, characterized by either a larger or a faster patient inflow that could potentially overwhelm the healthcare system. We include a set of potential endogenous variables Z i . First, one may be interested in whether the cross-country difference in the intensity of the COVID-19 outbreak explains the cross-country difference in the empirical shape of the mortality curve. We include log(Early Mortality i ), the logged cumulative mortality rate in the first week after the first death, and Early Mortality Growth i , the growth rate of daily mortality rate in the first week after the first death, to control for the cross-country difference in the initial level and growth of the mortality rate. Second, one may be interested in whether proactive stringency (i.e., stringency policies in place before the first reported death) policy interventions influence the cross-country difference in the Economics of Disasters and Climate Change empirical shape of the mortality curve. We also include Early SI i , the average of the Stringency Index (SI) from its first non-zero value to the first death, accounting for how strict government interventions were before the first confirmed death, and Days from First SI to First Death i , the number of days from the first non-zero SI to the first death, accounting for how early government interventions are implemented. Additionally, to account for cross-country differences in how aggressively countries respond to the pandemic and increase their policy intensities, we include a variable which we refer to as the SI Delta i (or Stringency Delta), calculated as the difference between a country's maximum level of SI and its initial level of SI (SI i, 0 ), normalized by the number of days (T i ) between them: A higher maximum SI and/or a shorter time to the maximum SI will yield a higher SI Delta. Third, we include Early Mobility i , the weekly average mobility index in terms of walking (reported by Apple) in the week before the first death. We emphasize that while these variables are important to investigate, all of them are endogenous, as they are calculated over part of the first wave period of COVID-19. We also include a set of country-specific control variables X i that we take as exogenous, including the proportion of the elderly population (people aged 65 and over), the proportion of the urban population, proportion of employment in vulnerable sectors, population density, GNI per capita, health expenditure (% of GDP), level of democracy, and country location measured with latitude and longitude. We then proceed with a survival analysis studying the association between the pandemic duration to the first peak (PD) and a set of explanatory variables. We focus on the survival function of the mortality peaking. which is defined as the probability that the PD is later than date t, which is the probability that the mortality peaks after date t. A higher probability that the PD is later than a certain date implies a longer PD, which could have ambiguous implications. On the one hand, it suggests a slower surge in hospitalization that could ease the burden on the healthcare system. On the other hand, a longer timeto-peak may imply a longer-lived, poorly managed pandemic. 15 Our benchmark specification is the Cox proportional hazards model (Cox 1972) , which examines the relationship between the hazard function and a set of explanatory variables. 16 The hazard function is defined as the probability that the peak is on date t conditional on that the peak is reached until date t or later, 15 Without accounting for the peak mortality rate or the cumulative mortality rate, a longer pandemic duration to the first peak could be accompanied by a higher daily new mortality rate at the peak. 16 The assumption of the Cox Model is that the variable of interest has a time-invariant multiplicative effect on the hazard of COVID-19 deaths (see George et al. (2014) for a review of survival methods), which will be verified in section IV. where λ(t| X i ) is the hazard function for country i on date t, conditioning on a set of endogenous variables Z i and exogenous variables X i , and λ 0 (t) is the baseline hazard function. In addition to the same set of endogenous and exogenous variables as in the cross-country regression analysis, we also include the endogenous logged peak new mortality rate, log(Peak Mortality i ), to control for the cross-country difference in the peak level of mortality rates. The specification of the Cox Model implies that the effect on the hazard function of a one-unit increase in one covariate w ∈ {Z, X} with coefficient δ is to multiply the hazard function by e δ , and that the effect on the survival function is to raise it to a power given by the effect on the hazard function where S 1 (t) is the survival function on date t for a group with a one-unit higher value of the covariate w, all other variables held constant. We wish to briefly call out the limitations of our research design. First, our estimates cannot (and should not) be interpreted as causal. What we are reporting, across all models, are associations. Some of our variables are endogenous, which may bias our estimates. Moreover, given our choice to investigate country factors one-by-one, our regression estimates may also be biased from omitted variables. To overcome these challenges, we are in the process of collecting additional data at varying levels of detail to help deal with these issues, with the aim of achieving cleaner identification going forward. Nonetheless, under such data constraints, we believe our approach strikes a balance between parsimony, robustness, and informativeness. Our data relate prudential and reactionary government interventions to COVID-19 mortality rates per capita, controlling for country-specific characteristics. We construct the mortality rates using the John Hopkins Center for Systems Science and Engineering (CSSE) COVID-19 data repository, which details confirmed cases and deaths across our sample period. The seven-day rolling average of the new mortality rate per capita is calculated as the ratio of the cumulative total of country deaths by population, while the seven-day rolling average of the cumulative mortality rate is the ratio of new daily deaths by population. 17 This data, which is provided as a global panel at the provincial level, was aggregated to the country-level for the purposes of our estimation. Our full sample considers 59 Advanced Economy (AE) and Emerging Market (EM) countries. 18 We choose to exclude developing nations from this first sample for two reasons. First, data quality and availability in these nations is poor over our sample period. Additionally, many of these developing nations experienced lags in the contagion of the virus, thus hitting their peaks later than global hubs like Europe and the United States. 19 The central covariates of interest were pulled from Oxford's COVID-19 Government Response Tracker. As of April 29th, 2020, Oxford provides country-level indicators on containment and closure, economic responses; the quality of health systems; and unorthodox responses to the COVID-19 pandemic. 20 Of the 59 countries in our sample, government response data is reported for 55. We focus our estimation on the Stringency Index (SI), which mainly captures variation in government policies related to containment and closure. Each nation is scaled with a composite score from 0 to 100, with higher scores indicating more stringent policy interventions. In our dynamic panel estimation, we lag these interventions by between 2 and 5 weeks to account for delayed implementation or latent effects. In addition to the Oxford data, we control for various country-specific features using a wide range of publicly available data. We pull country coordinates from Google and integrate data on recent mobility trends from Apple. 21 From the World Bank, we gather World Development Indicators for the proportion of the population above age 65; the proportion of the population which is urban; total population density (people per 100 sq. km. of land area); tourist arrivals and departures; the proportion of vulnerable employment; gross national income per capita (calculated in current USD using the Atlas method); the number of cellular subscriptions; current health expenditures; and micrograms of PM2.5 air pollution per capita. 22 We also include cross-country indicators of the strength of democracy, from the Economist and Freedom House (2020). 23 From the Correlates of War Project, we aggregate data on military expenditures and personnel, iron and steel production, and energy consumption. 24 Lastly, we collect data on prior infections and deaths by the disease through the World Health Organization's International Classifications of Diseases. 25 This data is collapsed by country, across the years 2015-2018. The analysis is restricted to deaths by respiratory, endocrine, or high blood pressure conditions. We choose not to incorporate testing covariates, given discrepancies, and 19 With the notable exception of Iran (63 days), developing countries excluded from the sample had all recorded their first one hundred cases less than 50 days since the end of the sample on April 28th. Only eleven out of sample developing nations had recorded more than one hundred deaths by the end of the sample period. Future research will integrate developing nations into a larger global sample, as data availability and quality become more robust in later phases of the pandemic. 20 A detailed review of Oxford's dataset may be found at https://www.bsg.ox.ac.uk/sites/default/files/2020-04 /BSG-WP-2020-032-v5.0_0.pdf. 21 The data from Google can be found at https://developers.google.com/public-data/docs/canonical/countries_ csv; while the data from Apple was pulled from Mobility, apple.com, Apple (2020) https://www.apple. com/COVID19/mobility. The Apple data is calculated at a base of 100, where reduced daily mobility results in a lower score (<100), while higher mobility results in a higher score (>100). Mobility is measured across walking, driving, and public transit. 22 All of these data sources can be found through the World Development Indicators, World Bank (2020), https://data.worldbank.org/indicator. 23 Data was pulled from the Economist Intelligence Unit (2019); and from Freedom House's (2020) Freedom in the World scores. 24 Data was pulled from the National Material Capabilities (2017) index. 25 Data from the International Classification of Diseases (V10) can be found at https://www.who. int/classifications/icd/icdonlineversions/en/. measurement error in currently available data. 26 Depending on the quality of future data, we may choose to integrate these covariates in our future analyses. We take a multi-faceted approach to understand the cross-country dynamics of pandemic diffusion. As mentioned, the dynamic panel enables us to study to what degree stringency policies are associated with weekly mortality growth rates, and whether such policies start taking effect with a lag. We include the Stringency Index at date t and one-day lagged mortality growth (date t − 1), and then estimate a set of specifications with the independent variable being weekly mortality growth two to four weeks from date t via LSDV. 27 Given that residuals within countries are likely correlated, we employ robust standard errors clustered by country. Then, moving to cross-country analysis we focus on explaining differences in mortality outcomes at the peak of the first quasi-bell curve. For this, we rely on OLS estimation, adjusting the standard errors for heteroscedasticity. Finally, we conduct a cross-country survival analysis to better understand what drives variation in pandemic duration to the first peak (PD) across countries. Initially, we employ a cross-country Kaplan-Meyer analysis on the number of days to the first new mortality peak, subsequently extending it under a Cox proportional hazard model to introduce additional controls, estimated via maximum (partial) likelihood. For robustness, we also test and report whether the proportional hazards assumption is satisfied. Table 1 reports the baseline results from the local projection regressions (dynamic panel analysis). Notice that across local projection regressions at horizons two to four weeks, Stringency is statistically significant and negative. 28 These cursory results suggest that countries administering more stringent pandemic policies realized, on average, significantly lower mortality growth rates, with a lag of two to four weeks (Fig. 6) . The estimated effect sizes are also economically significant: countries with a Stringency Index (SI) 10 units higher than average in week W realized mortality growth rates that were on average 22, 17, and 13 percentage points lower in weeks W + 2, W + 3 and W + 4, respectively. While the reductions in growth rates seem quite large (and certainly should not be ignored), it is important to put these numbers in perspective. Given the exponential nature of disease spread, weekly mortality growth rates are highly volatile, potentially reaching anywhere from +50% to +100% (or greater). In the baseline analysis (Table 1) we do not explore cross-country heterogeneity in mortality dynamics. However, it is of great interest to understand whether and to what degree stringency policies were more effective in some countries than others in slowing down mortality growth. In the Appendix, we report local projection regression results which allow 26 Testing data was originally pulled from the Our World in Data (https://ourworldindata.org/coronavirus-testing) project, but is excluded in our first round of analysis. 27 Estimates from the first specification, LSDV estimation of a model where mortality growth is lagged by 1 week with time/country fixed effects are potentially biased (Nickell 1981 ) under small T. However, given the large time dimension of our panel, the LSDV estimator performs comparatively well to bias-corrected approaches (Judson and Owen 1999) . 28 The significance disappears at horizons more than 30 days, so we only report results at horizons two to four weeks. for the impact of SI to depend on each of our country fundamentals following Eq. 3. We take these country fundamentals as exogenous, as they are from the period pre-COVID-19 and are highly persistent (e.g., population, democracy). Due to the number of country fundamentals we are interested in exploring, we estimate regressions one-by-one. While this yields parsimonious interpretations and preserves degrees of freedom, a key drawback of this approach is that it may suffer from omitted variables. That the estimates of policy stringency became statistically insignificant when its interaction with some of the country fundamentals may suggest that the level of stringency applied was correlated with several pre-conditioning variables in the first-wave of COVID-19. Figure 7 characterizes the results into impulse response plots, highlighting the impact of a 10 unit rise in SI for countries in the 25th percentile of a given characteristic against those in the 75th percentile of the same characteristic. For example, notice that countries in the 75th percentile in terms of proportion elderly (65 or older) saw a much stronger reduction in mortality growth rates from the same 10 unit rise in SI, compared to countries with relatively low proportions of elderly (25th percentile): Countries in the 75th (25th) percentile saw mortality growth rates fall about −25% (−10%) two weeks later. These results are consistent with the proportion of the elderly being a risk factor for the pandemic. Countries may be more susceptible to the pandemic conditional on a greater proportion of the elderly, and for this same reason, stringency policies may be more beneficial in these countries. We also find that in countries with a higher proportion of urban population, SI measures had a stronger impact on mortality growth than in countries with a lower proportion of urban population. This heterogeneity is be consistent with the proportion of urban population being a risk factor for the pandemic as urban population may face greater exposure possibility. SI policies are also more strongly associated with lower mortality growth in countries with greater international travel arrivals, greater GNI per capita, higher health expenditures, and greater EIU Democracy. While greater international flow is quite an intuitive risk factor for a pandemic like COVID-19, the role of democratic freedom is an ongoing topic of debate. Our results may support the view that greater ex ante individual liberties may be detrimental in the present situation in the sense that they increase the risk of pandemic severity (thereby leading to a stronger effect of mandated SI policies). 29 While our results suggest that stringency policies were more impactful among countries that are considered more democratically free, the effectiveness of governments at responding to Covid-19 also depends on other complementary factors, notably leadership and administrative capacity (Ang 2020) . The specific channel the association captures is not clear at this point, but it's possible that countries with higher EIU Democracy scores also exhibit greater ex-ante mobility. Alternatively, it's possible that there is greater variance in publicly available information among high EIU countries, inducing differential beliefs related to the seriousness of the pandemic. Both 29 Civil liberties is one of the five categories based on which the EIU Democracy Index is calculated. Average Impact Fig. 6 The mortality rate is negatively associated with the intensity of government response. Note: Pooled estimates from local projections are represented as gray circles. Error bars reflect 95% confidence intervals based on HAC-robust standard errors clustered by country scenarios can potentially lead to greater pandemic risk, and therefore greater stringency policy effectiveness. Kaplan-Meyer Curves of Pandemic Duration-to-Peak across Countries Figure 8 reports Kaplan-Meyer (KP) curves on PD (pandemic duration to the first peak). The Y-axis can be interpreted as the 'probability of peaking later than day (t)' or equivalently, as the 'probability the first peak is yet to come by day (t).' The LHS plots the unconditional duration probability and shows how the probability the peak is yet to come is decreasing with the number of days since the first death increases (X-axis). About 35 days into the pandemic, the probability of peaking reaches 50%; 45 days into the pandemic, the probability rises to 75%. The center chart plots the KP curves of countries stratified by early SI policy (above versus below average). According to the KP curve, pandemic durations were significantly shorter in countries which proactively issued more stringent policies. That is, for high Early SI countries, the probability of peaking reached 25% by about 35 days while it took close to 50 days for countries with low Early SI policies. The probability of peaking after 40 days is close to 0% for high Early SI countries (virtually all of these countries experienced their first local peak by Not only is there wide variation in mortality dynamics, but we also observe significant heterogeneity in mortality peaks and curvatures across countries. Cross-section regression results, where the dependent variable is the logged peak new mortality rate (calculated as the deaths/population at the first peak of daily new deaths) are reported in Table 2 , column [1]. In addition, we include results for the outcome logged peak new mortality rate to PD (column [2]), and results from the survival analysis via Cox regression on PD itself (column [3]). Several broad patterns emerge. First, early stringency measures are significantly negatively associated with mortality peaks: Greater early stringency is associated with lower peak mortality rates (column [1], an estimate of −0.12), implying that countries with Early SI 1 unit higher than average realized first mortality peaks that were about −12% lower. Quantitatively, greater policy stringency early on was also associated with flatter mortality curves (column [2]), but the estimate is statistically insignificant. However, there is a significant negative association between logged peak mortality rates and the PD, which implies that more aggressive policy interventions early on may still be associated with a shorter PD through its negative impacts on mortality peaks. Additionally, we find no evidence suggesting that the timing of initial policy interventions is correlated with the peak or slope of mortality curves (columns [1] and [2]). Unsurprisingly, higher early mortality growth rates are significantly associated with higher peak mortalities and steeper mortality curves (columns [1] and [2]), while higher levels of the mortality rate early on were associated with flatter overall mortality curves, but longer pandemic durations to the first peak (columns [2] and [3]). Greater early mobility was significantly associated with higher mortality peaks (columns [1]) and faster response, that is a larger number of days between the first non-zero SI and the first death, was significantly associated with shorter pandemic duration (column [3] ). Countries with greater elderly populations, greater population density, and/or greater shares employed in vulnerable sectors saw, on average, significantly longer PDs (column [3]). At the same time, countries scoring higher in EIU Democracy and Fig. 8 Time-to-peak duration analysis of mortality rates. Note: Y-axis indicates the probability the peak mortality/case is 'yet to come'. The higher y-axis implies a lower probability of peaking. X-axis reflects the number of days since the first mortality/case was realized. Shaded areas represent 95% confidence intervals those further from the equator realized shorter PDs. Higher population density is, somewhat surprisingly, associated with flatter mortality curves (column [2]). However, this association may be driven by the fact that many of the high-density countries are in Asia. These countries Table 2 Explaining cross-country differences in the empirical shape of mortality rate curve from the onset to the local peak of the first pandemic phase contained the spread of COVID-19 relatively effectively, given their preparedness in light of battling SARS in 2003. We find no cross-sectional evidence of increasing SI policies during the pandemic (Stringency Delta) affecting peak mortalities or PD. The cross-sectional regression and survival analysis results, taken together with the evidence from local projections in our dynamic panel analysis, indicate (but do not assert) that 1) aggressive policy responses may have helped slow down the growth rate of new mortalities with a lag, 2) having stringent policies in place early may have helped to lower peak mortality rates and PD (through reducing peak mortality rates). In contrast, evidence on the effectiveness of increasing the SI after the fact (via Stringency Delta) on peak mortality is much weaker. To summarize, we investigate the impact of stringency policies on mortality growth dynamics, along with cross-country patterns in the empirical shape of the mortality curve. First, we find that higher SI levels were significantly associated with lower mortality growth rates with a lag of 2, 3, and 4 weeks. The effect of a 10-unit rise in the SI levels on future mortality growth was stronger in countries that appear to be more vulnerable ex-ante to COVID-19 type breakouts: countries with greater elderly populations, cooler temperatures, more international travel flow, and higher levels of the EIU Democracy index. In terms of peak mortality rates across countries, we find that proactive stringency policies (Early SI)those already in place prior to the first reported mortality -were significantly associated with lower mortality peaks, while higher mobility levels were significantly associated with higher mortality peaks. As for pandemic duration-to-first peak, higher mobility levels were associated with shorter average pandemic durations (despite higher mortality peaks), while greater elderly population and greater shares of employment in vulnerable sectors were associated with longer average pandemic durations. Estimates from a Kaplan-Meyer analysis indicate that proactive stringency policies early on were associated with significantly shorter PD. We have conducted several robustness checks. First, to explore the potential nonlinear effect of SI on weekly new mortality growth, we replace the level variable of SI in the local projection method by the logarithm transformation of SI. Results are similar to our main findings: there were significantly negative associations between weekly new mortality growth and the log(SI) at horizons two to four weeks. Specifically, a 10 % increase in SI was associated with 7.1, 5.8, and 4.7 percentage points decrease in weekly new mortality growth rates, at horizon two, three and four weeks, respectively. These results suggest nonlinearity in the effects of SI on weekly new mortality growth rates: an increase in SI when the SI level is low is likely to have a greater effect on reducing mortality growth rates. Such nonlinearity is also consistent with the finding in the cross-sectional analysis that countries had a greater average SI before the first death saw a lower new mortality rate at the peak. Second, in addition to closure and containment measures, the SI constructed by Oxford COVID-19 Government Response Tracker considers public information campaigns, which may be argued as different from closure and containment measures. Hence, we check this by constructing an alternative SI excluding the public information campaigns. Additionally, we construct a second alternative SI which only consists of domestic closure and containment measures, 30 as many countries implemented international travel restrictions early and aggressively while not carrying out domestic closure and containment policies. Results using these alternative SIs are similar to our main findings: alternative measures of stringency were 30 Specifically, indicators with ID C1-C7 in Oxford COVID-19 Government Response Tracker. significantly and negatively associated with weekly new mortality growth, at horizons two to four weeks. Third, we experiment with including lagged SI as controls and find similar results. 31 In the Appendix, we also report a robustness check using growth in cumulative mortalities (as opposed to growth in new mortalities) and a first pass at a residual analysis based on our benchmark estimation, identifying countries which appear to be statistical outliers. 32 There are countries oft-cited in the public discussion as exceptionally better or worse than predicted in terms of their infection and mortality rates, several of which are consistent with the residuals in our estimation. The contributing factors are likely to be many, including that several of these countries previously experienced recent outbreaks, including SARS-CoV-1 (2002) (2003) (2004) . These include: (i) mortality rates worse than expected: Brazil, Ecuador, Iran, Italy, Peru, United States; (ii) mortality rates better than expected: Cambodia, Greece, India, Indonesia, Iraq, Japan, Lebanon, Myanmar, South Korea, Sub-Saharan Africa (including Ebola-affected areas), Thailand, Venezuela, Vietnam. Japan, for example, offers an interesting case study with its proximity to the initial hotspots (China and South Korea), but is currently cited as successfully managing the infection and mortality rates without lockdown or mass testing (Bloomberg 2020) . Additionally, the degree of enforcement of policy intervention rules may be of great importance in explaining the cross-country difference in the empirical shape of the mortality curve, which, however, lacks data to measure. As the pandemic continues to unfold and as new data comes in, we'll continue the process of analyzing these exceptional cases for country-specific factors underlying their mortality rates from COVID-19. Avery et al. (2020) provide an overview of the modeling of the spread of COVID-19. By and large, the ongoing challenges surround data on the infection rates. As noted by Manski and Molinari (2020) , because of missing data on tests for infection and imperfect accuracy of tests, reported rates of population infection by the SARS CoV-2 virus are lower than actual rates of infection, resulting in infection fatality rates that are lower than reported. In addition, as argued by Atkeson (2020) , in the presence of effective mitigation measures, the model with a high initial number of active cases and a low fatality rate gives the same predictions for the evolution of the number of deaths in the early stages of the pandemic as the same model with a low initial number of active cases and a high fatality rate. Modeling the COVID-19 infection also face with uncertainty about the evolution of infection rates, due to parameter uncertainty and the realization of future shocks (Liu et al. 2020) . To ameliorate these empirical challenges, our study uses the mortality rate per capita as the main variable of interest. In this section, we synthesize our estimates with findings from related studies in the field. Notwithstanding different empirical specifications across studies, we find our main results are consistent and communicable with the literature and offer valuable new evidence. 31 All results of robustness checks are shown in the Appendix. 32 Results of the panel analysis on cumulative mortality growth rates are broadly consistent with that on new mortality growth: Greater policy stringency is associated with lower cumulative growth, with a lag, and the association is more persistent (significant up to five weeks later). Also, we observe heterogeneity in the associations: Such associations are more pronounced in countries with a greater proportion of elderly population, a greater proportion of urban population, further away from equator, a higher level of GNI per capita, a greater share of health expenditure, and a higher level of democracy. Closely related to our study are The Economist (2020), Stojkoski et al. (2020) and Castex et al. (2020) . The Economist's review uses a sample of U.S. states focusing on the case-fatality rate. It finds the following associations between the case-fatality rate and the listed variables, ordered by their significance (one standard deviation increase): median age (pos.), ICU beds per 100,000 people (neg.), population density (pos.), and median income (neg.), the prevalence of heart disease, diabetes and smoking (pos.), the share of the population that is African-American (pos.), and amount of social distancing three weeks prior (neg.). Stojkoski et al.' s study finds that greater government intervention, measured by a slight variant of the stringency index used in our research, is associated lower mortality per capita. Our results are consistent with this. Given our most comprehensive specification in the cross-country estimates, we find a significant association between the stringency index and the mortality rate peak levels ( Table 2, column [1] ). An increase of 1 unit of the stringency index is statistically associated with −12% decrease in the first peak mortality rate (Table 2, column [1]). Our cross-country results revealing heterogeneity effects of government interventions on COVID-19 mortality growth rates also relate to those findings of Castex et al. (2020) on the COVID-19 transmission rates; that the effectiveness of non-pharmaceutical interventions (such as school and workplace closures) is increasing with lower air pollution and higher health expenditure. They also find that the effectiveness of such interventions is declining with higher proportion of elderly population, which seems to contradict with our results. One potential reason could be the dependent variables are different, with us focusing on mortality rates and them focusing on transmission rates. 33 Our study finds supportive evidence for the role of the aging population, urbanization, preexisting conditions, mortality from high-blood pressure, obesity, diabetes, and trust. The findings are consistent with case studies on the association between social networks and COVID-19 infection. Kuchler et al. (2020) find in Facebook data that areas with stronger social ties to two early COVID-19 "hotspots" (Westchester County, NY, in the U.S. and Lodi province in Italy) generally have more confirmed COVID-19 cases, after controlling for geographic distance to the hotspots as well as for the income and population density of the regions. Allcott et al. (2020) find in location data from a large sample of smartphones to show that areas with more Republicans engage in less social distancing (controlling for other factors including state policies, population density, and local COVID cases and deaths, pointing to significant gaps between Republicans and Democrats in beliefs about personal risk and the future path of the pandemic. Fetzer et al. (2020) find that the perception of a weak government and public response is associated with higher levels of worries and depression, and those strong government reactions correct misperceptions and reduce worries and depression. Our findings on the association between trust, the stringency index, and mortality from COVID-19 yield consistent messages along this line. We find that the mortality growth rates are associated with the impact of stringency policies interacted with distance from the equator: countries in higher latitude appear to benefit more from stringency policies, as they are associated with lower mortality growth from an equivalent rise in stringency relative to countries which lie closer to the equator. This evidence suggests that geographic location may be a potential COVID-related risk factor such that stringency measures are more effective in higher risk, high latitude countries. The finding is consistent with evidence linking temperatures to the spread of influenza. Slusky and Zeckhauser (2019) find that sunlight strongly protects against influenza, a relationship driven by sunlight in late summer, and early fall (when there are sufficient quantities of both sunlight and influenza activity). Across countries, our estimates suggest that government policy, institutions, and the intensity of government response to COVID-19 are negatively associated with the mortality per capita. This evidence reflects in the coefficient estimates of stringency index, government effectiveness, democracy, health expenditures, vulnerable employees, capita GNI, and level pollution. These cross-country findings are consistent with several case studies of COVID-19. Dave et al. (2020) find in the data of U.S. states that approximately three weeks following the adoption of a Shelter-in-Place Orders (SIPO), cumulative COVID-19 cases fell by 44% and, in an eventstudy analysis, that SIPO-induced case reductions grew larger over time, with the early adopters and high population density states appear to reap larger benefits from their SIPOs, though the estimated mortality reduction effects were imprecisely estimated. This finding is consistent with our estimates of the stringency index under the local projection approach. Our finding in the cross-country data on the importance of the extent of vulnerable employees is largely consistent with the evidence from influenza. Markowitz et al. (2019) find in the U.S. data that a one percentage point increase in the employment rate increases the number of influenza-related doctor visits by about 16%; these effects are highly pronounced in the retail sector and healthcare sector, the sectors with the highest levels of interpersonal contact. Clay et al. (2018) find evidence on the link between air pollution and influenza infection and suggest that poor air quality was an important cause of mortality during the pandemic. The empirical issues remain on the endogeneity of mobility and in-bound/out-bound travels to government response in the estimation. For instance, Gupta et al. (2020) find that mobility fell substantially in all the U.S. states, even ones that have not adopted major distancing mandates. They find that there is little evidence, for example, that stay-at-home mandates induced distancing; in contrast, early and information-focused actions have had bigger effects. Their event studies show that first case announcements, emergency declarations, and school closures reduced mobility by 1-5% after 5 days and 7-45% after 20 days. We are in the process of addressing endogenous regressors in our estimation (policy stringency, lagged mortality, and mobility) with more updated data. We conclude by cautioning that our results are subject to data limitations, including undercounts of COVID-19 infections and mortality. 'Better or worse performance' of a country in the first phase of the pandemic does not guarantee similar future outcomes. The levels of COVID-19 policy stringency applied in the first wave were correlated with several pre-conditions [aging population, urbanization, travels, income, health spending, and democracy], the correlations that might influence the policy efficiency in the following wave(s) if any (Table 3) . Flattening the mortality and infection curves may shift mortality and painful adjustment forwards. Premature opening of the economy without proper testing, contact-tracing, and selective quarantines of vulnerable or impacted segments of the population may induce future acceleration of the pandemic (Acemoglu et al. 2020 ). More medical research and advances towards better treatment and possible vaccinations, the quality of local and global public policies, and adjustment capabilities of countries will determine future dynamics of the pandemic (Lipsitch et al. 2020) . Average Impact Fig. 9 The cumulative mortality growth rate is negatively associated with the intensity of government response. Note: Results of panel analysis on cumulative mortality growth rates. Pooled estimates from local projections are represented as gray circles. Error bars reflect 95% confidence intervals based on HAC-robust standard errors clustered by country A multi-risk SIR model with optimally targeted Lockdownm Polarization and public health: partisan differences in social distancing during the coronavirus pandemic When COVID-19 meets centralized, personalized power How deadly is COVID-19? Understanding the difficulties with estimation of its fatality rate Policy implications of models of the spread of coronavirus: perspectives and opportunities for economists Did Japan just beat the virus without lockdowns or mass testing? COVID-19: cross-country heterogeneity in effectiveness of nonpharmaceutical interventions Pollution, infectious disease, and mortality: evidence from the 1918 Spanish influenza pandemic Regression models and life-tables When do shelter-in-place orders fight COVID-19 best? Policy heterogeneity across states and adoption time Global behaviors and perceptions at the onset of the COVID-19 pandemic Proportional hazards tests and diagnostics based on weighted residuals Survival analysis and regression models Tracking public and private responses to the COVID-19 epidemic: evidence from state and local government actions International Monetary Fund (2020) Policy responses to COVID19 Estimation and inference of impulse responses by local projections Estimating dynamic panel data models: a guide for macroeconomists The Geographic Spread of COVID-19 Correlates with Structure of Social Networks as Measured by Facebook Defining the epidemiology of COVID-19-studies needed Panel forecasts of country-level Covid-19 infections The incubation period of COVID-19: a rapid systematic review and meta-analysis of observational research Estimating the COVID-19 infection rate: anatomy of an inference problem The effects of employment on influenza rates National Material Capabilities (2017) Correlatesofwar.org, Correlates of War Biases in dynamic models with fixed effects Government spending multipliers in good times and in bad: evidence from US historical data Sunlight and protection against influenza The socio-economic determinants of the coronavirus disease Democracy index -world democracy report The south is likely to have America's highest death rate from COVID-19 Using social and behavioural science to support COVID-19 pandemic response Estimates of the severity of coronavirus disease 2019: a model-based analysis Acknowledgments We gratefully acknowledge the insightful comments by two anonymous referees, and the financial support by the Dockson Chair and the Center of International Studies at USC.Data Availability and Code Availability Data and codes available at Github Repository: https://github. com/snairdesai/COVID-19. Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.