Factors influencing the COVID-19 daily deaths peak across European countries ============================================================================ * Katarzyna Jabłońska * Samuel Aballéa * Mondher Toumi ## Abstract **OBJECTIVES** The purpose of this study was to determine predictors of the height of COVID- 19 daily deaths peak and time to the peak, in order to explain their variability across European countries. **STUDY DESIGN** For 34 European countries, publicly available data were collected on daily numbers of COVID-19 deaths, population size, healthcare capacity, government restrictions and their timing, tourism and change in mobility during the pandemic. **METHODS** Univariate and multivariate generalised linear models using different selection algorithms (forward, backward, stepwise and genetic algorithm) were analysed with height of COVID-19 daily deaths peak and time to the peak as dependent variables. **RESULTS** The proportion of the population living in urban areas, mobility at the day of first reported death and number of infections when borders were closed were assessed as significant predictors of the height of COVID-19 daily deaths peak. Testing the model with variety of selection algorithms provided consistent results. Total hospital bed capacity, population size, number of foreign travellers and day of border closure, were found as significant predictors of time to COVID-19 daily deaths peak. **CONCLUSIONS** Our analysis demonstrated that countries with higher proportions of the population living in urban areas, with lower reduction in mobility at the beginning of the pandemic, and countries which closed borders having more infected people experienced higher peak of COVID-19 deaths. Greater bed capacity, bigger population size and later border closure could result in delaying time to reach the deaths peak, whereas a high number of foreign travellers could accelerate it. Keywords * COVID-19 * mortality * healthcare capacity * modelling ## Introduction The coronavirus infectious disease 2019 (COVID-19) outbreak was announced as a pandemic by the World Health Organization (WHO) on 11 March 2020.1 By the end of March, Europe exceeded Asia and became the region experiencing the highest percentage mortality from the virus across the world,2 until 10 June 2020, when the rapidly growing COVID-19 mortality rate in the Americas exceeded all other continents.3 According to the WHO report from 5 July 2020, 38% of the global mortality due to COVID-19 was from Europe.4 Since incidence and mortality rates varied between countries, numerous studies have recently been published investigating factors associated with COVID-19 infection and death rate across countries. A variety of potential predictors have been assessed in the literature, such as country-specific demographic and health characteristics, economic and social indicators,5-8 mobility scores and social-distancing measures,9, 10 as well as ecological and environmental perspectives.8, 11 To assess the relationship between covariates and COVID-19 incidence or mortality, the most common approach used by authors was to analyse multivariate regression models of total number of infections or deaths up to a given time point, using log-transformed data or without any transformation, as well as daily data on infections or deaths as outcomes. In this study, we used data on numbers of deaths, not infections, since the former has a much higher degree of reliability than the latter, being better monitored and less dependent on the number of tests done. Since all European countries seem to already reach the peak of deaths by 3 June 2020 from the first wave of COVID-19 disease, our idea was to use height of daily deaths peak as a primary outcome of interest, and time to the peak as a secondary outcome. To the best of our knowledge, this perspective has not been investigated so far. When raw and cumulative daily numbers of infections and deaths are subject to deviations from between-country differences in reporting and depend on the date up to which the analysis is performed, maximum number of daily deaths can be assessed as an interesting new indicator of disease mortality magnitude. In addition, analysis of the height of daily deaths peak enables us to assess the overall capacity of healthcare systems. This study aims to detect significant drivers of COVID-19 mortality with the use of multivariate generalised linear models (GLM) and distinct selection algorithms, to explain the variability of height of and time to the deaths peak among European countries. This will enable us to draw conclusions about how governments and societies can improve the future response on similar pandemic or global life-threatening situations. ## Methods ### Data collection A total of 34 European countries were included in the analysis. The height of COVID-19 daily deaths peak was the primary outcome of interest, defined as maximum daily reported number of people who died due to COVID-19 per country up to 3 June 2020. The value was then divided by the number of inhabitants of a given country reported in January 2020 and multiplied by 1 million to obtain the deaths peak height per 1 million inhabitants. The secondary outcome of interest was time to COVID-19 daily deaths peak, defined as number of days from the day when the first death was reported in a given country, up to the day of reaching the COVID-19 daily deaths peak. A set of explanatory variables used to predict outcomes consisted of the following: * Healthcare capacity * ○ All beds capacity (number of hospital bed units) * ○ Intensive care unit (ICU) beds capacity (number of ICU) * ○ Number of tests conducted up to the time of the peak * Government restrictions and associated factors * ○ ‘Stay at home’ order date * ○ Educational facilities closure date * ○ Gathering restriction date * ○ Businesses closure date * ○ Border closure date * ○ Total number of COVID-19 infected people when borders were closed * ○ Total number of COVID-19 deaths when borders were closed * Indicators of the population size * ○ Country population size (in January 2020) * ○ Percentage of population living in urban areas * ○ Percentage of population living in metropolitan cities with more than 1 million inhabitants * Median age * Tourism: * ○ Number of travellers that arrived at airports in 2018 * ○ Number of foreign tourists in 2018 (arrivals at any touristic accommodation) * Mean mobility score at the day of first reported COVID-19 death (as a difference from the average mobility reported before the pandemic), calculated across mobility scores at * ○ Retail and recreation places * ○ Workplaces * ○ Transit stations All dates were considered as number of days after the day of first reported COVID-19 death in each country. Variables indicating healthcare capacity and tourism were considered in relation to the population size of a given country (per inhabitant or 1 million inhabitants). #### Data sources Data on COVID-19 deaths, infections, number of tests, all and ICU beds capacity, dates of government restrictions, population size and urban population size were taken from Worldometer.com.12 Missing dates of government restrictions, if officially issued, were found on Wikipedia.13 Mobility scores were uploaded from Google COVID-19 Community Mobility Reports.14 Scores were reported as percentage changes from a usual mobility calculated before the pandemic in Europe, between 3 January and 6 February 2020. Mobility scores were subject to oscillations due to daily reporting and presence of weekends and holidays; therefore, we used smoothed mobility scores produced with a nonparametric technique that uses local weighted regression and fits a smooth curve through points in a scatter plot, called *loess regression*.15 Data on number of passengers arrived at airports, tourism and population living in metropolitan areas were downloaded from Eurostat,16 except for Russia, Turkey and Ukraine, for which data were gathered from other sources. 17-20 #### Missing data and imputation If government restrictions were not officially set, the date was imputed with date of the peak in each country. Missing data on mobility for Cyprus and Iceland were imputed with average scores of remaining countries. ### Statistical analysis Descriptive statistics on outcomes and explanatory variables were produced. Since the distribution of deaths peak height per population size is right-skewed, a logarithmic transformation was applied. Factors influencing the height of deaths peak were analysed using univariate and multivariate GLM with a normal distribution function and logarithmic link function. The base case scenario analysis was performed using data available for 34 countries and explanatory variables with a p-value <0.1 in univariate GLM models. To avoid multicollinearity issues between independent variables, Pearson correlations were investigated to detect highly correlated pairs of variables. The correlation was assessed as high if the absolute value exceeded 0.7 and as moderate if the absolute value was contained in a range 0.5–0.7.21 Due to relatively low sample size, a risk of bias could appear for models with too high number of variables included. To overcome the problem and limit the number of covariates, sensitivity analyses were performed using variety of selection algorithms: stepwise, backward, forward and the genetic algorithm. For the backward, forward and stepwise algorithms, a criterion of having p-value lower than 0.1 was applied for each variable to stay in a model (backward and stepwise) and to enter a model (forward and stepwise). For the genetic algorithm, the best model was fitted based on the value of Akaike’s Information Criterion (AIC) corrected for small sample sizes (AICC). Models with only main effects were considered. Additional sensitivity analysis was performed removing countries for which imputation of government restriction dates was needed. Similar methods were used to analyse time to deaths peak but without any transformation since it seemed to follow a normal distribution. GLM models with a normal distribution function and identity link function were analysed. For all analyses, a p-value lower than 0.05 was considered as statistically significant. For each GLM model, fit statistics such as AIC and its equivalent AICC were produced with lower values indicating better fit. Analyses were performed using SAS 9.4 software. R 3.6.2 software was used to apply the genetic algorithm with the package *glmulti*.22 ## Results ### Descriptive analysis Characteristics of the countries included in the study are presented in Table 1, and a histogram of the height of COVID-19 deaths peak is depicted in Figure 1. The graph presenting deaths peak height by country is provided in Figure 2. View this table: [Table 1.](http://medrxiv.org/content/early/2020/11/05/2020.11.04.20225656/T1) Table 1. Characteristics of countries included in the analysis. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/05/2020.11.04.20225656/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/11/05/2020.11.04.20225656/F1) Figure 1. Histogram of the height of COVID-19 daily deaths peak per 1 million inhabitants with a fitted log-normal curve Histogram presents the distribution of height of COVID-19 deaths peak per 1 million inhabitants across 34 countries. A log-normal curve was fitted to the distribution with scale parameter sigma=1.12 and location parameter zeta=1.38, assuming logarithm of height of the deaths peak is normally distributed [Normal(zeta, sigma)]. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/05/2020.11.04.20225656/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/11/05/2020.11.04.20225656/F2) Figure 2. Height of COVID-19 daily deaths peak per 1 million inhabitants across countries The plot presents the height of COVID-19 deaths peak per 1 million inhabitants across 34 countries, with Belgium having outstandingly highest peak per population size (>40 deaths per 1 million inhabitants). Median height of the peak per 1 million inhabitants per day was equal to 3.48 deaths ([lower quartile; upper quartile] = [1.68; 12.78]), with Belgium reaching outstandingly higher peak than other countries. As can be seen from these statistics and the histogram, the height of the deaths peak does not follow a normal distribution, but it seems it can be well approximated using a log-normal distribution. Median time to the peak equalled 31 days from the first death reported in each country with a comparable mean (31.32, SD=13.94). Average numbers of days when closing schools/universities and when gatherings were banned were both negative (−0.26 and −1.76, respectively). ### Height of the deaths peak Univariate GLM models of height of the peak were analysed (Table 2). Eight factors turned out to be significant: educational facilities closure day, gathering restrictions day, businesses closure day, proportion living in urban areas, mobility score at the day of first reported death, border closure day, number of COVID-19 infections when borders were closed, and number of COVID-19 deaths when borders were closed. All of them have positive estimates. Stay-at-home order day, as well as the proportion living in metropolitan cities >1 million inhabitants were close to reaching the significance level (p<0.1). Pearson correlations between variables reaching or close to reaching significance in univariate models were verified. Results can be found in Supplementary Materials. View this table: [Table 2.](http://medrxiv.org/content/early/2020/11/05/2020.11.04.20225656/T2) Table 2. Results of univariate and multivariate GLM of COVID-19 deaths peak height, with normal distribution and logit link function A base case multivariate GLM model included six covariates. The proportion living in urban areas was found significant (p<0.001) and the mobility score at the day of first reported death was very close to reaching significance (p=0.052). Results are presented in Table 2. Considering small sample size (N=34), selection algorithms were applied to the base case model to limit the number of covariates, increase model precision and improve the model fit. All selection algorithms were consistent and selected the model with three significant parameters (“final” model), indicating its best fit. The proportion of the population living in urban areas (6.848, p<0.001), mobility score at the day of first reported death (0.049 p<0.001) and number of infections when borders were closed per 1 million inhabitants (0.0002, p=0.016) were all significantly associated with the deaths peak height (Table 2). Analysis of data removing countries with imputed dates of businesses or borders closure (N=29) provided consistent results regarding the significance of all three covariates from the final models in relation to the deaths peak height (details are provided in Supplementary Materials). ### Time to the deaths peak Univariate GLM models of time to the peak were analysed (Table 3). Nine factors turned out to be potentially significant (with p<0.1): all beds capacity, ICU beds capacity, educational facilities closure day, gathering restrictions day, population size, arrivals at airports per inhabitant, number of foreign tourists per inhabitant, mobility score at the day of first reported death, as well as day of border closure. The majority of these factors have positive estimates, except the number of arrivals and tourists. As for the next step, linear correlations between them were verified. Results can be found in Supplementary Materials. View this table: [Table 3.](http://medrxiv.org/content/early/2020/11/05/2020.11.04.20225656/T3) Table 3. Results of univariate and multivariate GLM of time to COVID-19 deaths peak, with normal distribution and identity link function A base case multivariate GLM model was run using six explanatory variables. All beds capacity, population size in millions and the number of foreign tourists per inhabitant were found significantly related with time to COVID-19 deaths peak. Results are presented in Table 3. Considering the small sample size (N=34), selection algorithms were applied to the above multivariate model. The “final” model had four covariates and was selected by the backward and the genetic algorithm (details in Supplementary Materials). All beds capacity per 1 million inhabitants (0.003, p=0.004), population size in millions (0.142, p=0.008), the number of foreign tourists per inhabitant (−2.651, p=0.037) and border closure day (0.297, p=0.008) were all significantly associated with time to the deaths peak (Table 3). Analysis of data removing countries with imputed dates of borders closure (N=31) provided consistent results regarding the significance of all beds capacity and population size in relation to time to deaths peak, whereas number of foreign tourists and borders closure day were close in reaching significance, both with p<0.1 (details provided in Supplementary Materials). ## Discussion This study provides some evidence about factors associated with COVID-19 mortality peak and time to the peak. One of the strongest predictors of the COVID-19 mortality peak identified in our analysis was the proportion of population living in urban areas. The relationship between urbanisation and COVID-19 infections ratio was earlier outlined by the United Nations Association,23 estimating that 90% of all reported COVID-19 cases (by July 2020) came from urban areas, becoming the epicentre of the pandemic. High population density increases the propensity of viruses spread by increasing the contact rates of individuals. However, some authors24 warn readers against putting too much weight on urban density, arguing that large cities just faced the coronavirus earlier (due to the higher number of incoming people) and that the timing of epidemy start was of bigger interest than population density itself. This observation seems to be on the contrary to our study, since neither proportion of population living in urban areas, nor proportion living in metropolitan cities were related with the time to deaths peak. Number of infections when borders were closed was found to be another important factor associated with the COVID-19 deaths peak height, whereas a positive association between borders closure day and time to reach the peak was observed in this study. It shows that stopping arrivals to the country at the earlier stage of epidemy can be crucial in reducing the peak height, stopping the increase of daily number of deaths earlier and, consequently, to flatten the deaths curve. These findings are on the contrary with Chaudhry who showed no association between rapid border closures and COVID-19 mortality per million people, using cumulative data available as of 1 April 2020.25 However, the peak height and mortality may not be strongly correlated. For example some countries experience a point peak and dropping fast, (France, United Kingdom, Spain, Italy) while other countries may experience a flattened peak with a mortality staying high over a long period of time (US, Brazil, Mexico).12 The difference in average mobility before the pandemic and at the day when first death was reported in a country was found to be associated with deaths peak height, which suggests that rapid reduction in mobility was important in terms of reducing overall mortality. These findings are consistent with other authors opinion that changes in social-distancing covariates can flatten the curve by changing the peak death rate. 9, 10 Tourism indicators were found not to be associated with deaths peak height, but to be significant drivers of time to reach the peak. It suggests that the magnitude of inbound tourism can have an impact on accelerating the moment of reaching the peak. However, Aldibasi5 and Garcia de Alcaniz8 found tourism to be a significant predictor of COVID-19 mortality and infection rates. Ostig and Askin26 also found a significant positive relationship between number of airline passengers and number of COVID-19 infections. Another conclusion resulting from the study is that timing of government restrictions, especially border closure, can be found as an important factor in terms of COVID-19 mortality. What can be observed for many countries is that closing borders (and other government restrictions as well) took place several days before any death was reported in a country. Therefore, increase or decrease in mortality can be seen as a result of undertaking preventive action by the government. There is no general conclusion on the association between COVID-19 fatality and hospital beds capacity in literature. Our study showed lack of association between beds capacity and deaths peak height, but higher beds capacity was related with longer time to peak. Garcia de Alcaniz showed no association between hospital beds density and both the number of infections and number of deaths at any moment of the pandemic.8 On the contrary, Sorci showed that case fatality rate was negatively associated with number of hospital beds per inhabitants.7 High bed capacity may be associated with more intensive treatment delaying time to death for a number of patients thus delaying time to peak but does not affect the peak height. ## Limitations Our study has several limitations. Data on a limited number of countries were used. We decided not to include countries for which reliable and comparative evidence could not be found, nor countries with very low mortality. Our idea was to focus only on Europe, not to merge it with countries from other continents and avoid data incomparability issues and eventual difference in strain as the virus mutate rapidly. However, the number of observations can be assessed as high enough to draw conclusions based on multivariate GLM models by using selection algorithms. The stability of results was tested with models based on different variable selection algorithms, being consistent with base case findings. The set of included explanatory variables can be viewed as non-exhaustive. We decided to focus on variables important from a social perspective, assuming that mobility factors, tourism, urbanisation, and government decisions are more of interest to inform upcoming decisions than other approaches observed in the literature; for example, ecological or environmental perspectives.8, 11 ## Conclusions This study demonstrated the significant drivers of COVID-19 mortality magnitude in Europe and shed a light on reasons behind the variability between different countries. Countries with a higher proportion of population living in urban areas, without a rapid reduction in mobility and with delayed government restrictions, correlated with a higher COVID-19 deaths peak. These findings can help improve future response in similar situations or in case of a second wave. ## Supporting information Supplementary Materials [[supplements/225656_file02.pdf]](pending:yes) ## Data Availability The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials. ## Contributors MT conceived the study, KJ collected and analysed the data, KJ, MT and SA interpreted the results, KJ wrote the first draft of the manuscript, KJ, MT and SA revised the manuscript and approved the final version. ## Competing interests The authors declare no competing interests. ## Acknowledgements The authors wish to thank Saima Khan (Apothecom) for proof reading and editing support. * Received November 4, 2020. * Revision received November 4, 2020. * Accepted November 5, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.World Health Organisation. Coronavirus disease 2019 (COVID-19) Situation Report – 51 (11 March 2020) 11 March 2020. 2. 2.World Health Organisation. Coronavirus disease 2019 (COVID-19) Situation Report – 71 (31 March 2020) 31 March 2020. 3. 3.World Health Organisation. Coronavirus disease (COVID-19) Situation Report – 142 (10 June 2020) 10 June 2020. 4. 4.World Health Organisation. Coronavirus disease (COVID-19) Situation Report – 167 (5 July 2020) 5 July 2020. 5. 5.Aldibasi OS, Alharbi NK, Alkelya M, Zowawi H, Alghnam S. The Association of Country-Level Factors with Outcomes of COVID-19: Analysis of the pandemic after one million cases. Research Square. 2020. 6. 6.Burden SJ, Rademaker J, Weedon BD, Whaymand L, Dawes H, Jones A. Associations of Global Country Profiles and Modifiable Risk Factors with COVID-19 Cases and Deaths. medRxiv. 2020:doi:2020.06.17.20133454. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=2020.06.17.20133454&link_type=DOI) 7. 7.Sorci G, Faivre B, Morand S. Why Does COVID-19 Case Fatality Rate Vary Among Countries? The Lancet. 2020. 8. 8.Garcia de Alcaniz JG, Romero-Lopez J, Martinez RP, Lopez-Rodas V, Costas E. What variables can better predict the number of infections and deaths worldwide by SARS-CoV-2? Variation through time. medRxiv. 2020:doi:2020.06.04.20122176. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=2020.06.04.20122176&link_type=DOI) 9. 9.Woody S, Garcia Tec M, Dahan M, Gaither K, Lachmann M, Fox S, et al. Projections for first-wave COVID-19 deaths across the US using social-distancing measures derived from mobile phones. medRxiv. 2020:2020.04.16.20068163. 10. 10.Courtemanche C, Garuccio J, Le A, Pinkston J, Yelowitz A. Strong Social Distancing Measures In The United States Reduced The COVID-19 Growth Rate. Health Affairs. 0:10.1377/hlthaff.2020.00608. 11. 11.Bray I, Gibson A, White J. Coronavirus disease 2019 mortality: a multivariate ecological analysis in relation to ethnicity, population density, obesity, deprivation and pollution. Public Health. 2020; 185:261–3. 12. 12.Worldometer. Covid-19 coronavirus pandemic. [cited 2020, 25 June 2020]; Available from: [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). 13. 13.Wikipedia. National responses to the COVID-19 pandemic. [cited 2020, 25 June 2020]; Available from: [https://en.wikipedia.org/wiki/National\_responses\_to\_the\_COVID-19\_pandemic](https://en.wikipedia.org/wiki/National\_responses\_to_the_COVID-19_pandemic). 14. 14.Google. COVID-19 Community Mobility Reports. [cited 2020, 25 June 2020]; Available from: [https://www.google.com/covid19/mobility/](https://www.google.com/covid19/mobility/). 15. 15.Wicklin R. What is loess regression? SAS Blogs: SAS; 2016 [cited 2020 29 June 2020]; Available from: [https://blogs.sas.com/content/iml/2016/10/17/what-is-loess-regression.html](https://blogs.sas.com/content/iml/2016/10/17/what-is-loess-regression.html). 16. 16.Eurostat. Eurostat Database. [cited 2020 25 June 2020]; Available from: [https://ec.europa.eu/eurostat/data/database](https://ec.europa.eu/eurostat/data/database). 17. 17.CEIC data. Visitor Arrivals. [cited 2020, 25 June 2020]; Available from: [https://www.ceicdata.com/en/indicator/russia/visitor-arrivals](https://www.ceicdata.com/en/indicator/russia/visitor-arrivals), [https://www.ceicdata.com/en/indicator/turkey/visitor-arrival](https://www.ceicdata.com/en/indicator/turkey/visitor-arrival), [https://www.ceicdata.com/en/indicator/ukraine/visitor-arrivals](https://www.ceicdata.com/en/indicator/ukraine/visitor-arrivals). 18. 18.Analytical Credit Rating Agency ACRA. Air traffic increases, finances in crisis [cited 2020, 25 June 2020]; Available from: [https://www.acra-ratings.cn/research/1282](https://www.acra-ratings.cn/research/1282). 19. 19.Hurriyet Daily News. Air passenger traffic in Turkish airports rose 8.8 percent in 2018. [cited 2020, 25 June 2020]; Available from: [https://www.hurriyetdailynews.com/air-passenger-traffic-in-turkish-airports-rose-8-8-percent-in-2018-140322](https://www.hurriyetdailynews.com/air-passenger-traffic-in-turkish-airports-rose-8-8-percent-in-2018-140322). 20. 20.Kyiv Post. Record-breaking 20.5 million passengers at Ukrainian airports in 2018. [cited 2020, 25 June 2020]; Available from: [https://www.kyivpost.com/ukraine-politics/record-breaking-20-5-million-passengers-passed-through-ukrainian-airports-in-2018.html](https://www.kyivpost.com/ukraine-politics/record-breaking-20-5-million-passengers-passed-through-ukrainian-airports-in-2018.html). 21. 21.Mukaka MM. Statistics corner: A guide to appropriate use of correlation coefficient in medical research. Malawi Med J. 2012; 24:69–71. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23638278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F05%2F2020.11.04.20225656.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309277200006&link_type=ISI) 22. 22.Calcagno V, de Mazancourt C. glmulti: An R Package for Easy Automated Model Selection with (Generalized) Linear Models. Journal Of Statistical Software. 2010; 34. 23. 23.Guterres A. COVID-19 in an Urban World. United Nations,; 2020 [cited 2020 13th August 2020]; Available from: [https://www.un.org/en/coronavirus/covid-19-urban-world](https://www.un.org/en/coronavirus/covid-19-urban-world). 24. 24.Barr J, Tassier T. Are Crowded Cities the Reason for the COVID-19 Pandemic? Placing too much blame on urban density is a mistake. Scientific American; 2020 [cited 2020 13th August 2020]; Available from: [https://blogs.scientificamerican.com/observations/are-crowded-cities-the-reason-for-the-covid-19-pandemic/](https://blogs.scientificamerican.com/observations/are-crowded-cities-the-reason-for-the-covid-19-pandemic/). 25. 25.Chaudhry R, Dranitsaris G, Mubashir T, Bartoszko J, Riazi S. A country level analysis measuring the impact of government actions, country preparedness and socioeconomic factors on COVID-19 mortality and related health outcomes. EClinicalMedicine. 2020; 25. 26. 26.Oztig LI, Askin OE. Human mobility and coronavirus disease 2019 (COVID-19): a negative binomial regression analysis. Public Health. 2020; 185:364–7.