Next Article in Journal
Promoting or Inhibiting? The Impact of Environmental Regulation on Corporate Financial Performance—An Empirical Analysis Based on China
Next Article in Special Issue
Impacts of COVID-19 on Trade and Economic Aspects of Food Security: Evidence from 45 Developing Countries
Previous Article in Journal
Perceived Versus Real Swimming Skills of Adolescents under Standard and Challenging Conditions: Exploring Water Competencies as an Approach to Drowning Prevention
Previous Article in Special Issue
The Exponential Phase of the Covid-19 Pandemic in Central Italy: An Integrated Care Pathway
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Spread of the Covid-19 Pandemic in Time and Space

Louvain Institute of Data Analysis and Modelling in Economics and Statistics, Institute of Statistics, Biostatistics and Actuarial Sciences, Université Catholique de Louvain, Voie du Roman Pays, 20, 1348 Louvain-la-Neuve, Belgium
Int. J. Environ. Res. Public Health 2020, 17(11), 3827; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17113827
Submission received: 11 May 2020 / Revised: 23 May 2020 / Accepted: 25 May 2020 / Published: 28 May 2020

Abstract

:
As the COVID-19 pandemic has had a profound impact on public health and global economies in 2020; it is crucial to understand how it developed and spread in time and space. This paper contributes to the growing literature by considering the dynamics of country-wise growth rates of infection numbers. Low-order serial correlation of growth rates is predominantly negative with cycles of two to four days for most countries. The results of fitted spatial autoregressive models suggest that there is high degree of spillover between countries. Forecast variances of many countries, in particular those with a high absolute number of infections, can to a large extent be explained by structural innovations of other countries. A better understanding of the serial and spatial dynamics of the spread of the pandemic may contribute to an improved containment and risk management.

1. Introduction

On December 31, 2019, the first case of a pneumonia caused by a new type of coronavirus was reported to the World Health Organization (WHO) office in Wuhan, China. Since then, the new disease, named COVID-19 on February 11, 2020, by the WHO, has developed and spread quickly throughout the world, touching essentially all countries by the end of April 2020, despite numerous containment strategies, social distancing, and the closing of borders.
The lockdown in many countries led to economic recessions associated with a slump of economic activity, rise of unemployment, and high uncertainty on financial markets. As shown, e.g., by Li et al. (2020), Wang et al. (2020), and Zhang and Ma (2020) [1,2,3], the pandemic had enormous psychological consequences for the population. Risk management has become essential during times of crises, and McAleer (2020) emphasizes the importance of prevention relative to the cure [4]. The outbreak and evolution of the COVID-19 pandemic has been analyzed from different perspectives. For example, Li and Linton (2020) fit country-wise quadratic regressions to estimate the peak periods [5]. Yue et al. (2020) analyze the impact of the pandemic on China’s economy [6], and Wang et al. (2020) consider the risk management of COVID-19 by universities in China [7]. Yue et al. (2020) propose solutions and recommendations related to early warning, identification, and monitoring of risks [8]. Liu et al. (2020) analyze the Chinese experience and its implications for other countries [9]. Chang et al. (2020) present a charter for a sustainable travel, tourism, and hospitality industry for the time after COVID-19 [10]. Chang and McAleer (2020) critically evaluate the Global Health Security Index (2019) [11], which provides data before the discovery of COVID-19 and makes it possible to evaluate how countries might have been prepared for a pandemic and acted accordingly. Further research analyzing the containment strategies of individual countries and global analysis of the COVID-19 situation includes Wang and Wang (2020), Zhao and Chen (2020), Han et al. (2020), and Di Gennaro et al. (2020) [12,13,14,15].
Despite the growing literature, there seems to be a gap in the understanding how the pandemic evolved over time, and how it spread to other countries. Rather than merely providing descriptive statistics, we would like to fit a joint model for the temporal and spatial dynamics that allows to measure the degrees of contagion between countries and yield predictions for the numbers of infected people. In this paper, we use spatial autoregressions to model daily growth rates of the number of incidences per country. Day-to-day relative changes in the number of incidences are informative about short-term effects of containment strategies and spillovers between countries. We refrain from modeling long-term trends in the data, which are notoriously difficult to capture, but rather focus on short-term dynamics, assuming stationarity for the investigated sample horizon.
It turns out that, perhaps surprisingly, there is a strong informational content in incidence growth rates, in the sense that they reveal negative autocorrelations and cyclical behavior. In a preliminary univariate analysis, we find strong negative low-order (i.e., mostly first order, but many even higher) autocorrelations of growth rates, as well as stochastic cycles of 2 to 4 days, meaning that fitted autoregressive models of order two have complex roots for the vast majority of countries. We embed these findings into a multivariate context, allowing for spillover (i.e., Granger causality in econometric terminology) between countries, as well as for spatial autocorrelation. The modeling framework is that of a spatial vector autoregression.
Spatial vector autoregressions have been used, e.g., by Wei (2015) to forecast influenza incidence rates of US states [16]. In a panel data context, similar models have been proposed by Elhorst (2003), Beenstock and Felsenstein (2007) and Lee and Yu (2012) [17,18,19]. Spatial vector autoregressions impose a structure on the covariance matrix of the error term that is motivated, e.g., by geographical factors such as distances. This structure enables a structural analysis of the vector autoregressive model (VAR) model that is not depending on a particular decomposition of the covariance matrix, a common problem in classical VAR models. In our case, we find that the best specification of the spatial weights matrix is based on inverse distances, modulated by population sizes.
We follow the approach of Diebold and Yilmaz (2014) to consider forecast error variance decompositions [20], interpreted as directed weighted networks. In our case, these decompositions depend on both the serial and spatial autocorrelation matrices. We find strong spillovers between countries in the sense that forecast error variances of many countries can, to a large extent, be explained by structural innovations of other countries. This seems to be particularly strong for countries with high absolute number of infections such as the USA, Italy, China, Germany, and Great Britain.
Our results indicate that geographical contagion of the COVID-19 pandemic has been and continues to be an important factor. The estimated spatial correlation is positive and highly significant, meaning that contagion is particularly strong between close and highly populated countries. These findings may have consequences for local and global containment and mitigation strategies.

