Next Article in Journal
Influence of Financial Variables on the Development of Rural Communes of Eastern Poland in 2009–2018
Previous Article in Journal
Management of Distribution Risks and Digital Transformation of Insurance Distribution—A Regulatory Gap in the IDD
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multivariate Analysis of Energy Commodities during the COVID-19 Pandemic: Evidence from a Mixed-Frequency Approach

1
Scuola Normale Superiore, 56126 Pisa, PI, Italy
2
MEMOTEF Department, Sapienza University of Rome, 00161 Rome, RM, Italy
3
Department of Statistical Sciences, Sapienza University of Rome, 00161 Rome, RM, Italy
*
Author to whom correspondence should be addressed.
Risks 2021, 9(8), 144; https://doi.org/10.3390/risks9080144
Submission received: 8 July 2021 / Revised: 29 July 2021 / Accepted: 3 August 2021 / Published: 11 August 2021

Abstract

:
This paper shows the effects of the COVID-19 pandemic on energy markets. We estimate daily volatilities and correlations among energy commodities relying on a mixed-frequency approach that exploits information from the number of weekly deaths related to COVID-19 in the United States. The mixed-frequency approach takes advantage of the MIxing-Data Sampling (MIDAS) methods. We compare our results to those obtained by employing two well-known models that do not account for the COVID-19 low-frequency variable, namely the Dynamic EquiCorrelation (DECO) and corrected Dynamic Conditional Correlation (cDCC). Moreover, we consider four possible specifications of the volatility: GARCH, GJR, GARCH-MIDAS, and Double-Asymmetric GARCH-MIDAS. The empirical results show that our approach is statistically superior to other models and represents a valuable methodology that can be used for risk managers, investors, and policy makers to assess the effects of the pandemic on spillovers effects in energy markets.

1. Introduction

During the last two decades, the interest of investors in energy commodities has increased, especially during financial and economic recessions. This process, referred to as the financialization of energy commodities, along with the deregulation of the over the counter markets, has created significant shifts in commodity returns volatility (Falkowski 2011; Masters and White 2008). Given the relevant role of energy commodities on the inflation and the industrial production (Cunado and de Gracia 2003; Hamilton 1983) and their employment as hedging and speculative assets ((Arouri et al. 2012; Kang et al. 2015); Sadorsky 1999, 2003), many contributions to the econometric and financial literature have proposed models aimed at studying the volatility of energy commodities, especially for risk management purposes (Aloui and Mabrouk 2010; Denton et al. 2003; Fan et al. 2008; Hung et al. 2008; Krehbiel and Adkins 2005).
Popular GARCH (Bollerslev 1986) models have been extensively used by practitioners and scholars to account for well-known stylized facts of energy commodities such as asymmetry, leptokurtosis, leverage effect, and volatility clustering. For instance, Lv and Shan (2013) compare the performance of different specifications of GARCH models to forecast natural gas volatility, whereas Laporta et al. (2018) investigate the performance of different GARCH specifications for modeling the Value-at-Risk of seven energy commodities. However, univariate models are not able to explain dependences in the dynamics of the volatility of energy commodities. For this reason, many studies employ multivariate GARCH models (surveyed by Bauwens et al. 2006; Silvennoinen and Teräsvirta 2009) to study the comovements of commodities returns, capturing time-varying volatilities and correlations. The Dynamic Conditional Correlation (DCC) class of models (Engle 2002) achieves such a goal, and it has been widely applied to energy markets. For instance, Sadorsky (2012) finds that the DCC model performs better than other multivariate GARCH models in capturing the interdependencies among the West Texas Intermediate (WTI) oil returns and two weighted returns’ indexes of clean energy and technology companies. Wang and Wu (2012) use the DCC model to forecast the volatility and correlations among the returns of crude oil, gasoline, heating oil, and jet fuel. Behmiri et al. (2019) investigate the relationships among natural gas, oil and other non-energy commodities, Mishra (2019) studies the correlation between oil and foreign exchange markets, whereas Basher and Sadorsky (2016|) focus on the connection between stock indices and oil returns.
The rapid spread of COVID-19 favored the transmission of shocks from the real economy to the financial markets, emphasizing the importance of developing accurate analytic models to assess and quantify such shocks. Recently, many contributions investigated the effects of the recent COVID-19 outbreak to understand the real effects of the pandemic on energy commodity returns. For instance, Ghorbel and Jeribi (2021) investigate the interconnections between the volatilities of crude oil, gas, energy index, and financial assets’ returns during the COVID-19 pandemic through a Markov Switching multivariate GARCH model; Narayan (2020) evaluates the effect of COVID-19 infections and oil price news on oil prices; Gil-Alana and Monge (2020) study the persistence of the outbreak shock on oil prices; Iyke (2020) and Foglia and Angelini (2020) identify the volatility spillovers between oil and clean energy companies; Gil-Alana and Monge (2020) and Bouazizi et al. (2021) study the effect of the pandemic on the volatility of oil returns; Corbet et al. (2020), Lin and Su (2021), and De Blasis and Petroni (2021) study the connected relationship between oil and other energy markets, such as renewables and coal, employing the DCC model. This strand of literature is particularly relevant since it studies both the behavior of energy commodities in the extraordinary and unexpected circumstances of the pandemic, which caused a strong rise in volatility levels, as well as in the unprecedented situation of negative prices. In April 2020, crude oil recorded a 300% drop in price, assuming a negative price for the first time in the history. The rise in volatility of these commodities is found to be mainly due to government restrictions (such as the lockdown and travel restrictions), which lowered the global demand of energy commodities, such as crude oil and Kerosene-Type Jet Fuel, and, in turn, significantly affected their prices.
The aim of our paper is to investigate the time-varying correlations of energy commodities during the whole pandemic period in order to understand the deep impact of the COVID-19. The period considered starts from the first wave in March 2020 until the midst of the vaccination campaign in June 2021, corresponding to the beginning of the economic upturn and the consequent recovery of energy commodities’ prices. We investigate daily returns of six of the main energy commodities: WTI crude oil, Brent crude oil, Heating oil #2, Propane, New York Harbor Conventional Gasoline Regular, and Kerosene-Type Jet Fuel, by means of a set of DCC-GARCH models, which is the corrected DCC (cDCC, Aielli 2013), the DCC-MIDAS (Colacito et al. 2011), and the Dynamic Equicorrelation (DECO, Engle and Kelly 2012). A distinctive feature of the DCC class of models is the possibility to split the estimation phase into two steps: one for the univariate volatilities and another for the correlations.
In order to assess the impact of the COVID-19 outbreak on the volatility of energy markets, we include in the analysis the number of weekly deaths in the United States (US) caused by the COVID-19, which is the variable most used by governments to impose restrictive measures. The use of such a variable sampled at a lower frequency than daily energy prices leads us to rely on the class of MIDAS models (Ghysels et al. 2007). In particular, we exploit the number of US weekly deaths in the univariate specifications of the GARCH-MIDAS (GM, Engle et al. 2013) and Double Asymmetric GARCH-MIDAS (DAGM, Amendola et al. 2019). The inclusion of the DAGM model allows us to evaluate the effects of positive and negative variations of the COVID-19-related deaths in the US on the daily volatility of commodity returns. To the best of our knowledge, the assessment of the potential effect of the COVID-19 pandemic on the interdependences of daily commodity returns through low-frequency variables has not yet been explored in the literature.
We evaluate the performance of our approach comparing the DCC models employed using the GM and DAGM specifications with GARCH model and the GJR (Glosten et al. 1993) specification, which includes an asymmetric term related to negative daily returns. In order to choose the best model able to detect the phenomena under investigation, we perform a model selection procedure. In particular, we consider the Model Confidence Set of Hansen et al. (2011) to create a Set of Superior Models (SSM) relying on a particular loss function. Due to the correlation structure of the models, we consider three robust loss functions (Laurent et al. 2013) based on the distance between the predicted conditional covariance matrix and its proxy, built as the cross-products of the daily log-returns. The same loss functions have also been used in Amendola et al. (2020), among others.
Regardless of the models applied in the multivariate step, our empirical results show that for two out of three loss functions only models with the DAGM specification enter in the SSM. The time-varying conditional correlations estimated with the models belonging to the SSM remain constant through the whole pandemic period, ranging between 0.6 and 0.9, which indicates strong comovements among commodities. Considering the importance of energy commodities in stock markets and macroeconomic fundamentals, the presence of a strong correlation between energy markets is translated in high risks of spillovers effects from commodities to financial markets and to the real economy. The results obtained in this paper confirm that the proposed methodology represents a valuable approach, which is statistically superior to others to discover fragilities in the economic and financial systems that, regardless of the economic upturn, still need to be addressed at the present time.
The contribution of our paper to the literature is two-fold. First, our paper considers for the first time the impact of the COVID-19 on six of the main energy commodities. The choice of energy commodities depends on the fact that they are exploited in industrial applications and used for hedging financial risks as discussed in (Abid et al. 2020). As a matter of fact, the COVID-19 epidemic threatened the stability of the oil and gas industries, which were already weakened by the increasing use of renewables. Moreover, the mobility restrictions imposed during the first and second COVID-19 waves (Arenas et al. 2020; Nouvellet et al. 2021; Sadowski et al. 2021) shocked the aviation and transportation sectors, which represent 60% of the demand of oil and other commodities, including Propane and Kerosene-Type Jet Fuel used in our analysis (IEA 2020). Previous studies, instead, mainly focus on the effects of COVID-19 on oil and natural gas (De Blasis and Petroni 2021; Gil-Alana and Monge 2020; Lin and Su 2021; Narayan 2020). In particular, Nyga-ukaszewska and Aruga (2020) investigate the impact of the COVID-19 cases in the US and Japan on the crude oil and natural gas markets by means of Auto-Regressive Distributive Lag approach. Wang and Su (2021) examine the asymmetric causal relationship between COVID-19 and fossil energy prices and find that COVID-19 influences oil, natural gas, and coal prices. Second, we apply the mixed-frequency methodology, as in Candila (2021b), to the energy sector to evaluate the impact of COVID-19. Specifically, we investigate the effects of weekly deaths in the US due to COVID-19 on the daily energy commodities volatilities. We, therefore, contribute to the existing literature regarding the application of multivariate GARCH (Ku et al. 2007; Yousfi et al. 2021) to energy commodities (Chang et al. 2011; Chkili et al. 2014; Silvennoinen and Thorp 2013) where time series data are sampled only at the daily frequency. To the best of our knowledge, this is the first attempt to study the impact of COVID-19 to energy sector using a mixed-frequency approach.
The rest of the paper is organized as follows. Section 2 presents the models and different specifications used. Section 3 shows the empirical study and discusses the results obtained. Section 4 concludes.

2. Methodology

We denote with r i , t = [ r i , t , 1 , , r i , t , n ] the vector of daily log-returns observed at time i of period t, where t has a lower frequency than i. We consider T low-frequency periods, each composed of N T days. Denoting with M V N the multivariate normal distribution1, we model the returns as in Engle (2002):
r i , t = H i , t 1 / 2 ξ i , t
r i , t M V N ( 0 , H i , t )
ξ i , t i i d M V N ( 0 , I n )
where H i , t is the n × n conditional covariance matrix of returns, and H i , t 1 / 2 is obtained by the Cholesky decomposition of H t . In addition, ξ i , t is the ( n × 1 ) -vector of standardized residuals obtained through a univariate volatility model, and I n is an n × n identity matrix. Within the DCC-GARCH model, the conditional covariance matrix H i , t can be decomposed as follows:
H i , t = E i 1 , t ( r i , t r i , t ) = D i , t R i , t D i , t ,
where D i , t is the diagonal matrix containing the conditional standard deviations, and R i , t is the conditional correlation matrix of returns. The standardized residuals can be then expressed in terms of D i , t :
ξ i , t = D i , t 1 r i , t .
From Equation (4), R i , t is estimated as:
R i , t = E i 1 , t [ ξ i , t ξ i , t ] = D i , t 1 H i , t D i , t 1 .
In the DCC-GARCH model, the estimation procedure of the time-varying D i , t and R i , t is performed in two steps. First, a univariate volatility model is used to estimate the standardized residuals ξ i , t , and then R i , t is estimated according to a pre-specified correlation model. As a result, the model is governed by the parameter space Θ = Θ v , Θ c , where Θ v is the parameter space of the volatility model and Θ c that of the correlation model. The estimation of Θ is based on the global log-likelihood function L ( Θ ) :
L ( Θ ) = L ( Θ v ) + L ( Θ c ) = t = 1 T i = 1 N t n log ( 2 π ) + 2 log | D i , t | + r i , t D i , t 2 r i , t
t = 1 T i = 1 N t log | R i , t | + ξ i , t R i , t 1 ξ i , t + ξ i , t ξ i , t .
In this paper, we consider: three specifications of the correlation model (DECO, cDCC, and DCCMIDAS) and four univariate specifications of the volatility model (GARCH, GJR, GM, and DAGM) for a total of twelve parametric models. Regarding the univariate specifications of the volatility model, the GARCH model is described as follows:
σ i , t , j 2 = ω + α r i 1 , t , j 2 + β σ i 1 , t , j 2 , with j = 1 , , n .
The GJR specification is aimed at modeling the so-called leverage effect for which negative shocks have a higher impact on return volatility than positive shocks. Specifically:
σ i , t , j 2 = ω + α σ i 1 , t , j 2 + β ϵ i 1 , t , j + γ ϵ i 1 , t , j 2 𝟙 ϵ i 1 , t , j < 0 , with j = 1 , , n .
In the GM and DAGM models, the volatility is decomposed into a long-term and short-term component. The long-term component is denoted with τ t , j and varies at the same frequency as the low-frequency variable X t , j . The short-term component is g i , t , j and varies on a daily basis. The difference between the GM and the DAGM model relies on the specifications of the long-term component. Formally, r i , t , j and the short-term component for both the GM and DAGM models are:
r i , t , j = g i , t , j × τ t , j z i , t , j , with j = 1 , , n
g i , t , j = 1 α γ / 2 β + ( α + γ 𝟙 r i 1 , t , j < 0 ) r i 1 , t , j 2 τ t , j + β g i 1 , t , j
where z i , t , j N ( 0 , 1 ) , and γ is the parameter that accounts for the skewness. The long-run term for the GM model is:
τ t , j = e x p m + θ k = 1 K δ k ( ω ) X t k , j ,
whereas, in the DAGM model, the contributions of positive and negative values of X t are estimated separately, decomposing θ in Equation (13) into two parameters θ + and θ :
τ t , j = e x p m + θ + k = 1 K δ k ( ω ) + X t k , j 𝟙 X t k , j 0 + = θ k = 1 K δ k ( ω ) X t k , j 𝟙 X t k , j < 0 .
Regarding the correlation models, the functional forms of the cDCC, DCC-MIDAS, and DECO are shown in Table 1.

Model Selection