2. Methodology and Results

We consider the daily number of newly infected people per country across time, as published by the European Centre for Disease Prevention and Control (ECSD), an agency of the European Union. Data are available for download at https://opendata.ecdc.europa.eu/covid19/casedistribution. In the early period of the epidemic, i.e., January and February of 2020, only few countries had significant numbers of infections, so we discard this period from the sample and consider the period of March and April, i.e., a sample size of 61 days. Population data are provided by the World Bank and are also available in the data sets provided by ECSD.
There are many alternatives to model the time trends and dynamic properties of the data, including quadratic trends as in Li and Linton (2020) [5]. We follow Yue et al. (2020) in considering autoregressions for the growth rates [6]. Let Y i t be the number of newly infected people in country i , i = 1 , , N , at time t, t = 1 , , T . We model growth rates y i t : = l o g Y i t + 1 l o g Y i , t 1 + 1 via spatial autoregressive processes of order p, i.e.,
y i t = μ i + j = 1 p α i j y i , t j μ i + ε i t
ε i t = ρ j = 1 N w i j ε j t + u i t , u i t N 0 σ 2
where w i j are spatial weights with w i i = 0 , i = 1…, N, and μ i , α i j are coefficients. The weights could, for example, be inversely related to the distance d i j between countries i and j. Beenstock and Felsenstein (2007) suggest to modulate inverse distances with the relative population size of both countries [18], that is
w i j = 1 d i j Z j Z i + Z j
where Z i is the population of country i. Hence, bigger countries receive more weight than smaller countries. The growth rates of the nine countries with highest number of cases until end of April, 2020, are depicted in Figure 1.
The error term can be written in vector form, stacking it into an (N × 1) vector ε t , as
ε t = ρ W ε t + u t
where W is the spatial weight matrix W with w i i = 0 , i = 1…, N. An advantage of the spatial error model is that, once ρ is estimated, the structural shocks u t can readily be obtained by the transformation I N ρ W ε t . This does not require decomposing a variance–covariance matrix as in classical VAR models, for which different methods are often a matter of debate in the empirical literature, see e.g., the discussion in Lütkepohl (2005) [21].
The dynamics of y t can be generalized to a vector autoregressive model (VAR) of the form
y t = μ + j = 1 p A j y t j μ + ε t
where ε t is a spatial error term as in Equation (4), and A j are N × N parameter matrices. For large N, restrictions are necessary for A j as otherwise there would be too many parameters to estimate for reasonable sample sizes. Indeed, in our case, the number of countries is larger than the number of days in the sample, so that methods are required that penalize model complexity.
To gain insight into the dynamics of infection growth rates, we will first look at univariate autoregressive models because it reveals some interesting dynamic features of the data. We will then consider spatial vector autoregressions to model the spatial and temporal relationships between countries. Of the more than 200 countries that have reported Covid-19 infections, we have selected the 100 countries with the highest number of infections as of end of April 2020. Thus, our cross-section dimension is given by N = 100. The reason to exclude countries with smaller number of cases is the highly erratic behavior of corresponding growth rates that could potentially bias the estimation of spatial correlation patterns.

2.1. Univariate Analysis

We fit autoregressive models of order p to the growth rate of infections for each country, where the lag order p is selected by the Akaike information criterion (AIC). The maximum lag chosen is 10, but for no country the selected p is larger than 2. Results are reported in Table A1 in the Appendix.
The empirical results suggest that:
  • First order serial correlations of growth rates are negative.
  • For the countries with AR(2) dynamics, there are stochastic cycles in growth rates with an average length of about three days.
  • The spatial autocorrelation coefficient ρ is strongly significant.
It is quite remarkable that except for one country (Andorra), all first order autocorrelation coefficients are negative with an average of −0.46. The second order autocorrelation coefficients are closer to zero, but still quite significant for many.
It turns out that among the 100 countries with the highest number of reported cases, 72 have growth rates with AR(2) dynamics, and of these all but one (Denmark) have roots of the characteristic equations that are complex, indicating the presence of stochastic cycles. The average length of these cycles is calculated as
k i = 2 π a r c c o s α i 1 / 2 α i 2
All average cycle durations happen to be between 2 and 4 days, with an average of 2.8 days across countries. We do not have an explanation for the presence of such cycles, and whether they are genuine to the flux of new infections within a country, or artificially generated by the reporting practice. However, given that almost all autocorrelations of order larger than two tend to be small and often negligible, we can exclude seasonal or weekend effects as potential explanation. We next move on to the multivariate analysis in order to understand potential spillovers between countries.

2.2. Multivariate Analysis