We select the best-performing models in terms of correlation estimation exploiting the Model Confidence Set (MCS) methodology (Hansen et al. 2003), which allows us to select an SSM from those considered in the analysis. Starting from a set of models M 0 , the MCS procedure allows us to select a subset of best models, the SSM, denoted with M 1 δ ٭ of dimension m ٭ m , with a given confidence level 1 δ . The MCS procedure consists of three steps: (i) collect the initial set of models M 0 used in the analysis; (ii) select the loss function and perform the test of equal predictive ability (EPA) for each model in M 0 with respect to that loss function; (iii) select the SSM according to the result of the EPA test. If the null hypothesis of the EPA test for a given subset of models is accepted, that subset is the SSM M 1 δ ٭ . If the null hypothesis is rejected, rule out the worst performing model and return to step (ii) until the set M 1 δ ٭ is defined.
More formally, recall that H i , t is the n × n conditional covariance matrix estimated through the multivariate GARCH in Equation (4) and denote with Σ i , t the true conditional variance matrix. Since Σ i , t is unobservable, the proxy Σ ^ i , t is considered and computed as the cross-products of r i , t . We apply the MCS procedure to select the SSM according to the loss function i , t l for the model l, with l = 1 , , M 0 , for the day i of the period t:
i , t l = ( Σ ^ i , t , H i , t l ) .
The evaluation of the models is done based on three robust loss functions (Laurent et al. 2013): Frobenius (FROB), Root Mean Squared Error (RMSE), and Eclidean (EUCL), defined as follows:
FROB : t r [ ( Σ ^ i , t H i , t l ) ( Σ ^ i , t H i , t l ) ] ;
RMSE : [ n o r m ( H i , t l Σ ^ i , t ) ] 1 / 2 ;
EUCL : [ v e c h ( Σ ^ i , t H i , t l ) v e c h ( Σ ^ i , t H i , t l ) ] ,
where t r ( · ) and n o r m ( · ) , respectively, denote the trace and norm of the matrix ( · ) , and v e c h ( · ) is the operator that stacks the lower triangular half of ( · ) into a single vector of length n ( n + 1 ) / 2 .
For ease of notation, we here suppress the double time index to only report the day i. In order to perform the EPA test, we define d l k , i , which is the loss differential of model l with respect to model k, and d l · , i , which is the average loss differential between model l and every other competing model contained in the set M, with M M 0 , at time i:
d l k , i = l , i k , i l , k M
and
d l · , i = 1 m 1 k M d l k , i , l , k M .
Model l has a better forecasting performance to model k when d l k , i < 0 . The null and alternative hypothesis of the EPA test are:
H 0 : E [ d l · ] = 0 , l M and H 1 : E [ d l · ] 0 , for some l M .
H 0 against H 1 is tested according to the Hansen et al. (2011) t-statistic:
t l , · = d ¯ l , · v a r ^ ( d ¯ l , · ) , l M ,
In particular, d ¯ l , · = 1 m 1 k M d ¯ l k , i represents the sample loss of model l related to the average losses across every other model k, d ¯ l k = 1 n d l k , i is the average loss between model l and k and v a r ^ ( d ¯ l , · ) is the bootstrap estimate of v a r ( d ¯ l , · ) . The EPA hypothesis is tested in Equation (21) with the statistic:
T max , M = max l M t l , · .
Large values of t l , · indicate that the distance between the estimates of the model and the actual realizations is greater for model l than any other model k M . This results in the rejection of the l-th model from the SSM according to an elimination rule. Following Hansen et al. (2011), an elimination rule coherent with the statistic (23) is:
E M = arg max l M d ¯ l , · v a r ^ ( d ¯ l , · ) = arg max l M t l · .
The aim of Equation (24) is to reject those models increasing t l · the most. At each iteration of the MCS procedure, if the null hypothesis of Equation (21) is rejected at the fixed significance level δ , the elimination rule in Equation (24) eliminates the model that delivers the largest standardized excess loss with respect to the average loss (computed across all other models in M). Then, the statistic from Equation (23) is recomputed for all l , k M M 0 . The procedure stops when the null hypothesis in Equation (21) is accepted, and the SSM M ^ 1 δ ٭ is determined.

3. Empirical Application

In this work, six energy commodities used in industrial applications are considered: WTI crude oil, Brent crude oil, Heating oil #2, Propane, New York Harbor Conventional Gasoline Regular, and Kerosene-Type Jet Fuel. These commodities are widely traded in financial markets and represent a valid sample to investigate the role of the COVID-19 pandemic on the fluctuations of energy commodities. They are used in portfolio analysis and are highly sought after by academics and practitioners to investigate the transmissions of commodity risks, as discussed in Naeem et al. (2021) and Chevallier and Ielpo (2013), which demonstrate that oil is a net transmitter of risk to other commodities. Furthermore, Karali and Ramirez (2014) discover volatility spillover effects between natural gas and crude oil and between natural gas and heating oil, making the transmission of commodity risks through the financial system more likely. The times series of daily data have been collected from the US Energy Information Administration (EIA) site and span from January 2020 (one month before the start of the COVID-19 pandemic) to June 2021, leading to a total of 276 daily observations. All the estimates have been computed in R, using the package dccmidas (Candila 2021a). The summary statistics of the considered commodities, and the related plots are reported in Table 2 and Figure 1.
Lin and Su (2021) points out that the connectedness of energy commodities experienced a peak during the pandemic outbreak. Along this line of research, our goal is to provide a quantitative measure of the impact of the COVID-19 weekly deaths on such commodities. We, therefore, consider the weekly number of deaths in the US related to COVID-19 infections from March 2020 to June 2021 (Figure 2) as a low-frequency variable. In particular, this variable is used in Equations (13) and (14) by taking the first differences.
The empirical application is developed in two parts: the first concerns the detection of structural breaks to investigate the effects of the COVID-19 outbreak on energy commodities. To do so, we apply the Bai–Perron test (Bai and Perron 1998) to detect both the presence and the number of structural breaks occurred at unknown dates during the period considered. The second part is related to the estimation and the comparison through the MCS procedure of twelve different DCC model specifications. We recall that in the analysis, we employed three specifications of the correlation model (DECO, cDCC, and DCC-GARCH-MIDAS) and four univariate specifications of the volatility model (GARCH, GJR, GM, and DAGM).
The presence of structural breaks coinciding with the beginning of the COVID-19 outbreak justifies the use of the low-frequency variable “COVID-19 deaths in the US” related to the effects of the pandemic in the univariate step of the DCC-GARCH estimation. The dates of the identified structural beaks are reported in Table 3, and the plots showing the structural breaks are reported in Figure 3. In these plots, the black lines are the volatilities estimated with a standard GARCH(1,1) model, the green lines indicate the structural breaks dates, and the red bars represent confidence intervals at 95% level for the dates considered. The plots show three main findings. First, every commodity experiences a strong rise in the volatility level in the midst of the first wave of the pandemic, spanning from March to July 2020. The second finding confirms the presence of structural breaks in all six commodities for the period considered. Third, Propane and Gasoline experienced additional structural breaks at the end of 2020. In particular, for Propane, a structural break in December 2020 marks a period of rising volatility that might be explained by two main factors. Given the main use of this commodity for domestic applications and space heating, it is natural to expect a spike in volatility in winter caused by the potential mismatch between supply and a higher demand of this commodity. Moreover, as reported in US Government Accountability Office (2003), a higher volatility in Propane returns might be traced back also to the limited capacity of the distribution system and the resulting bottlenecks in delivering the commodity. This issue, already present in the Propane delivering system, might have been worsened by the restrictions imposed during the second wave, which caused a spike in the demand of this commodity and a simultaneous negative impact on the delivery system. Gasoline instead experienced a structural break in October 2020, marking a period of low volatility that, differently from the other commodities, is maintained for the rest of the sample period. Moreover, the period of structural breaks common to all the six commodities coincides with a growth in deaths due to COVID-19 infections, as shown in Figure 2, confirming that this variable represents a good proxy for the COVID-19 pandemic trend.
The second part of the empirical application concerns the DCC models estimation and evaluation. In order to estimate the DCC model, the univariate models (GARCH, GJR, GM, and DAGM) are estimated. The estimated coefficients are reported in Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9.
For almost every commodity, the parameters related with the low-frequency variable are significant, meaning that the information conveyed by this variable improves the estimation of the commodities’ volatility. The correlation models (cDCC, DCC-MIDAS, DECO) are estimated using standardized residuals obtained from the univariate model’s estimation. The number of lagged realizations of COVID-19-related deaths in the US used in the long-run equation is K = 24 ; thus, the correlations are estimated starting from June 2020. The estimated coefficients for the DCC models are reported in Table 10.
In Table 11, we report the results of the MCS model based on three robust loss function, namely RMSE, FROB, and EUCL, measuring the distance between the estimated correlation matrix and the proxy of the covariance matrix, computed as the cross-products of r i , t . The models belonging to the SSM are highlighted in gray.
Except for the standard GARCH model entering in the SSM with the FROB loss function, the only models entering in the SSM are represented by the ones using the DAGM as the univariate model. Thus, the inclusion of the weekly low-frequency variable in the univariate model through a mixed-frequency approach improves the conditional covariance’s predictions.
Along with the results of the DCC models estimations, we also present a graphical representation of the correlations estimated with the models entering in the SSM. Plots in Figure 4, Figure 5 and Figure 6 report the time-varying pattern of the estimated correlations. Most of the commodities appear to be highly correlated regardless of the model used. In particular, the DCC-MIDAS model shows that, in general, the correlations increase over time, except for the correlation between Propane and the other commodities. In this case, the correlation reaches its maximum of 0.6 around the end of 2020 and then decreases to 0.3 at the end of the sample period. For the rest of the commodities, as well as for the estimates from the cDDC model (both with the GARCH and DAGM univariate specifications), the correlation is around 0.8. Again, the only exception is found for the Propane, whose correlation with other commodities ranges between 0.4 and 0.5.
Our results confirm that the COVID-19 pandemic increases the spillover effects and contagion among energy commodities, which, in turn, contribute to higher the level of systemic risk. These findings are in line with recent studies that show high correlations during the pandemic both in financial markets (Akhtaruzzaman et al. 2021; Chevallier 2020; Corbet et al. 2020; Yousfi et al. 2021; Zaremba et al. 2020; Zhang et al. 2020) and in energy markets (Hauser et al. 2020; Nyga-ukaszewska and Aruga 2020). Moreover, our results suggest that the DAGM specification used in the multivariate models allows us to obtain more accurate estimates of the correlations (as also found in Candila 2021b) among commodities, and this highlights the presence of an asymmetric effect of COVID-19 on energy prices, as shown in Wang and Su (2021).
Our analysis can be particularly useful in a risk management setting. The estimation of the matrix H i , t can be easily exploited to compute well-known risk measures. For example, models that possess good predictive ability of the (conditional) covariance matrix should also yield portfolios with lower variance with respect to less performing models (Engle and Kelly 2012). The evaluation of the resulting portfolios opens the door to the economic evaluation of the models (as done, for instance, by Amendola and Candila 2017). Moreover, the high correlation detected among the commodities provides useful information about their role in the portfolio diversification. Furthermore, multivariate GARCH models are widely studied to model the role of commodities for hedging financial risks (Chang et al. 2011; Chkili et al. 2014; Ku et al. 2007; Silvennoinen and Thorp 2013). Therefore, our methodology can be exploited as a valuable tool to mitigate commodity risks, especially during financial and economic recessions such as the 2008 crisis and the recent COVID-19 pandemic (Hauser et al. 2020). In particular, the inclusion of the number of weekly deaths caused by COVID-19 used to estimate the correlations allows us to discover the fragility regarding energy markets implied by the increased spillovers effects among commodities. This allows us to connect an increased exposure to spillovers effects from the energy markets to the financial markets. With this respect, our proposed approach can make the risk diversification during the pandemic more effective.