For the specification of the VAR model in Equation (5), we first choose the lag order of p = 2, which corresponds to the maximum of the selected lag orders found for the univariate models. We then estimate the model by minimizing equation by equation the following criterion
t = 1 T ε i t 2 + λ i j = 1 N k = 1 2 A i j k
with respect to the coefficients A i j k , where λ i is chosen by ten-fold cross-validation. This criterion is known in statistics as the least absolute shrinkage and selection operator (LASSO), see e.g., Tibshirani (1996) [22]. Figure 2 and Figure 3 show the parameter estimates for the 50 × 50 subset of countries with most infections. The number of non-zero coefficients in A 1 by country ranges from 1 to 32 (out of 100), while the total number is 1741 (out of 10,000). Similarly, the number of non-zero coefficients in A 2 ranges from 6 to 35, with a total number of 1614 non-zero coefficients.
In a second step, the coefficient ρ of the spatial error component model in Equation (4) is estimated by Gaussian maximum likelihood, see e.g., Bivand et al. (2013) [23]. The model obviously depends on the specification of the spatial weights matrix W, for which we used inverse distances modulated by the relative population sizes, as in Equation (3). Distances are measured as the geographic distance in kilometers on the WGS ellipsoid between the centroids of two countries. This specification led to the highest likelihood among a set of alternative specifications including unweighted inverse distances, k-nearest neighbor and contiguity matrices. We obtain an estimate of 0.2805 with associated standard error of 0.0518, which suggests that spatial correlation is highly significant and important to be included in the subsequent analysis. After estimation, the residuals u t have been tested for remaining spatial autocorrelation using the Baltagi et al. (2007) spatial autocorrelation test for panel models [24], which is an extension of the classical Moran test. The p-value of the test is 0.157, which suggests that spatial autocorrelation has been sufficiently captured by the model.
In the following, we study the implications of the estimated parameters for the decompositions of the variances of forecast errors, in the spirit of Diebold and Yilmaz (2014) [20]. This allows quantifying the network relationships between countries. Starting from the reduced form VAR (2) model in Equation (5) with spatial error terms ε t in Equation (4), we may obtain the infinite order vector moving average (VMA) representation,
y t = j = 0 Φ j ε t j
with coefficient matrices Φ j that, for our VAR (2) model, can be obtained recursively as Φ 0 = I N , Φ 1 = A 1 , Φ 2 = Φ 1 A 1 + A 2 , , Φ j = Φ j 1 A 1 + Φ j 2 A 2 , , see e.g., Lütkepohl (2005) [21].
From the spatial error model in Equation (4) we then obtain orthogonalized shocks as u t = I N ρ W ε t and rewrite the model as
y t = j = 0 Θ j u t j
where Θ j : = Φ j I N ρ W 1 . Thus, having orthogonalized the error terms, we can perform structural analysis, and this is independent of a particular decomposition of a variance-covariance matrix, which is a common problem in classical VAR models. The proportion of the h-step forecast error variance of country i accounted for by innovations in country j is then given by
ω i j , h = k = 0 h 1 θ i j , k 2 / j = 1 N k = 0 h 1 θ i j , k 2
where θ i j , k is the ij-th element of the parameter matrix Θ k . It is important to emphasize that ω i j , h depends on both the dynamics of the VAR system via the parameter matrices Φ j , and the spatial lag matrix, I N ρ W . Hence, forecast variance decompositions are determined by spatial and temporal dynamics jointly, which corroborates the importance of a joint spatio-temporal modeling approach.
The proportions ω i j , h are reported graphically in Figure 4 for the 50 countries having the highest number of infections. We chose a forecast horizon of h = 10, which for our parameter estimates is almost identical to long-run forecasts ( h = ) as ω i j , h converges quickly to a constant as h increases. Not surprisingly, most countries’ forecast errors are explained to a good extent by their own history, i.e., the diagonal elements of ω tend to be larger than the off-diagonal elements. The average of the diagonal elements is 0.47 with a maximum of 0.96 and a minimum of 0.1. This also implies that there are many countries for which the sum of the contributions of other countries is larger than that of their own history.
This analysis can be refined following Diebold and Yilmaz (2014) by viewing variance decompositions as weighted directed networks [20]. They define two directional connectedness measures, “to” and “from”, for each country. The “to” measure is defined as
ω j : = i = 1 , i j N ω i j
where we suppress the h-index for simplicity. The ω j measures the sum of the contributions of country j to all other countries’ forecast errors, and it can be viewed as a “to”-degree of a node (i.e., a country) of the network. The support of the univariate distribution of ω j is 0 , N . The countries’ “to”-degrees are visualized in Figure 5, and the corresponding distribution in Figure 6. The distribution of “to”-degrees has a bi-modal structure, where the majority of “to”-degrees is close to zero, and a small group of countries having “to”-degrees of 2 or larger. The countries with highest “to”-degrees are Ghana, United Arab Emirates, Philippines, Ecuador, Oman, Belarus, Singapore, Iraq, Spain, and Egypt. Note that each continent has at least one country with a high “to”-degree, i.e., a country whose innovations in infection rates help to predict infection rates of other countries.
Likewise, the “from” directional connectedness is defined as
ω i : = j = 1 , j i N ω i j
And measures the sum of the contributions of all other countries j to explain the forecast error variance of country i. The support of the univariate distribution of ω i is 0 , 1 . These measures are visualized in Figure 7, and the corresponding distribution in Figure 8. Note that, by construction, the means of the “to” and “from” distributions are the same, 0.5222, but that the “from” distribution is more concentrated around values close to one, meaning that there are many countries whose forecast variances can to a large proportion be explained by the innovations of other countries. The countries with highest “from”-degrees are Canada, Italy, USA, China, Belgium, Bulgaria, Mexico, Germany, Great Britain, and Denmark. Note that these countries tend to have high absolute number of infection cases. Note also that the proportion of developing and emerging countries among those with high “to”-degrees is larger than among those with high “from”-degrees.
Finally, the total connectedness measure for the network of countries is given by ω : = i = 1 N ω i = j = 1 N ω j , which corresponds to the mean degree of the network. We obtain a mean degree of 52.22.

3. Conclusions