4. Conclusions

The increasing use of energy commodities for speculation and risk diversification purposes has resulted in significant shifts in prices of these assets and requires reliable estimates of their daily volatility and correlations. In this paper, we employed the Dynamic Conditional Correlation (DCC) model with a mixed-frequency approach in the univariate specifications to assess the effects of the COVID-19 pandemic in determining the correlation among energy markets. As in Candila (2021b), the DCC model with the GARCH-MIDAS and Double Asymmetric GARCH-MIDAS (DAGM) models as univariate specifications has been applied to six energy commodities, including the low-frequency variable of “weekly deaths in the US related to COVID-19 infection”. This approach allowed us to consider the information deriving from the low-frequency variable in order to estimate daily volatilities eventually driven by the COVID-19 data. Correlations have been estimated relying on three well-known models, namely the corrected DCC (cDCC), DCC-MIDAS, and Dynamic EquiCorrelation (DECO).
Our contribution has been two-fold. On one hand, we extended the current literature on the effects of the COVID-19 pandemic on the energy sector considering six energy commodities (West Texas Intermediate crude oil, Brent crude oil, Heating oil #2, Propane, New York Harbor Conventional Gasoline Regular, and Kerosene-Type Jet Fuel), whereas previous studies mainly focused on oil and natural gas (De Blasis and Petroni 2021; Gil-Alana and Monge 2020; Lin and Su 2021; Narayan 2020; Wang and Su 2021). On the other hand, we employed multivariate GARCH models in the energy sector (Chang et al. 2011; Chkili et al. 2014; Ku et al. 2007; Silvennoinen and Thorp 2013; Yousfi et al. 2021) within a risk management setting that, for the first time in this strand of literature, considers time series sampled at mixed-frequency. This allowed us to include the number of weekly deaths caused by the COVID-19 in the US for the estimation of daily volatilities and correlations. Overall, we highlighted an increase in the correlations among the commodities studied during the COVID-19 pandemic providing evidence of the contribution of the COVID-19 to the level of systemic commodity risk.
The empirical results showed that correlation models based on the DAGM model deliver better results than those that do not account for the MIDAS component. Moreover, the DAGM specification also provided better correlations estimates with respect to the GARCH-MIDAS model, which does not account for the different effects of positive and negative values of the low-frequency variable on the volatility of the commodities.
In conclusion, the analysis showed that our approach is statistically valuable and superior to other models in estimating correlations among energy commodities. The results matter in a risk management framework, where they can be exploited by investors for portfolio diversification. The results are also relevant to accomplish the requirements of both the microeconomic and macroeconomic regulators. At the microeconomic level, our estimates of the conditional covariance matrix can be exploited for obtaining reliable estimates of well-known risk measures. Moreover, the correct estimation of correlation in energy markets and volatilities has implications in terms of policy targeting for importing and exporting nations as well as in the evaluation of investments through asset pricing techniques. Our analysis suggests to conjecture increasing spillovers effect from the energy markets to the financial markets. This investigation is left for further studies.

Author Contributions

The authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the findings of this study are openly available at the US Energy Information Administration (EIA) website: https://www.eia.gov.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
USUnited States
WTIWest Texas Intermediate
DCCDynamic conditional correlation (Engle 2002)
cDCCCorrected DCC (Aielli 2013)
DECODynamic equicorrelation (Engle and Kelly 2012)
MIDASMIxing-Data Sampling
GMGARCH-MIDAS (Engle et al. 2013)
DAGMDouble Asymmetric GARCH-MIDAS model (Amendola et al. 2019)
MCSModel Confidence Set (Hansen et al. 2011)
SSMSet of Superior Models

Note

1
The literature has proposed alternative distributions to capture fat tails and skewness of returns, such as the multivariate Student’s t distribution. However, according to Pesaran and Pesaran (2010), the use of the M V N allows us to maintain the two-step process used to estimate the DCC model. Moreover, in order to take into consideration possible deviations from the normality assumption, all the reported standard errors are based on Quasi-Maximum Likelihood (Bollerslev and Wooldridge 1992) methods.

References

  1. Abid, Ilyes, Abderrazak Dhaoui, Stéphane Goutte, and Khaled Guesmi. 2020. Hedging and diversification across commodity assets. Applied Economics 52: 2472–92. [Google Scholar] [CrossRef]
  2. Aielli, Gian Piero. 2013. Dynamic conditional correlation: On properties and estimation. Journal of Business & Economic Statistics 31: 282–99. [Google Scholar]
  3. Akhtaruzzaman, Md, Sabri Boubaker, and Ahmet Sensoy. 2021. Financial contagion during COVID-19 crisis. Finance Research Letters 38: 101604. [Google Scholar] [CrossRef]
  4. Aloui, Chaker, and Samir Mabrouk. 2010. Value-at-risk estimations of energy commodities via long-memory, asymmetry and fat-tailed GARCH models. Energy Policy 38: 2326–39. [Google Scholar] [CrossRef]
  5. Amendola, Alessandra, Manuela Braione, Vincenzo Candila, and Giuseppe Storti. 2020. A model confidence set approach to the combination of multivariate volatility forecasts. International Journal of Forecasting 36: 873–91. [Google Scholar] [CrossRef]
  6. Amendola, Alessandra, and Vincenzo Candila. 2017. Comparing multivariate volatility forecasts by direct and indirect approaches. Journal of Risk 19: 33–57. [Google Scholar] [CrossRef]
  7. Amendola, Alessandra, Vincenzo Candila, and Giampiero M. Gallo. 2019. On the asymmetric impact of macro–variables on volatility. Economic Modelling 76: 135–52. [Google Scholar] [CrossRef]
  8. Arenas, Alex, Wesley Cota, Jesús Gómez-Gardeñes, Sergio Gómez, Clara Granell, Joan T. Matamalas, David Soriano-Paños, and Benjamin Steinegger. 2020. Modeling the spatiotemporal epidemic spreading of COVID-19 and the impact of mobility and social distancing interventions. Physical Review X 10: 041055. [Google Scholar] [CrossRef]
  9. Arouri, Mohamed El Hedi, Jamel Jouini, and Duc Khuong Nguyen. 2012. On the impacts of oil price fluctuations on European equity markets: Volatility spillover and hedging effectiveness. Energy Economics 34: 611–17. [Google Scholar] [CrossRef]
  10. Bai, Jushan, and Pierre Perron. 1998. Estimating and testing linear models with multiple structural changes. Econometrica 66: 47–78. [Google Scholar] [CrossRef]
  11. Basher, Syed Abul, and Perry Sadorsky. 2016. Hedging emerging market stock prices with oil, gold, VIX, and bonds: A comparison between DCC, ADCC and GO-GARCH. Energy Economics 54: 235–47. [Google Scholar] [CrossRef] [Green Version]
  12. Bauwens, Luc, Sébastien Laurent, and Jeroen V. K. Rombouts. 2006. Multivariate GARCH models: A survey. Journal of Applied Econometrics 21: 79–109. [Google Scholar] [CrossRef] [Green Version]
  13. Behmiri, Niaz Bashiri, Matteo Manera, and Marcella Nicolini. 2019. Understanding dynamic conditional correlations between oil, natural gas and non-energy commodity futures markets. The Energy Journal 40: 55–76. [Google Scholar] [CrossRef]
  14. Bollerslev, Tim. 1986. Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics 31: 307–27. [Google Scholar] [CrossRef] [Green Version]
  15. Bollerslev, Tim, and Jeffrey M. Wooldridge. 1992. Quasi-maximum likelihood estimation and inference in dynamic models with time-varying covariances. Econometric Reviews 11: 143–72. [Google Scholar] [CrossRef]
  16. Bouazizi, Tarek, Mongi Lassoued, and Zouhaier Hadhek. 2021. Oil Price Volatility Models during Coronavirus Crisis: Testing with Appropriate Models Using Further Univariate GARCH and Monte Carlo Simulation Models. International Journal of Energy Economics and Policy 11: 281–92. [Google Scholar] [CrossRef]
  17. Candila, Vincenzo. 2021a. Dccmidas: A Package for Estimating DCC-Based Models in R. R Package Version 0.1.0. Available online: https://www.researchgate.net/publication/350064227_dccmidas_A_package_for_estimating_DCC-based_models_in_R (accessed on 1 February 2021).
  18. Candila, Vincenzo. 2021b. Multivariate analysis of cryptocurrencies. Econometrics 9: 28. [Google Scholar] [CrossRef]
  19. Chang, Chia-Lin, Michael McAleer, and Roengchai Tansuchat. 2011. Crude oil hedging strategies using dynamic multivariate GARCH. Energy Economics 33: 912–23. [Google Scholar] [CrossRef] [Green Version]
  20. Chevallier, Julien. 2020. COVID-19 pandemic and financial contagion. Journal of Risk and Financial Management 13: 309. [Google Scholar] [CrossRef]
  21. Chevallier, Julien, and Florian Ielpo. 2013. Volatility spillovers in commodity markets. Applied Economics Letters 20: 1211–27. [Google Scholar] [CrossRef]
  22. Chkili, Walid, Chaker Aloui, and Duc Khuong Nguyen. 2014. Instabilities in the relationships and hedging strategies between crude oil and US stock markets: Do long memory and asymmetry matter? Journal of International Financial Markets, Institutions and Money 33: 354–66. [Google Scholar] [CrossRef]
  23. Colacito, Riccardo, Robert F. Engle, and Eric Ghysels. 2011. A component model for dynamic correlations. Journal of Econometrics 164: 45–59. [Google Scholar] [CrossRef] [Green Version]
  24. Corbet, Shaen, John W. Goodell, and Samet Günay. 2020. Co-movements and spillovers of oil and renewable firms under extreme conditions: New evidence from negative WTI prices during COVID-19. Energy Economics 92: 104978. [Google Scholar] [CrossRef] [PubMed]
  25. Corbet, Shaen, Charles Larkin, and Brian Lucey. 2020. The contagion effects of the COVID-19 pandemic: Evidence from gold and cryptocurrencies. Finance Research Letters 35: 101554. [Google Scholar] [CrossRef]
  26. Cuñado, Juncal, and Fernando Pérez de Gracia. 2003. Do oil price shocks matter? Evidence from some European countries. Energy Economics 25: 137–54. [Google Scholar] [CrossRef] [Green Version]
  27. De Blasis, Riccardo, and Filippo Petroni. 2021. Price leadership and volatility linkages between oil and renewable energy firms during the COVID-19 pandemic. Energies 14: 2608. [Google Scholar] [CrossRef]
  28. Denton, Michael, Adrian Palmer, Ralph Masiello, and Petter Skantze. 2003. Managing market risk in energy. IEEE Transactions on Power Systems 18: 494–502. [Google Scholar] [CrossRef]
  29. Engle, Robert F. 2002. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. Journal of Business & Economic Statistics 20: 339–50. [Google Scholar]
  30. Engle, Robert F., Eric Ghysels, and Bumjean Sohn. 2013. Stock market volatility and macroeconomic fundamentals. Review of Economics and Statistics 95: 776–97. [Google Scholar] [CrossRef]
  31. Engle, Robert F., and Bryan Kelly. 2012. Dynamic equicorrelation. Journal of Business & Economic Statistics 30: 212–28. [Google Scholar]
  32. Falkowski, Michał. 2011. Financialization of commodities. Contemporary Economics 5: 4–17. [Google Scholar] [CrossRef] [Green Version]
  33. Fan, Ying, Yue-Jun Zhang, Hsien-Tang Tsai, and Yi-Ming Wei. 2008. Estimating ‘Value at Risk’of crude oil price and its spillover effect using the GED-GARCH approach. Energy Economics 30: 3156–71. [Google Scholar] [CrossRef]
  34. Foglia, Matteo, and Eliana Angelini. 2020. Volatility connectedness between clean energy firms and crude oil in the COVID-19 era. Sustainability 12: 1–22. [Google Scholar] [CrossRef]
  35. Ghorbel, Achraf, and Ahmed Jeribi. 2021. Volatility spillovers and contagion between energy sector and financial assets during COVID-19 crisis period. Eurasian Economic Review 1: 1–9. [Google Scholar]
  36. Ghysels, Eric, Arthur Sinko, and Rossen Valkanov. 2007. MIDAS regressions: Further results and new directions. Econometric Reviews 26: 53–90. [Google Scholar] [CrossRef]
  37. Gil-Alana, Luis A., and Manuel Monge. 2020. Crude Oil Prices and COVID-19: Persistence of the Shock. Energy Research Letters 1: 1–4. [Google Scholar] [CrossRef]
  38. Glosten, Lawrence R., Ravi Jagannathan, and David E. Runkle. 1993. On the relation between the expected value and the volatility of the nominal excess return on stocks. The Journal of Finance 48: 1779–801. [Google Scholar] [CrossRef]
  39. Hamilton, James D. 1983. Oil and the Macroeconomy since World War II. Journal of Political Economy 91: 228–48. [Google Scholar] [CrossRef]
  40. Hansen, Peter R., Asger Lunde, and James M. Nason. 2003. Choosing the best volatility models: The model confidence set approach. Oxford Bulletin of Economics and Statistics 65: 839–61. [Google Scholar] [CrossRef]
  41. Hansen, Peter R., Asger Lunde, and James M. Nason. 2011. The Model Confidence Set. Econometrica 79: 453–97. [Google Scholar] [CrossRef] [Green Version]
  42. Hauser, Philipp, Carl-Philipp Anke, J. Gutiérrez López, Dominik Möst, Hendrik Scharf, David Schönheit, and S. Schreiber. 2020. The impact of the COVID-19 crisis on energy prices in comparison to the 2008 financial crisis. In IAEE Energy Forum/Covid-19 Issue. Cleveland: International Association for Energy Economics. [Google Scholar]
  43. Hung, Jui-Cheng, Ming-Chih Lee, and Hung-Chun Liu. 2008. Estimation of Value-at-Risk for energy commodities via fat-tailed GARCH models. Energy Economics 30: 1173–91. [Google Scholar] [CrossRef]
  44. IEA. 2020. Oil Market Report—April 2020. Paris: IEA. [Google Scholar]
  45. Iyke, Bernard N. 2020. COVID-19: The reaction of US oil and gas producers to the pandemic. Energy Research Letters 1: 1–7. [Google Scholar] [CrossRef]
  46. Kang, Wensheng, Ronald A. Ratti, and Kyung Hwan Yoon. 2015. The impact of oil price shocks on the stock market return and volatility relationship. Journal of International Financial Markets, Institutions and Money 34: 41–54. [Google Scholar] [CrossRef] [Green Version]
  47. Karali, Berna, and Octavio A. Ramirez. 2014. Macro determinants of volatility and volatility spillover in energy markets. Energy Economics 46: 413–21. [Google Scholar] [CrossRef]
  48. Krehbiel, Tim, and Lee C. Adkins. 2005. Price risk in the Nymex energy complex: An extreme value approach. Journal of Futures Markets: Futures, Options, and Other Derivative Products 25: 309–37. [Google Scholar] [CrossRef]
  49. Ku, Yuan-Hung Hsu, Ho-Chyuan Chen, and Kuang-hua Chen. 2007. On the application of the dynamic conditional correlation model in estimating optimal time-varying hedge ratios. Applied Economics Letters 14: 503–9. [Google Scholar] [CrossRef]
  50. Laporta, Alessandro G., Luca Merlo, and Lea Petrella. 2018. Selection of Value at Risk models for energy commodities. Energy Economics 74: 628–43. [Google Scholar] [CrossRef]
  51. Laurent, Sébastien, Jeroen V.K. Rombouts, and Francesco Violante. 2013. On loss functions and ranking forecasting performances of multivariate volatility models. Journal of Econometrics 173: 1–10. [Google Scholar] [CrossRef] [Green Version]
  52. Lin, Boqiang, and Tong Su. 2021. Does COVID-19 open a Pandora’s box of changing the connectedness in energy commodities? Research in International Business and Finance 56: 101360. [Google Scholar] [CrossRef]
  53. Lv, Xiaodong, and Xian Shan. 2013. Modeling natural gas market volatility using GARCH with different distributions. Physica A: Statistical Mechanics and its Applications 392: 5685–99. [Google Scholar] [CrossRef]
  54. Masters, Michael W., and Adam K. White. 2008. The Accidental Hunt Brothers–Act 2: Index Speculators Have Been a Major Cause of the Recent Drop in Oil Prices, Masters Capital Management and White Knight Research and Trading. Index Speculators Have Been a Major Cause of the Recent Drop in Oil Prices. Special Update September 10: 4. [Google Scholar]
  55. Mishra, Amritkant. 2019. Crude oil, stock market, and foreign exchange return volatility and spillover: A GARCH DCC analysis of Indian and Japanese financial market. International Journal of Business Innovation and Research 20: 25–46. [Google Scholar] [CrossRef]
  56. Naeem, Muhammad A., Saqib Farid, Safwan M. Nor, and Syed J.H. Shahzad. 2021. Spillover and drivers of uncertainty among oil and commodity markets. Mathematics 9: 441. [Google Scholar] [CrossRef]
  57. Narayan, Paresh K. 2020. Oil price news and COVID-19—Is there any connection? Energy Research Letters 1: 13176. [Google Scholar] [CrossRef]
  58. Nouvellet, Pierre, Sangeeta Bhatia, Anne Cori, Kylie E. C. Ainslie, Marc Baguelin, Samir Bhatt, Adhiratha Boonyasiri, Nicholas F. Brazeau, Lorenzo Cattarino, Laura V. Cooper, and et al. 2021. Reduction in mobility and COVID-19 transmission. Nature Communications 12: 1–9. [Google Scholar] [CrossRef] [PubMed]
  59. Nyga-ukaszewska, Honorata, and Kentaka Aruga. 2020. Energy prices and COVID-immunity: The case of crude oil and natural gas prices in the US and Japan. Energies 13: 6300. [Google Scholar] [CrossRef]
  60. Pesaran, Bahram, and M. Hashem Pesaran. 2010. Conditional volatility and correlations of weekly returns and the VaR analysis of 2008 stock market crash. Economic Modelling 27: 1398–416. [Google Scholar] [CrossRef] [Green Version]
  61. Sadorsky, Perry. 1999. Oil price shocks and stock market activity. Energy Economics 21: 449–69. [Google Scholar] [CrossRef]
  62. Sadorsky, Perry. 2003. The macroeconomic determinants of technology stock price volatility. Review of Financial Economics 12: 191–205. [Google Scholar] [CrossRef]
  63. Sadorsky, Perry. 2012. Correlations and volatility spillovers between oil prices and the stock prices of clean energy and technology companies. Energy Economics 34: 248–55. [Google Scholar] [CrossRef]
  64. Sadowski, Adam, Zbigniew Galar, Robert Walasek, Grzegorz Zimon, and Per Engelseth. 2021. Big data insight on global mobility during the covid-19 pandemic lockdown. Journal of Big Data 8: 1–33. [Google Scholar] [CrossRef]
  65. Silvennoinen, Annastiina, and Timo Teräsvirta. 2009. Multivariate GARCH Models. Handbook of Financial Time Series; Berlin and Heidelberg: Springer, pp. 201–29. [Google Scholar]
  66. Silvennoinen, Annastiina, and Susan Thorp. 2013. Financialization, crisis and commodity correlation dynamics. Journal of International Financial Markets, Institutions and Money 24: 42–65. [Google Scholar] [CrossRef] [Green Version]
  67. US Government Accountability Office (GAO). 2003. Propane: Causes of Price Volatility, Potential Consumer Options, and Opportunities to Improve Consumer Information and Federal Oversight. Washington: US Government Accountability Office. [Google Scholar]
  68. Wang, Kai-Hua, and Chi-Wei Su. 2021. Asymmetric link between COVID-19 and fossil energy prices. Asian Economics Letters 1: 18742. [Google Scholar]
  69. Wang, Yudong, and Chongfeng Wu. 2012. Forecasting energy market volatility using GARCH models: Can multivariate models beat univariate models? Energy Economics 34: 2167–81. [Google Scholar] [CrossRef]
  70. Yousfi, Mohamed, Abderrazak Dhaoui, and Houssam Bouzgarrou. 2021. Risk spillover during the COVID-19 global pandemic and portfolio management. Journal of Risk and Financial Management 14: 222. [Google Scholar] [CrossRef]
  71. Zaremba, Adam, Renatas Kizys, David Y. Aharon, and Ender Demir. 2020. Infected Markets: Novel Coronavirus, Government Interventions, and Stock Return Volatility around the Globe. Finance Research Letters 35: 101597. [Google Scholar] [CrossRef]
  72. Zhang, Dayong, Min Hu, and Qiang Ji. 2020. Financial markets under the global pandemic of COVID-19. Finance Research Letters 36: 101528. [Google Scholar] [CrossRef]
Figure 1. Log−returns of commodities. Notes: Plots of log-returns of commodities Sample period: April 2020–June 2021. Number of observations: 276.
Figure 1. Log−returns of commodities. Notes: Plots of log-returns of commodities Sample period: April 2020–June 2021. Number of observations: 276.
Risks 09 00144 g001
Figure 2. Weekly deaths related to COVID-19 infections in the US. Notes: Plot of the weekly deaths related to COVID-19 infections in the US. Sample period: March 2020–June 2021. The labels of the x-axis indicate the week and the year of the observations. The time series starts the first week of March 2020 (9 March 2020).
Figure 2. Weekly deaths related to COVID-19 infections in the US. Notes: Plot of the weekly deaths related to COVID-19 infections in the US. Sample period: March 2020–June 2021. The labels of the x-axis indicate the week and the year of the observations. The time series starts the first week of March 2020 (9 March 2020).
Risks 09 00144 g002
Figure 3. Volatility structural breaks. Notes: The figure shows the GARCH(1,1) estimated volatility (black lines), the structural breaks (green lines), and the related 95% confidence intervals (red lines).
Figure 3. Volatility structural breaks. Notes: The figure shows the GARCH(1,1) estimated volatility (black lines), the structural breaks (green lines), and the related 95% confidence intervals (red lines).
Risks 09 00144 g003
Figure 4. The plots show the estimated volatilities on the main diagonal and correlation estimated with the model GARCH-cDCC. Sample period: June 2020–June 2021.
Figure 4. The plots show the estimated volatilities on the main diagonal and correlation estimated with the model GARCH-cDCC. Sample period: June 2020–June 2021.
Risks 09 00144 g004
Figure 5. The plots show the estimated volatilities on the main diagonal and correlation estimated with the model DAGM-cDCC. Sample period: June 2020–June 2021.
Figure 5. The plots show the estimated volatilities on the main diagonal and correlation estimated with the model DAGM-cDCC. Sample period: June 2020–June 2021.
Risks 09 00144 g005
Figure 6. The plots show the estimated volatilities on the main diagonal and the correlation estimated with the model DAGM-DCC-MIDAS. Red lines represent the estimated long-run correlation. Sample period: June 2020 to June 2021.
Figure 6. The plots show the estimated volatilities on the main diagonal and the correlation estimated with the model DAGM-DCC-MIDAS. Red lines represent the estimated long-run correlation. Sample period: June 2020 to June 2021.
Risks 09 00144 g006
Table 1. Functional forms of DCC correlation models.
Table 1. Functional forms of DCC correlation models.
ModelFunctional Form
cDCC (Aielli 2013) H i , t = D i , t R i , t D i , t
D i , t = d i a g ( σ i , t , 1 , , σ i , t , j , , σ i , t , n )
R i , t = ( diag ( Q i , t ) ) 1 / 2 Q i , t ( diag ( Q i , t ) ) 1 / 2
Q i , t = ( 1 a b ) Ψ + a ( ξ i 1 , t ξ i 1 , t ) + b Q i 1 , t
ξ i , t = D i , t 1 r i , t
Ψ = E ξ i , t ξ i , t
DCC-MIDAS (Colacito et al. 2011) H i , t = D i , t R i , t D i , t
D i , t = d i a g ( σ i , t , 1 , , σ i , t , j , , σ i , t , n )
R i , t = ( diag ( Q i , t ) ) 1 / 2 Q i , t ( diag ( Q i , t ) ) 1 / 2
Q i , t = ( 1 a b ) R ¯ i , t ( ω ) + a ( ξ i 1 , t ξ i 1 , t ) + b Q i 1 , t
ξ i , t = D i , t 1 r i , t
DECO (Engle and Kelly 2012) H i , t = D i , t R i , t D E C O D i , t
D i , t = d i a g ( σ i , t , 1 , , σ i , t , j , , σ i , t , n )
R i , t D E C O = ( 1 ρ i , t ) I n + ρ i , t J n
ρ i , t = 1 n ( n 1 ) ι R i , t ι n
Notes: The table shows the different specifications of the correlation models used in this paper. R ¯ i , t ( ω ) is the long-run correlation of Colacito et al. (2011), J n and ι are, respectively, a matrix and a ( n × 1 )-vector of ones.
Table 2. Summary statistics.
Table 2. Summary statistics.
Min.Max.MeanSDSkew.Kurt.
WTI−0.0820.2460.0060.0332.4515.718
Europe Brent−0.0870.2220.0060.0291.58211.723
Heating Oil−0.1390.1120.0040.029−0.0043.686
Propane−0.0820.1720.0040.030.7854.604
Gasoline−0.0870.1140.0050.0280.1312.039
Kerosene−0.130.1530.0050.0330.5533.889
Notes: The table presents the main statistics (the minimum (Min.) and maximum (Max.), the mean, standard deviation (SD), skewness (Skew.), and excess kurtosis (Kurt.)) for the close-to-close log-returns. Sample period: March 2020–June 2021. Number of observations: 271.
Table 3. Structural break dates of the six commodities.
Table 3. Structural break dates of the six commodities.
WTIBrent OilHeating OilPropaneGasolineKerosene
Mar. 2020Mar. 2020Mar. 2020Mar. 2020Mar. 2020Mar. 2020
June 2020June 2020June 2020June 2020June 2020June 2020
Dec. 2020Oct. 2020
Table 4. The estimated coefficients of the GARCH model.
Table 4. The estimated coefficients of the GARCH model.
EstimateStd. Errort ValuePr (>|t|)Sig.
WTI
ω 0.0000.0002.0570.040**
α 0.1380.0443.1690.002***
β 0.7860.05314.8140.000***
Europe Brent
ω 0.0000.0002.3160.021**
α 0.1380.0532.5910.010***
β 0.8070.04119.6460.000***
Heating Oil
ω 0.0000.0002.7440.006***
α 0.1360.0433.1330.002***
β 0.7930.04617.4080.000***
Propane
ω 0.0000.0000.7360.462
α 0.0690.0501.3780.168
β 0.8670.1227.1200.000***
Gasoline
ω 0.0000.0001.5580.119
α 0.1000.0551.8190.069*
β 0.8560.06613.0050.000***
Kerosene
ω 0.0000.0002.4780.013**
α 0.1080.0293.7860.000***
β 0.8260.04120.2000.000***
Notes: The table reports the estimates of the GARCH model. The column “Std. Error” shows the Quasi-Maximum Likelihood-based standard errors. Sample Period: March 2020–June 2021, 271 daily observations. *, **, and *** indicate the significance at levels 10%, 5%, and 1%, respectively.
Table 5. The estimated coefficients of the GJR model.
Table 5. The estimated coefficients of the GJR model.
EstimateStd. Errort ValuePr (>|t|)Sig.
WTI
ω 0.0000.0002.3360.019**
α 0.1040.0442.3800.017**
β 0.7970.04019.7050.000***
γ 0.1010.0771.3030.193
Europe Brent
ω 0.0000.0002.3080.021**
α 0.1330.0681.9560.050*
β 0.8070.04119.7060.000***
γ 0.0180.0750.2380.812
Heating Oil
ω 0.0000.0002.5430.011**
α 0.1250.0562.2440.025**
β 0.7810.05215.0830.000***
γ 0.0570.0910.6240.532
Propane
ω 0.0000.0000.9190.358
α 0.0460.0421.1030.270
β 0.8760.05914.7880.000***
γ 0.0990.0811.2280.219
Gasoline
ω 0.0000.0001.2170.224
α 0.1000.0601.6750.094*
β 0.8840.07611.6480.000***
γ −0.0550.047−1.1690.242
Kerosene
ω 0.0000.0001.9190.055*
α 0.0890.0362.4670.014**
β 0.8350.04219.6960.000***
γ 0.0580.0680.8600.390
Notes: The table reports the estimates of the GJR model. The column “Std. Error” shows the Quasi-Maximum Likelihood-based standard errors. Sample Period: March 2020–June 2021, 271 daily observations. *, **, and *** indicate the significance at levels 10%, 5%, and 1%, respectively.
Table 6. The estimated coefficients of the DAGM model for WTI, Europe Brent, and Heating Oil.
Table 6. The estimated coefficients of the DAGM model for WTI, Europe Brent, and Heating Oil.
EstimateStd. Errort ValuePr (>|t|)Sig.
WTI
α 0.0001.6630.0001.000
γ 0.2110.1371.5460.122
β 0.8482.0260.4190.675
m−8.3831.518−5.5240.000***
θ + 4.6443.8011.2220.222
ω 2 + 1.0511.2230.8590.390
θ −14.01411.070−1.2660.206
ω 2 2.6987.5870.3560.722
Europe Brent
α 0.0000.1900.0011.000
γ 0.2460.2181.1300.258
β 0.8750.1874.6770.000***
m−5.1642.723−1.8960.058*
θ + −1.10721.944−0.0500.960
ω 2 + 5.2269.2510.5650.572
θ −8.34016.254−0.5130.608
ω 2 5.0782.5921.9590.050*
Heating Oil
α 0.0000.0430.0020.998
γ 0.0870.0302.8740.004***
β 0.9510.05018.9120.000***
m−7.8650.669−11.7540.000***
θ + 2.7611.1952.3100.021**
ω 2 + 1.0010.5971.6780.093*
θ −11.2383.981−2.8230.005***
ω 2 3.5801.5852.2580.024**
Notes: The table reports the estimates of the DAGM model for WTI, Europe Brent, and Heating Oil. The column “Std. Error” shows the Quasi-Maximum Likelihood-based standard errors. Sample Period: March 2020–June 2021, 271 daily observations. *, **, and *** indicate the significance at levels 10%, 5%, and 1%, respectively.
Table 7. The estimated coefficients of the GM model for WTI, Europe Brent, and Heating Oil.
Table 7. The estimated coefficients of the GM model for WTI, Europe Brent, and Heating Oil.
EstimateStd. Errort ValuePr (>|t|)Sig.
WTI
α 0.0000.0590.0020.999
γ 0.3200.1033.1000.002***
β 0.8380.0849.9470.000***
m−4.0090.788−5.0880.000***
θ −3.3321.459−2.2830.022**
ω 2 5.7151.9842.8810.004***
Europe Brent
α 0.0000.0330.0030.998
γ 0.2200.0892.4750.013**
β 0.8880.06413.8420.000***
m−4.8550.916−5.3010.000***
θ −4.9021.454−3.3720.001***
ω 2 5.5991.4153.9580.000***
Heating Oil
α 0.0790.2260.3500.726
γ 0.4710.0964.9000.000***
β 0.6850.2602.6360.008***
m−2.7990.533−5.2550.000***
θ −3.2004.227−0.7570.449
ω 2 2.0921.5621.3390.180
Notes: The table reports the estimates of the GM model for WTI, Europe Brent, and Heating Oil. The column “Std. Error” shows the Quasi-Maximum Likelihood-based standard errors. Sample Period: March 2020–June 2021, 271 daily observations. **, and *** indicate the significance at levels 10%, 5%, and 1%, respectively.
Table 8. The estimated coefficients of the GM model for Propane, Gasoline, and Kerosene.
Table 8. The estimated coefficients of the GM model for Propane, Gasoline, and Kerosene.
EstimateStd. Errort ValuePr (>|t|)Sig.
Propane
α 0.0500.0520.9700.332
γ 0.1500.0781.9200.055*
β 0.8610.06712.8010.000***
m−5.9661.141−5.2290.000***
θ −0.8512.370−0.3590.720
ω 2 3.1391.2052.6060.009***
Gasoline
α 0.0000.0280.0040.997
γ 0.0620.0252.5070.012**
β 0.9680.03130.7500.000***
m−6.3170.261−24.1590.000***
θ −4.4780.876−5.1150.000***
ω 2 7.0421.0546.6800.000***
Kerosene
α 0.0400.0890.4450.656
γ 0.1780.0503.5280.000***
β 0.8690.0969.0610.000***
m−4.8110.581−8.2770.000***
θ −1.1590.929−1.2470.212
ω 2 23.798612.7920.0390.969
Notes: The table reports the estimates of the GM model for Propane, Gasoline, and Kerosene. The column “Std. Error” shows the Quasi-Maximum Likelihood-based standard errors. Sample Period: March 2020–June 2021, 271 aily observations. *, **, and *** indicate the significance at levels 10%, 5%, and 1%, respectively.
Table 9. The estimated coefficients of DAGM model for Propane, Gasoline, and Kerosene.
Table 9. The estimated coefficients of DAGM model for Propane, Gasoline, and Kerosene.
EstimateStd. Errort ValuePr (>|t|)Sig.
Propane
α 0.1830.1791.0220.307
γ 0.0150.4720.0320.974
β 0.4220.6630.6370.524
m−7.6920.399−19.2770.000***
θ + 2.0245.9730.3390.735
ω 2 + 1.0017.0150.1430.887
θ −7.4482.914−2.5560.011**
ω 2 1.0015.8490.1710.864
Gasoline
α 0.0000.0240.0040.997
γ 0.0640.0302.1600.031**
β 0.9670.02144.9810.000***
m−6.4980.349−18.6130.000***
θ + −2.6951.786−1.5090.131
ω 2 + 6.3971.4064.5500.000***
θ −6.1362.476−2.4780.013**
ω 2 7.8031.2166.4160.000***
Kerosene
α 0.0000.1110.0010.999
γ 0.0750.0531.4010.161
β 0.9460.1645.7550.000***
m−7.9290.549−14.4400.000***
θ + 3.5861.3782.6020.009***
ω 2 + 1.0010.8381.1940.232
θ −12.2613.098−3.9570.000***
ω 2 1.9300.8192.3570.018**
Notes: The table reports the estimates of the DAGM model for Propane, Gasoline, and Kerosene. The column “Std. Error” shows the Quasi-Maximum Likelihood-based standard errors. Sample Period: March 2020 to June 2021, 271 daily observations. **, and *** indicate the significance at levels 10%, 5%, and 1%, respectively.
Table 10. Correlation estimates.
Table 10. Correlation estimates.
GARCHGJRGMDAGM
cDCCa0.141***0.155 ***0.139 ***0.132 ***
(0.045)(0.045)(0.039)(0.033)
b0.579 ***0.586 ***0.45 ***0.473 ***
(0.097)(0.092)(0.154)(0.156)
DCCMIDASa0.0210.0230.0230.025 **
(0.022)(0.023)(0.017)(0.012)
b0.956 ***0.955 ***0.958 ***0.961 ***
(0.024)(0.025)(0.018)(0.014)
ω 2 1.0011.0011.0011.001
(0.747)(0.652)(0.718)(0.896)
DECOa0.216 ***0.227 ***0.197 *0.165 **
(0.069)(0.072)(0.118)(0.071)
b0.439 ***0.424 ***0.3530.286
(0.155)(0.15)(0.281)(0.415)
Notes: The table reports the estimated coefficients of the correlation models (first column) according the univariate specifications reported in columns from three to six. Numbers in parentheses represent the Quasi-Maximum Likelihood-based standard errors. *, **, and *** indicate the significance at levels 10%, 5%, and 1%, respectively.
Table 11. A set of Superior Models according to the Model Confidence Set.
Table 11. A set of Superior Models according to the Model Confidence Set.
Univ.GARCHGARCHGARCHGJRGJRGJRGMGMGMDAGMDAGMDAGM
Mult.DCCMcDCCDECODCCMcDCCDECODCCMcDCCDECODCCMcDCCDECO
RMSE5.7445.7995.8595.755.8175.8745.8275.935.9715.6555.7365.784
FROB3.383.273.353.373.273.353.383.323.383.273.193.27
EUCL2.262.202.242.252.202.242.272.252.282.202.162.19
Notes: The table reports the average losses multiplied by 100,000. Label “DCCM” stands for DCC-MIDAS. Shades of gray denote inclusion in the SSM at significance level α = 0.25 .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Andreani, M.; Candila, V.; Morelli, G.; Petrella, L. Multivariate Analysis of Energy Commodities during the COVID-19 Pandemic: Evidence from a Mixed-Frequency Approach. Risks 2021, 9, 144. https://doi.org/10.3390/risks9080144

AMA Style

Andreani M, Candila V, Morelli G, Petrella L. Multivariate Analysis of Energy Commodities during the COVID-19 Pandemic: Evidence from a Mixed-Frequency Approach. Risks. 2021; 9(8):144. https://doi.org/10.3390/risks9080144

Chicago/Turabian Style

Andreani, Mila, Vincenzo Candila, Giacomo Morelli, and Lea Petrella. 2021. "Multivariate Analysis of Energy Commodities during the COVID-19 Pandemic: Evidence from a Mixed-Frequency Approach" Risks 9, no. 8: 144. https://doi.org/10.3390/risks9080144

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