Our findings suggest that daily growth rates of new infections per country are negatively autocorrelated and show stochastic cycles of about two to four days. In a spatial vector autoregression, we find high spillovers between countries, motivated by highly significant serial and spatial autocorrelations. Forecast error variances are to a large extent explained by other countries’ innovations. This holds especially for countries with high absolute numbers of infections.
High contagion degrees between countries requires risk and crisis management that takes the evolution of the pandemic in other countries into account. We have seen the pandemic starting to develop in China, then spreading to Europe and the USA, and finally, developing strongly in countries such as Russia and Brazil. This paper contributes to a better understanding of how the crisis has evolved both over time and across countries, which may be fruitful for an improved risk management of future crises.

Funding

This research received no external funding.

Conflicts of Interest

The author declares has no conflict of interest.

Appendix

Table A1. Estimation results for the univariate AR(p) models.
Table A1. Estimation results for the univariate AR(p) models.
CountryMuAlpha1Alpha2KAcf1Acf2
1USA0.1190NANANA−0.58120.1025
2ESP0.0551−0.6357−0.30702.879772−0.4078−0.0964
3ITA0.0197−0.2187NANA−0.3080−0.1204
4GBR0.0994−0.2747−0.31703.457548−0.4656−0.1669
5DEU−0.0483NANANA−0.75260.4425
6FRA0.0591−0.5001NANA−0.51800.0404
7TUR0.12900.2093NANA−0.46030.1770
8RUS0.1454−0.6229−0.30602.897091−0.2381−0.1603
9IRN0.0153−0.7193−0.34992.824746−0.3826−0.2127
10BRA0.1456−0.5581NANA−0.74320.4734
11CAN0.0951−0.3906NANA−0.4648−0.0175
12BEL0.0951NANANA−0.3717−0.1870
13NLD0.0705NANANA−0.13790.0532
14PER0.1315−0.7798−0.37402.777552−0.4865−0.0313
15IND0.1246−0.7751−0.19022.357376−0.4758−0.1840
16CHE0.0532−0.6605−0.18182.557440−0.55810.1938
17PRT0.1035−0.4378NANA−0.2475−0.1414
18ECU0.0618−0.3215−0.35903.410295−0.68460.3053
19SAU0.1182−0.7888−0.50282.908032−0.59230.0215
20SWE0.0980−0.2257NANA−0.39060.0849
21IRL0.0965−0.3648NANA−0.42580.0670
22MEX0.1077−0.3284−0.19603.220894−0.4257−0.0989
23PAK0.1142−0.6143−0.31252.919015−0.53960.1344
24SGP0.0764−0.6688−0.33312.870741−0.44100.0195
25CHL0.1154−0.4258NANA−0.56790.1749
26ISR0.0548−0.8392−0.26852.498557−0.41770.1557
27AUT0.0481−0.2929−0.23003.340180−0.52340.3204
28JPN0.0409−0.3556−0.30843.312745−0.68710.3177
29BLR0.1105−0.5120−0.37743.140573−0.2509−0.0071
30QAT0.0925−0.8442−0.33082.623628−0.22740.2722
31POL0.0897−0.4912NANA−0.71590.4200
32ARE0.1035−0.9655−0.45362.651255−0.49710.0442
33ROU0.0913−0.4175NANA−0.2366−0.2829
34IDN0.0779−0.7072−0.33862.825275−0.4945−0.1337
35UKR0.0000−0.5586−0.37373.071974−0.37390.0234
36DNK0.0709−0.17450.2325NaN−0.5398−0.0478
37SRB0.0927−0.5346−0.34793.078214−0.50010.1388
38PHL0.0922−0.8493−0.33432.622693−0.0336−0.0254
39NOR0.0357−0.7009−0.42922.942658−0.65480.2291
40CZE0.0534−0.2509NANA−0.57570.3656
41BGD0.1039−0.6032−0.29782.913854−0.65750.1833
42KOR−0.0693−0.7663−0.32362.720150−0.53320.0641
43DOM0.0833−0.6311−0.26952.825041−0.55450.0766
44AUS0.0237−0.4603NANA−0.4901−0.0438
45PAN0.0827−0.3611NANA−0.65130.3146
46COL0.0933−0.6590−0.22122.677165−0.52830.0350
47MYS0.0402−0.8145−0.45812.834799−0.53290.0334
48ZAF0.0934−0.7034−0.22442.609739−0.65580.2550
49EGY0.0804−0.7420−0.50072.959918−0.36480.2310
50FIN0.0510−0.7983−0.47872.874698−0.66160.2868
Notes: Countries are indexed by their ISO3 code. The number mu is the average growth rate. AR orders were selected by the AIC criterion. If p = 0, then both α1 and α2 are marked as NA, if p = 1, then α2 is marked as NA. The number k gives the average length of the stochastic cycles, in days, for those AR(2) models that have complex roots. For all other models, k is marked as NA; acf1 and acf2 are autocorrelation coefficients of order 1 and 2 for the growth rates.

References

  1. Li, S.; Wang, Y.; Xue, J.; Zhao, N.; Zhu, T. The Impact of COVID-19 Epidemic Declaration on Psychological Consequences: A Study on Active Weibo Users. Int. J. Environ. Res. Public Health 2020, 17, 2032. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wang, C.; Pan, R.; Wan, X.; Tan, Y.; Xu, L.; Ho, C.S.; Ho, R.C. Immediate Psychological Responses and Associated Factors during the Initial Stage of the 2019 Coronavirus Disease (COVID-19) Epidemic among the General Population in China. Int. J. Environ. Res. Public Health 2020, 17, 1729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Zhang, Y.; Ma, Z.F. Impact of the COVID-19 Pandemic on Mental Health and Quality of Life among Local Residents in Liaoning Province, China: A Cross-Sectional Study. Int. J. Environ. Res. Public Health 2020, 17, 2381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. McAleer, M. Prevention Is Better Than the Cure: Risk Management of COVID-19. J. Risk Financ. Manag. 2020, 13, 46. [Google Scholar] [CrossRef] [Green Version]
  5. Li, S.; Linton, O. When Will the Covid-19 Pandemic Peak? Cambridge INET Working Paper Series 2020/11. Available online: https://www.inet.econ.cam.ac.uk/research-papers/wp-abstracts?wp=2011 (accessed on 27 May 2020).
  6. Yue, X.-G.; Shao, X.-F.; Li, R.Y.M.; Crabbe, M.J.C.; Mi, L.; Hu, S.; Baker, J.S.; Liu, L.; Dong, K. Risk Prediction and Assessment: Duration, Infections, and Death Toll of the COVID-19 and Its Impact on China’s Economy. J. Risk Financ. Manag. 2020, 13, 66. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, C.; Cheng, Z.; Yue, X.-G.; McAleer, M. Risk Management of COVID-19 by Universities in China. J. Risk Financ. Manag. 2020, 13, 36. [Google Scholar] [CrossRef] [Green Version]
  8. Yue, X.G.; Shao, X.F.; Li, R.Y.; Crabbe, M.J.; Mi, L.; Hu, S.; Baker, J.S.; Liang, G. Risk Management Analysis for Novel Coronavirus in Wuhan, China. J. Risk Financ. Manag. 2020, 13, 22. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, W.; Yue, X.-G.; Tchounwou, P.B. Response to the COVID-19 Epidemic: The Chinese Experience and Implications for Other Countries. Int. J. Environ. Res. Public Health 2020, 17, 2304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Chang, C.-L.; McAleer, M.; Ramos, V. A Charter for Sustainable Tourism after COVID-19. Sustainability 2020, 12, 3671. [Google Scholar] [CrossRef]
  11. Chang, C.-L.; McAleer, M. Alternative Global Health Security Indexes for Risk Analysis of COVID-19. Int. J. Environ. Res. Public Health 2020, 17, 3161. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, J.; Wang, Z. Strengths, Weaknesses, Opportunities and Threats (SWOT) Analysis of China’s Prevention and Control Strategy for the COVID-19 Epidemic. Int. J. Environ. Res. Public Health 2020, 17, 2235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhao, S.; Chen, H. Modeling the epidemic dynamics and control of COVID-19 outbreak in China. Quant. Biol. 2020, 8, 11–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Han, X.; Wang, J.; Zhang, M.; Wang, X. Using Social Media to Mine and Analyze Public Opinion Related to COVID-19 in China. Int. J. Environ. Res. Public Health 2020, 17, 2788. [Google Scholar] [CrossRef] [PubMed]
  15. Di Gennaro, F.; Pizzol, D.; Marotta, C.; Antunes, M.; Racalbuto, V.; Veronese, N.; Smith, L. Coronavirus Diseases (COVID-19) Current Status and Future Perspectives: A Narrative Review. Int. J. Environ. Res. Public Health 2020, 17, 2690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wei, X. Spatial Panel VAR and Application to Forecast Influenza Incidence Rates of US States, SSRN Discussion Paper 2015. Available online: https://ssrn.com/abstract=2646870 (accessed on 27 May 2020).
  17. Elhorst, J.P. Specification and estimation of spatial panel data models. Int. Reg. Sci. Review. 2003, 26, 244–268. [Google Scholar] [CrossRef]
  18. Beenstock, M.; Felsenstein, D. Spatial Vector Autoregressions. Spat. Econ. Anal. 2007, 2(2), 167–196. [Google Scholar] [CrossRef] [Green Version]
  19. Lee, L.; Yu, J. Spatial panels: Random components versus fixed effects. Int. Econ. Rev. 2012, 53, 1369–1412. [Google Scholar] [CrossRef]
  20. Diebold, F.X.; Yilmaz, K. On the network topology of variance decompositions: Measuring the connectedness of financial firms. J. Econom. 2014, 182, 119–134. [Google Scholar] [CrossRef] [Green Version]
  21. Lütkepohl, H. New Introduction to Multiple Time Series Analysis; Springer Verlag: New York, NY, USA, 2005. [Google Scholar]
  22. Tibshirani, R. Regression Shrinkage and Selection via the lasso. J. R. Stat. Soc. Ser. B. 1996, 58, 267–288. [Google Scholar] [CrossRef]
  23. Bivand, R.S.; Pebesma, E.; Gomez-Rubio, V. Applied Spatial Data Analysis with R; Springer Verlag: New York, NY, USA, 2013. [Google Scholar]
  24. Baltagi, B.H.; Song, S.H.; Jung, B.; Koh, W. Testing panel data regression models with spatial and serial error correlation. J. Econom. 2007, 140, 5–51. [Google Scholar] [CrossRef]
Figure 1. Growth rates of the nine countries with highest number of cases until end of April 2020: USA, Spain, Italy, Great Britain, Russia, Germany, Brazil, France, Turkey.
Figure 1. Growth rates of the nine countries with highest number of cases until end of April 2020: USA, Spain, Italy, Great Britain, Russia, Germany, Brazil, France, Turkey.
Ijerph 17 03827 g001
Figure 2. Estimated parameter matrix A1.
Figure 2. Estimated parameter matrix A1.
Ijerph 17 03827 g002
Figure 3. Estimated parameter matrix A2.
Figure 3. Estimated parameter matrix A2.
Ijerph 17 03827 g003
Figure 4. Decompositions of forecast variances.
Figure 4. Decompositions of forecast variances.
Ijerph 17 03827 g004
Figure 5. Directional “TO” connectedness measure, or “TO”-degrees.
Figure 5. Directional “TO” connectedness measure, or “TO”-degrees.
Ijerph 17 03827 g005
Figure 6. Kernel estimate of the density of “TO”-degrees.
Figure 6. Kernel estimate of the density of “TO”-degrees.
Ijerph 17 03827 g006
Figure 7. Directional “FROM” connectedness measure, or “FROM”-degrees.
Figure 7. Directional “FROM” connectedness measure, or “FROM”-degrees.
Ijerph 17 03827 g007
Figure 8. Kernel estimate of the density of “FROM”-degrees.
Figure 8. Kernel estimate of the density of “FROM”-degrees.
Ijerph 17 03827 g008

Share and Cite

MDPI and ACS Style

Hafner, C.M. The Spread of the Covid-19 Pandemic in Time and Space. Int. J. Environ. Res. Public Health 2020, 17, 3827. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17113827

AMA Style

Hafner CM. The Spread of the Covid-19 Pandemic in Time and Space. International Journal of Environmental Research and Public Health. 2020; 17(11):3827. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17113827

Chicago/Turabian Style

Hafner, Christian M. 2020. "The Spread of the Covid-19 Pandemic in Time and Space" International Journal of Environmental Research and Public Health 17, no. 11: 3827. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17113827

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop