Next Article in Journal
Scale Mixture of Maxwell-Boltzmann Distribution
Previous Article in Journal
Power System Transient Stability Assessment Based on Machine Learning Algorithms and Grid Topology
Previous Article in Special Issue
Risk-Taking, Financial Knowledge, and Risky Investment Intention: Expanding Theory of Planned Behavior Using a Moderating-Mediating Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantile Dependence between Crude Oil Returns and Implied Volatility: Evidence from Parametric and Nonparametric Tests

1
Department of Business Administration, College of Business Administration in Majmaah, Majmaah University, Majmaah 11952, Saudi Arabia
2
Faculty of Economics and Management in Nabeul, University of Carthage, Carthage 2085, Tunisia
3
BESTMOD Laboratory, Higher Institute of Management of Tunis, University of Tunis, Tunis 1002, Tunisia
4
School of Business, Lebanese American University, Beirut P.O. Box 13-5053, Lebanon
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(3), 528; https://doi.org/10.3390/math11030528
Submission received: 8 December 2022 / Revised: 7 January 2023 / Accepted: 16 January 2023 / Published: 18 January 2023
(This article belongs to the Special Issue Application of Mathematical Methods in Financial Economics)

Abstract

:
We examine the daily dependence and directional predictability between the returns of crude oil and the Crude Oil Volatility Index (OVX). Unlike previous studies, we apply a battery of quantile-based techniques, namely the quantile unit root test, the causality-in-quantiles test, and the cross-quantilogram approach. Our main results show evidence of significant bi-directional predictability that is quantile-dependent and asymmetric. A significant positive Granger causality runs from oil (OVX) returns to OVX (oil) returns when both series are in similar lower (upper) quantiles, as well as in opposite quantiles. The Granger causality from OVX returns to oil returns is only significant during periods of high volatility, although it is not always positive. The findings imply that the forward-looking estimate of oil volatility, reflecting the sentiment of oil market participants, should be considered when studying price variations in the oil market, and that crude oil returns can be used to predict oil implied volatility during bearish market conditions. Therefore, the findings have implications regarding predictability under various conditions for oil market participants.

1. Introduction

Large price fluctuations in the crude oil market have become a staple of the last two decades, subsequent to changes in supply–demand fundamentals [1], global economic activity [2], geopolitical risk [3], and various crisis periods such as the COVID-19 outbreak [4,5]. On 20 April 2020, crude oil prices experienced huge and unprecedented turbulence, during which the price of US benchmark crude oil West Texas Intermediate (WTI) futures plunged to below USD 0 a barrel for the first time in history. This was particularly due to two main factors. The first is the slowdown in the world economy following the rapid spread of the COVID-19 pandemic, which eventually spawned a sharp decrease in crude oil demand and an oversupply of crude oil. The second is the conflict between OPEC and non-OPEC producers on the oil production level, including the oil price war between Saudi Arabia and Russia.
Oil price volatility is a major concern for consumers, corporations, and governments. This is not surprising given that oil price fluctuations potentially have considerable effects on the investment decisions of investors and corporations as well as the global economy [6,7,8,9]. Usually, volatility is computed based on historical price data using realized volatility estimators, GARCH-based models, and stochastic volatility models. Another way to compute volatility is via the so-called implied volatility, based on data from the options market. Several studies indicate that implied volatility is superior to the former models [10,11,12,13]. In this regard, one of the indicators that investors and analysts regularly follow to assess the future state of crude oil volatility is the Chicago Board Options Exchange (CBOE) Crude Oil Volatility Index (OVX), which was introduced in 2007. OVX represents the 30-day forward-looking market estimate of volatility and is viewed as a barometer of investor sentiment about the crude oil market. Compared to historical volatility, it offers a more timely assessment of risk as it comprises both past market information and future anticipation [13,14]. Accordingly, OVX can reasonably transmit trading signals for market operators because it allows market participants to identify and link the movement of OVX to overall market confidence and panic in the crude oil market.
Interestingly, the risk-return trade-off is a crucial component of investment decisions and portfolio management. Despite the importance of this relationship, there is no clear agreement on its direction among researchers and analysts. Previous studies mainly concentrate on the relationship between returns and historical volatility in financial markets [15,16,17,18,19], while providing mixed results. Some show that the relationship is positive [20,21,22], while others find evidence of a negative relationship [23,24,25]. With the emergence of implied volatility and its role as a barometer of risk assessment [13,14], much attention has been paid to the return–implied volatility link in financial markets [26]. However, such research in the context of the oil market is relatively scarce. As noted by Ji and Fan [27], recommendations deduced in the stock market cannot be easily adopted to OVX in virtue of the more complex structure of crude oil options. Existing studies using OVX consider its predictive power for oil return volatility (e.g., [10,28,29,30,31,32,33]), or the spillover effect between oil and equity markets based on implied volatility indexes (e.g., [34,35,36,37]), whereas limited studies concentrate on the relationships between crude oil returns and OVX. In this regard, most studies examine the linear relationship between oil returns and OVX (e.g., [38,39,40]), whereas fewer consider the non-linear dependence. Agbeyegbe [41], Ji and Fan [27], Liu et al. [21], and Echaust and Just [42] use approaches based on copula models; Li et al. [43] apply a multiscale multivariate multifractal detrended fluctuation analysis; and Shaikh [44], Silva Junior [45], and Fousekis [46] use approaches based on quantile regression methods. Nevertheless, none of these studies assess the Granger causal relationship between crude oil returns and OVX based on quantile-based techniques accounting for various market conditions. Moreover, the disagreement in outcomes regarding the relationship between the two metrics indicates that the issue is still largely open to empirical investigation. In this paper, we examine the dependence and directional predictability between the returns of crude oil and OVX using daily data and various quantile-based techniques. Taking daily data from 10 May 2007 to 17 April 2020, we apply a battery of quantile-based techniques, specifically the quantile unit root test developed by Galvao [47], the causality-in-quantiles test of Troster [48], and the cross-quantilogram of Han et al. [49].
This paper contributes to the existing literature in the following ways. Firstly, to the best of our knowledge, this is the first study to cast new light on the Granger causal interaction between crude oil returns and changes in implied volatility using the quantile unit root test developed by Galvao [47], the causality-in-quantiles test of Troster [48], and the cross-quantilogram of Han et al. [49]. These tests can show whether variables are nonstationary or exhibit a heterogeneous causal structure across quantiles. Particularly, the test of Troster [48] has the merit of taking into account potential nonlinear causal links in conditional quantiles of the distribution, while the approach of Han et al. [49] is advantageous in at least two aspects. Notably, it is less prone to misspecification and takes into account large lags in the causality test. These approaches, based on a non-linear econometric framework, fit perfectly in our context, given that the causality links between oil prices and OVX are potentially non-linear and specific to market conditions and volatility states. Secondly, our analysis allows us to investigate not merely the magnitude but also the sign, duration, potential asymmetry, and Granger causal quantile of the dependence between the variables under consideration. Thirdly, the existence of causality in tails advocates the utility and significance of our modelling approach, and therefore casts doubt on the suitability of the conventional conditional mean-founded causality test, which gives only a partial picture by relying on the conditional mean.
The main findings reveal the existence of significant bi-directional predictability between oil returns and OVX. Notably, the relationship is found to be quantile-dependent and asymmetric across various market conditions. In fact, a significant positive dependence running from WTI (OVX) returns to OVX (WTI) is seen when both oil returns and oil implied volatility are in similar lower (upper) quantiles. Evidence of dependence is also reported when they are in opposite quantiles; for example, WTI returns Granger-cause OVX during extreme and moderate market conditions. For the reverse case, the relationship appears to be significant only during periods of high levels of OVX, but is not always positive. Furthermore, the quantile dependence from WTI returns to OVX is different from that from OVX to WTI returns, suggesting an asymmetry. Henceforth, investor sentiment in the oil market, as reflected by OVX, can be considered a driver of forthcoming variations in the returns in the oil market, and WTI returns affect OVX.
The rest of the paper is structured as follows. Section 2 provides a brief review of the related literature. Section 3 presents the econometric methodology. Section 4 defines the data and offers a preliminary investigation. Section 5 presents and discusses the empirical results. Section 6 concludes.

2. Literature Review

How the expected return of an investment is linked to the expected risk is widely documented in the academic literature [50]. Risk-return trade-off is an essential component of investment decisions and assessing portfolios. Interestingly, with the emergence of implied volatility [13,14], much attention has been paid to the return–implied volatility link in financial markets [26], but less evidence exists on the implied volatility in the context of the oil market.
The existing literature on OVX can be divided into three main strands. The first strand focuses on the predictive power of OVX for oil return volatility. For example, Szakmary et al. [10] find that the implied volatility (IV) is superior to the historical volatility (HV) as a predictor of the afterward realized volatility. They note that the IV embeds a significant amount of information that should be taken into consideration. Consequently, the authors propose a combined version of GARCH and IV models to improve oil volatility predictions. Likewise, Dutta [30] proves that OVX has a better predictive power for oil volatility than realized volatility (RV) measures. By incorporating GARCH models in the analysis, Agnolucci [51] shows that these models, mainly with Generalized Error Distribution (GED) residuals, outperform implied volatility (IV). In the same context, Lux et al. [29] evaluate the forecasting performance of various competing volatility models (GARCH with linear and non-linear dynamics, implied volatility, and multifractal models (MSM)), and show the superiority of MSM models over the other ones, across various forecasting horizons and sub-periods. This superiority is equally maintained in predicting the value-at-risk (VaR). Chen et al. [32] reveal that integrating OVX-based implied volatility into GARCH-type models generally ameliorates the predictive performance of volatility in crude oil market. Using the heterogeneous autoregressive-realized volatility (HAR-RV) models, Haugom et al. [28] show that econometric models can be enhanced by integrating implied volatility. Utilizing similar models, Lv [31] show evidence of the significant and positive effect of OVX on future realized volatility. Furthermore, the impact of a large OVX on future volatility is larger than a small OVX. Recently, Benedetto et al. [33] apply an entropy-based method and find a surge in the information flow between OVX and the spot variance of Brent returns and a corresponding drop in the information flow with WTI. Nevertheless, the causal flow running from the oil spot variance to OVX is larger in magnitude and more statistically significant for WTI than for Brent.
Previous studies also examine the impact of crude oil shocks on the job market [52,53], GDP [54], and the oil-food nexus [55]. On a different front, Ma et al. [56] forecast the realized volatility of oil futures. Notably, another strand of empirical studies focuses on the spillover effect between crude oil and equity markets based on implied volatility indexes. For example, Maghyereh et al. [34] investigate the directional connectedness between oil and equities based on implied volatility indexes. The authors prove the existence of two-way information flows between the two markets, although it is stronger from the oil market to equity markets than the other way around. Bouri et al. [35] use the ARDL cointegration technique and the causality approach and suggest the existence of strong nonlinear relationships between the implied volatilities of gold and oil and that of the Indian stock market. However, their outcomes indicate no evidence of feedback link from the Indian stock market to the gold and crude oil markets. Similarly, Choi and Hong [36] use autoregressive distributed lag (ARDL) bounds and Toda–Yamamoto Granger causality tests to show that OVX and VIX exhibit a two-way causality in the period that covers the shale gas revolution, but no causal link in the period that does not. They confirm that OVX Granger-causes VKOSPI in the pre-shale gas revolution phase but fail to ascertain any causal flow between the two variables in the post-revolution phase. Recently, Naeem et al. [37] investigate the dynamic connectedness between a set of implied volatility indexes (oil (OVX), gold (GVZ), currency (EVZ), and bond (TYNVI)), and conclude that OVX is the leading driver of the inherent VIX connectedness of the global equity markets during anxiety phases.
The third strand of literature concentrates on the relationships between crude oil and OVX returns. For example, Aboura and Chevallier [38] assess the role exercised by implied volatility on the WTI crude oil price and establish that the so-called “inverse leverage effect”, which means that volatility is rising more following positive returns than following negative returns, is the main effect driving WTI crude oil price evolutions. Specifically, they argue that this stylized evidence could reflect the anxiety of oil consumers facing rising crude oil prices. Chen et al. [39] employ a Kalman filter with time varying coefficients in the regression equations and find a feeble negative link between OVX changes and future crude oil price returns, and that extreme values (high or low) of OVX cannot correctly forecast future evolutions of returns (positive or negative). Using copula quantile regressions, Agbeyegbe [41] show an inverted U-shaped relationship between oil returns and OVX, suggesting a stronger connection at the tails of the distribution than at the central part. Similarly, Ji and Fan [27] investigate the causal links between WTI returns and OVX and find evidence of a strong unidirectional causal link in-mean (in-variance) from WTI returns (OVX changes) to OVX changes (WTI returns). They also show that the contemporaneous relation between them is negative and asymmetric, which discloses that OVX manifests as a gauge of investor fear than as risk preference. Liu et al. [21] use copula models and find that WTI returns and OVX are mostly negatively linked. They note the presence of significant risk flows from OVX to WTI. In an analogous manner, Silva Junior [45] studies the relationship between United States Oil Fund (USO) returns and OVX, employing ordinary least squares (OLS), quantile regressions (QR), and the nonparametric B-splines regression model. The results show a negative, asymmetric and nonlinear contemporary relationship between the two series, that the effects of negative shocks in volatility are more obvious than positive shocks, and, excitingly, the association has no homogeneous behavior across the various quantiles. Utilizing nonparametric quantile models, Fousekis [46] reveals that the association between oil price returns and OVX is contemporaneous, negative, and asymmetric, and indicates the presence of an inverted U-shaped dependence between the two measures, indicating that the pricing of implied volatility is higher for large fluctuations than for small ones, and are lighter for large positive variations than for large negative ones. Lin and Tsai [40] examine the long-run dynamics between oil prices and OVX. They reveal that increasing OVX exhibits a larger negative influence on oil prices than declining OVX, indicating the existence of a long-run asymmetric cointegrating relationship between them. It is worth mentioning that this result of asymmetry is also recognized in previous findings [41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57]. Furthermore, Shaikh [44] uses neural network and quantile approaches and suggests that crude oil prices are aligned with OVX. They find an asymmetric relationship and show evidence that the implied volatility is likely to be calm during global financial crises and increases after the crisis periods. Echaust and Just [42] study the tail dependence behavior of WTI oil returns and OVX changes. The results suggest that the strongest tail dependence of negative oil price shocks and OVX changes happens during the pandemic crisis. They establish the ability of the OVX index to be a sentiment indicator for the crude oil market. Nevertheless, they fail to ascertain its capability to ameliorate value-at-risk estimates. Li et al. [43] explore the characteristics of the cross-correlations between the crude oil market and OVX. The findings reveal that the cross-correlated comportments of smaller changes are persistent, although large changes are anti-persistent between crude oil and OVX. Furthermore, the results confirm that oil price variations and OVX change in a positive direction in the short-term, although they change in the inverse direction in the long-term.
More recently, Raggad [58] applies the Cross-Quantilogram approach to explore whether the OVX can reliably predict WTI returns, and finds that the information content of OVX plays a role in forecasting WTI returns once it is in the higher conditional quantile, but no evidence of predictability is uncovered once the OVX is at the low to medium quantile levels.
As we can see from the above review, the existing literature does not examine the directional predictability patterns between oil returns and oil implied volatility, based on a quantile-based methodology that can potentially enrich our understanding of the sign, magnitude, duration, and asymmetric Granger causal interactions between oil returns and oil implied volatility across upper, middle, and lower quartiles. This is where we aim to contribute.

3. Methodology

The first sub-section briefly introduces the quantile autoregressive (QAR) unit root tests of Koenker and Xiao [59] and Galvao [47]. The second sub-section presents the Granger causality test in the quantiles of Troster [48]. The third sub-section describes the nonparametric cross-quantilogram approach of Han et al. [49].

3.1. Quantile Unit Root Test

Researchers give attention to quantile autoregressive (QAR) unit root tests when examining the stationarity properties of time series. A major reason for the use of this kind of test is its ability to investigate the stationarity property of data across the quantiles of the conditional distribution beyond the limitation of the conditional mean [60]. This test is of immense importance for accurately identifying the dynamic behavior and distributed characteristics of the time series data [61]. QAR unit root tests were introduced by Koenker and Xiao [59] and later generalized by Galvao [47], who explicitly incorporated covariates and a linear time trend into the QAR specification.
To formalize the idea empirically, suppose that the time series process X t is strictly stationary and let t = X t 1 , , X t s d be an information set based on the history of Y t . Let us denote by F X . t X the conditional distribution function of X t given t X , so the QAR unit root test is specified by a quantile linear regression model. Specifically, the QAR model is described as:
Q τ X X t t X = λ 1 τ + λ 2 τ t + β τ X t 1 + j = 1 p β j τ Δ X t j + F u 1 τ
where Q τ X . t Y is the τ-quantile of F X . t X ,   λ 1 τ   defines the drift term, t   designates the linear time trend, β τ represents the persistence measure, and F u 1 denotes the inverse conditional distribution function of the errors based on various quantiles τ, τ Γ ⊂ [0,1].
By estimating Equation (1) for various quantiles, we obtain the estimation of the persistence measure β ^ corresponding to each quantile of the conditional distribution Y t . The null hypothesis of the QAR test is expressed as: H 0 : β τ = 1 , and its respective t-test is calculated through the t-statistic suggested by Koenker and Xiao [59] and Galvao [47] for the entire range of quantiles τ Γ .

3.2. Granger Causality in Quantiles

As noted, the test of Troster [48] has the merit of taking into account potential nonlinear causal relations in the conditional quantiles of the return distribution. It is sufficiently flexible to test for Granger causality in specific regions of the distribution, such as the center or tails (left or right).
To formally present the test, let Y t , Z t ,   t Z   be a jointly stationary and ergodic time series defined in some probability space Ω , A , P and Y t be a dependent variable; and let J t   d , d = s + q, be an explanatory vector, given by: J t J t Y , J t Z , where J t Y : = Y t 1 , , Y t s s   and J t Z : = Z t 1 , , Z t q q   .
The null hypothesis for the Granger causality test assumes that Z t does not Granger-cause Y t   and is given by:
H 0 Z Y : F Y y J t Y , J t Z = F Y y J t Y   f o r   a l l   y     ,      
where F Y . J t Y , J t Z   and F Y . J t Y   designate the conditional distribution function of Y t given J t Y , J t Z and J t Y , respectively. The null hypothesis defined by Equation (2) is called Granger non-causality in distribution.
Given the difficulty in estimating the conditional distribution, the formulation defined by Equation (2) is limited to the test for Granger non-causality in mean, which is barely a necessary condition for the expression given by Equation (3) [48]. It follows that Z t does not Granger-cause Y t in mean if:
E Y t J t Y , J t Z = E Y t J t Y   f o r   a l l   y            
where E Y t J t Y , J t Z and E Y t J t Y   represent the means of F Y . J t Y , J t Z and F Y . J t Y , respectively.
As noted by Troster [48], the preceding specifications have two major drawbacks. Firstly, the uniform causality in mean fails to accommodate potential links in the extreme conditional quantiles of the distribution. Secondly, the null hypothesis described in Equation (3) tells us nothing of the “level of the causality” in the case of rejecting the hypothesis. To overcome these two limitations, Troster [48] proposes Granger non-causality in conditional quantiles.
By substituting F Y . J t Y , I t Z and F Y y J t Y   by their τ-quantiles, noted by Q τ Y , Z . J t Y , J t Z and Q τ Y , Z . J t Y , respectively, in Equation (2), the following expression can be obtained:
H 0 Q C ,   Z Y : Q τ Y , Z y J t Y , J t Z = Q τ Y y J t Y   ,   a . s .   f o r   a l l   τ     Γ ,    
where Γ is a compact set such that Γ ⊂ [0,1], and the conditional τ-quantiles of Y t fulfil the subsequent constraints (conditions):
P r Y t Q τ Y y J t Y J t Y τ   f o r   a l l   τ     Γ
P r Y t Q τ Y , Z y J t Y , J t Z J t Y , J t Z τ   f o r   a l l   τ     Γ
With an explanatory vector J t , we get P r Y t Q τ Y y J t Y J t Y = E 1 Y t Q τ Y y J t Y J t Y , where 1 Y t y     denotes an indicator function of the event that Y is less than or equal to y. Therefore, Equation (4) can be written in its equivalent form:
E { 1 [ Y t Q τ Y , Z y J t Y , J t Z J t Y , J t Z ] } = E { 1 [ Y t Q τ Y y J t Y J t Y ] }   ,   a . s .   f o r   a l l   τ     Γ
where the left-hand side of Equation (5) defines the τ-quantile of F Y . J t Y , J t Z   by definition. Following Troster [48], a parametric model is required for the estimation of the τ-th quantile of F Y . J t Y . To this end, suppose that Q Y . J t Y is appropriately modelled by a parametric quantile specification s . , φ τ , which belongs to a set of functions S = s . , φ τ φ . : τ φ τ Φ p ,   f o r   τ   Γ 0 , 1 .
Next, under the null hypothesis given by Equation (5), the τ-conditional quantile, Q Y . J t Y , is appropriately modelled by a parametric quantile specification s   I t Y , φ 0 τ defined on the subset information J t Y . It follows that the null hypothesis of absence of Granger causality is reformulated as follows:
H 0 Q C ,   Z Y : E 1 Y t s J t Y , φ 0 τ J t Y , J t Z = τ   ,   a . s .   f o r   a l l   τ     Γ ,  
Versus
H A Q C ,   Z Y : E 1 Y t s J t Y , φ 0 τ J t Y , J t Z τ   ,   a . s .   f o r   a l l   τ     Γ ,    
where   s I t Y , φ 0 τ appropriately defines the true conditional quantile Q Y . I t Y , for all τ   Γ . Equation (6) can be reformulated as H 0 Q C ,   Z Y : E 1 ( Y t m I t Y , φ 0 τ 0 τ I t Y , I t Z = 0 , for all   τ     Γ . Furthermore, the null hypothesis of Equation (6) can be described by a sequence of unconditional moment restrictions such as:
E 1 ( Y t s I t Y , φ 0 τ 0 τ exp i ω I t = 0   ,   for   all     τ     Γ  
where exp i ω I t exp i ω 1 Y t 1 Z t 1 + + i ω r Y t r Z t r represents a weighting function for all ω r   with r d , and i = 1 is the imaginary root. Hence, a parametric test statistic can be formulated as the sample analogue of this expression   E 1 ( Y t s I t Y , φ 0 τ 0 τ exp i ω I t :
v T ω , τ = 1 T t = 1 T 1 ( Y t s I t Y , φ T τ 0 τ exp i ω I t                        
where φ T is a T -consistent estimator of φ 0 τ , for all τ     Γ . Then, the following test statistic suggested by Troster [48] can be deduced:
z T : = T W v T ω , τ 2 d F ω ω d F τ τ
where the function F ω . designates the conditional distribution of a d-variate standard normal vector, F τ . is a uniform discrete distribution defined on the set of T in n equally spaced random points, T n =   τ j j = 1 n , and ω   d is the vector of weights obtained from a standard normal distribution. The expression of Equation (10) can be assessed by its sample equivalent. Consider a T × n matrix with components ψ i j = Ψ τ j Y i s I t Y , φ T τ j , where Ψ τ j . is the function Ψ τ j ε 1 ε 0 τ j . The test is therefore based on the statistic:
z T : = 1 T n j = 1 n ψ . j W ψ . j
where W is the n × n matrix with components s w t , s = exp ( 0.5 J t J s 2 ) and   ψ . j   designates the j-th column of Ψ. The rejection of the null hypothesis of the test is conditioned by observing large values of z T in Equation (11).
Following the pioneering work of Troster [48], the subsampling procedure is applied to compute critical values for z T in Equation (11). The procedure therefore consists of computing z T for the whole subsample, and the p-values are deduced from the average value of the B subsamples test statistics.
Notably, the subsample is selected of size b = k T 2 / 5 , where . stands for the integer part of a number, and k is a constant term [62].
To obtain the z T test, as given in Equation (11), we follow Troster [48] and choose three QAR specifications s . , for different quantiles of order τ , τ   Γ 0 , 1 , under the null hypothesis of absence of Granger causality in Equation (8) as:
Q A R 1 :   s 1 J t Y , φ τ = λ 1 τ + λ 2 τ Y t 1 + β t F u 1 τ        
Q A R 2 :   s 2 J t Y , φ τ = λ 1 τ + λ 2 τ Y t 1 + λ 3 τ Y t 2 + β t F u 1 τ  
Q A R 3 :   s 3 J t Y , φ τ = λ 1 τ + λ 2 τ Y t 1 + λ 3 τ Y t 2 + λ 4 τ Y t 3 + β t F u 1 τ  
where the parameters φ τ = λ 1 τ , λ 2 τ , λ 3 τ , λ 4 τ , β t are estimated by a maximum likelihood method, and F u 1 denotes the inverse of a standard normal distribution.

3.3. The Cross-Quantilogram

We investigate the directional quantile dependence and spillover effect between WTI returns and OVX changes across quantiles by utilizing the cross-quantilogram method of Han et al. [49]. This method has various attractive features. Firstly, it captures the correlation and causal dependence between the two variables, irrespective of the distribution region. Secondly, it is free of any distributional assumptions regarding the variables under study to test causality. Thirdly, it allows us to deal with long lags in the Granger causality test.
Formally, the method is presented as follows. For any two-dimensional strictly stationary time series x i , t ,   t Z ,   i = 1 , 2 , we denote by F i . the cumulative distribution of x i , t with density functions f i . for i = 1 , 2 . The respective quantile of the series x i , t is given by q i α i = i n f u :   F i u   α i for α i   0 , 1 . In a bivariate setting, the quantiles of order α = α 1 , α 2 are given by q 1 α 1 , q 2 α 2 . To detect the quantile dependence between the two series, the cross-quantilogram, CQ, for the α-quantile with k lags can be found by applying the formula:
ρ k α = E Ψ α 1 x 1 , t q 1 α 1 Ψ α 2 x 2 , t k q 2 ,   t k α 2 E Ψ α 1 2 x 1 , t q 1 α 1   E Ψ α 2 2 x 2 , t k q 2 , t k α 2
For k = 0, ±1, ±2, …, where Ψ α u 1 u < 0 α ,   the term 1[·] indicates the indicator function, while 1[ x i , t q i α i ] is recognized as the quantile hit or quantile exceedance process. It follows that the sample form of CQ, based on observed data, is defined as:
ρ ^ α k = t = k + 1 T Ψ α 1 x 1 , t q ^ 1 α 1 Ψ α 2 x 2 , t k q ^ 2 , t k α 2 t = k + 1 T Ψ α 1 2 x 1 , t q ^ 1 α 1   t = k + 1 T Ψ α 2 2 x 2 , t k q ^ 2 , t k α 2
Note that the q ^ i α i parameter in Equation (16) designates the unconditional sample quantile of x i , t , as specified by Han et al. [49]. Furthermore, Han et al. [49] suggest a quantile form of the Ljung–Box–Pierce test statistic with the null hypothesis H0:   ρ ^ α k = 0 for all k 1 , , p . H1: ρ ^ α k 0 for all k 1 , , p :
Q ^ α p = T T + 1 k = 1 p ρ ^ α 2 k T k
where Q ^ α p represents the portmanteau test of directional predictability from one time series to another corresponding to p lags and quantile pair α = α 1 , α 2 . For the construction of confidence intervals for the cross-quantilogram, Han et al. [49] suggest using the stationary bootstrap inference procedure and a self-normalized approach.

4. Data and Preliminary Analysis

4.1. Data Description

We use WTI crude oil prices and OVX. WTI is a global benchmark for crude oil prices in the US. The sample period spans 10 May 2007 to 17 April 2020 (yielding 3376 daily observations), based on the availability of OVX data sourced from CBOE. Data on crude oil prices are collected from the online database of the US Energy Information Administration (EIA). Our empirical analysis employs daily WTI and OVX returns computed as log-difference R t = 100 ln ( S t S t 1 ) ; where S t = W T I , or OVX. Figure 1 shows the time evolution of the daily levels and returns series over the period of investigation, and Table 1 gives the summary statistics of the daily return series.
As shown in Figure 1, both level and return series exhibit various significant fluctuations that usually concur with major geopolitical, fundamental, or economic events that affect the oil market, such as the 2008 Global Financial Crisis (GFC), the 2009–2012 euro-zone debt crisis, the Libyan crisis of 2011, the 2014–2016 Oil Price Slump, and the oil market turmoil related to the COVID-19 health crisis. The descriptive statistics of daily returns shown in Table 1 indicate that neither WTI returns nor OVX returns are normally distributed, as shown by the excess kurtosis and evidence that the return distributions are leptokurtic with asymmetric tails. The Jarque-Bera test confirms the rejection of the null hypothesis of normality for both return series. This result is validated by the Q-Q normal plot of the series given in Figure 2, because the quantile points deviate from the reference line.
While the OVX returns have a positive skewness, the WTI returns exhibit negative skewness, demonstrating a fatter tail on the left side (losses) of the distribution. The correlation between the two variables is significantly negative (−0.325). To test for serial correlation, we use the Ljung-Box statistic. The result reveals that the null hypothesis of no autocorrelation up to 20 lags is rejected and proves the existence of conditional heteroskedasticity in the dynamics of both series. As a reaction to the influx of new information, the market adjusts the volatility forecast and, accordingly, options prices.
It is worth noting that, as a part of our preliminary investigation, we assess the presence of non-linearity by performing the Brock–Dechert-Scheinkman (BDS) test of Broock et al. [63] on the residuals of the oil and OVX returns equations in the vector autoregressive (VAR) (1) model. As Table 2 shows, the null hypothesis of i.i.d. residuals is rejected at different embedding dimensions (m), providing evidence of non-linearity in the data.
We also check for the existence of multiple structural breaks by conducting the methodology developed by Bai and Perron [64]. Unreported results provide evidence of non-linearity and the presence of at least one structural break in both WTI and OVX return series.

4.2. Unit Root Testing

We begin by examining the stationarity properties of our time series. This step is a key requirement for performing our statistical methodology. Using three familiar unit root tests, specifically the Augmented Dickey and Fuller [65] (ADF), the Phillips and Perron [66] (PP), and the Zivot and Andrews [67] (ZA) tests, Table 3 shows that the return series are stationary at the 1% significance level.
As a part of our full quantile-based statistical methodology, we also include in our analysis the QAR unit root test to check whether or not the unit root results persist across the quantiles of the conditional distribution, not only at the conditional mean. As noted by Koenker and Xiao [59], this kind of test has higher power than conventional unit root tests. Furthermore, the quantile-based unit root test outperforms conventional unit root tests in case of a deviation from the hypothesis of Gaussian residuals. The results of the quantile unit root test are given in Table 4, which shows the persistence measure estimates as well as the t-statistics at the 5% significance levels, under the null hypothesis H0: β(τ) = 1, for various quantiles τ ∈ [0.05,0.95].
Both return series are stationary, irrespective of the level of significance, for all the considered quantiles of the conditional distribution. These findings support the previous results obtained with the conventional unit root tests (Table 3). Having this in hand, we proceed to the implementation of the quantile-based methodology.

4.3. Testing for Cointegration

We conduct the linear cointegration test of Johansen [68,69] and the quantile cointegration test suggested by Xiao [70] to check whether the series are cointegrated. Panel A of Table 5 indicates that there exists a linear cointegration between the series at the 5% significance level. Panel B, the results of the quantile cointegration test proposed by Xiao [70], shows evidence of a significantly non-linear cointegration relationship between the quantiles of both series at any level of significance, indicating that the cointegration link between the variables changes over the distribution.
Overall, the non-normality of the unconditional distribution, the non-linearity, and the presence of structural breaks in our data provide preliminary support for the use of the causality in quantiles test and the cross-quantilogram approach.

5. Empirical Results

5.1. Granger Causality Test Results

To explore the Granger causal relationship between WTI returns and OVX returns, we begin by conducting the standard test of Granger causality, for the purpose of comparison. The results of the Granger causality test between WTI and OVX, with three lag-lengths, k = 1, 2, 3, are given in Table 6. They show that the null of no Granger causality existing from WTI returns to OVX changes, cannot be rejected at any level of significance. However, OVX returns are found to Granger-cause in mean WTI returns at the 5% significance level. Notably, these results should be taken with caution because the standard test of Granger causality does not account for potential non-linearity.
After verifying the stationarity of our data, we are justified in the application of the parametric test of Troster [48] to examine the Granger causality in quantiles. While the standard Granger causality test shows evidence of a mean unidirectional causal relationship running from OVX changes to WTI returns, it is important to evaluate the validity of such a relationship at each of the quantiles considered. To this end, we follow Troster [48] and estimate three different QAR models, given by Equations (11)–(13), for each dependent variable in the z T test, given by Equation (11). Following Sakov and Bickel [62], a subsample is chosen of size b = k n 2 / 5 , where . corresponds to the integer part of a number, the constant parameter is given by k = 5, and n = 3373 is the sample size. The subsequent subsample size is b = 129 , but the unreported results show that the findings are not very sensitive to other choices of subsample size.
With this quantile approach test, we can distinguish between causality affecting the various portions of the conditional distribution (i.e., lower, median, or upper quantiles). Table 7 shows the subsample p-values of the z T test defined by Equation (11). WTI returns oil prices Granger-cause OVX returns at the 1% significance level, for all quantiles considered, except some quantiles in the middle part of the distribution (i.e., 0.4 and 0.45). For the reverse direction, there is evidence of Granger causality from OVX returns to WTI returns. However, we fail to confirm evidence for Granger causality either at the extreme quantiles or in the near centrally located quantiles of the distribution. Notably, these results are not sensitive to lag-length defining the QAR specification, under the null hypothesis of no Granger causality. These findings are partially in line with the result of Lin et al. [71], who report that oil prices have a bi-directional relation with OVX. Similarly, in a different application, Liu et al. [17] find a significant bi-directional implied volatility transmission between oil and US stock markets.
While the test of Troster [48] provides some insight into the causal relationships in conditional quantiles, it has one major drawback, which is common for parametric approaches, which is that it might be subject to ambiguous results owing to misspecification and consequently may lose its full ability to capture the true causal relationships between the return series. Additionally, the test of Troster [48] does not indicate the extent or sign of the relationship for the variables in question. Subsequently, to overcome these potential problems and for reasons of comparability, we apply the nonparametric cross-quantilogram approach of Han et al. [49].

5.2. Cross-Quantilogram Analysis

The results of the CQ are graphic, and in each plot we report the sample cross-quantilograms for lag k = 1, 2, …, 30, in conjunction with the 95% bootstrapped confidence intervals for no directional predictability with 1000 bootstrapped replicates, distinguished by red-dashed lines. The findings for the relationship between WTI and OVX returns are given in Figure 3, Figure 4, Figure 5 and Figure 6, when both markets are in the same quantiles, i.e., for the grid of 19 quantiles τ = 0.05 ;   0.95 , and in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 when both markets are in opposite quantiles. Note that if ρ k α = 0 , there is no dependence or directional predictability from an event x 2 , t k ,     q 2 , t k τ 2 : k = 1 , , p to an event x 1 , t k ,   q 1 , t k τ 1   : k = 1 , , p . In contrast, if ρ k α 0 , there is a quantile dependence or directional predictability between the two events. In our framework, the lead/lag parameter k, defined in terms of days, governs the delay in the predictability from one return time series to another.

5.2.1. Results When Both Return Series Are in the Same Quantiles

Figure 3 and Figure 4 show the cross-quantilogram and Box–Ljung test statistics for the directional predictability from WTI returns to OVX returns. The results show that WTI returns have a positive and significant effect on OVX returns when both are at their lower quantiles (less than 0.4). However, in general, this OVX quantile reliance on WTI is dispelled for larger lags, with a few non-substantial exceptions for some lags, suggesting that the relationship weakens over time. On the other hand, there is no sufficient evidence of dependence between the two return series when they are at intermediate to higher quantiles, except for some quantiles such as 0.45 and 0.5. We note that the quantile dependence is asymmetric as the negative spillover is stronger than the positive spillover.
In the case of the spillover from OVX to WTI, as shown in Figure 5 and Figure 6, OVX exhibits a positive and significant impact on WTI when both return series are at the upper quantiles. We, however, do not see any sufficient evidence of causality when the two return series are in lower to intermediate quantiles. These findings indicate that oil prices surge as OVX, which includes investor sentiment or a fear component, surges. In other words, when there is increasing uncertainty in the crude oil market, as reflected by an increasing OVX, both crude oil prices and OVX tend to leap upwards together. This result is in line with previous studies such as Jagerson (2008), who postulates that OVX can forecast oil prices, and Aboura & Chevallier [38], who establish that a surge in the oil price is followed by a surge in volatility. The results support the findings of Liu et al. [21], who show that WTI returns are more sensitive to OVX upsurges than OVX declines. Equally, we note that the quantile dependence is asymmetric as the spillover between right tails is higher than the spillover between left tails.
Based on the above analysis, we reveal the quantile dependence between the two return series when both are in the same quantiles, and generally find that OVX (WTI) tends to co-crash (boom) with WTI (OVX). Furthermore, the dependence between both return series is essentially found to be significant when both return series are in the tails of the distribution.
We next examine the interdependence between OVX and WTI returns when the returns series are in opposite quantiles.

5.2.2. Results When the Return Series Are in Opposite Quantiles

To capture the potential quantile dependence and directional predictability of one series to the other when the return series are in opposite quantiles, we present the cross-quantilogram of lag k = 1, 2, …, 30 when the distributions of the series returns are in bear, bull, and normal states, i.e., three quantiles are considered: low (τ = 0.1), high (τ = 0.9), and median (τ = 0.5). The full results for other quantiles and the Box–Ljung test statistic results are available from the authors upon request.
We first provide the results (Figure 7) of the cross-quantilogram when the WTI returns are in the lower quantile, i.e., q 2 α 2 for α 2 = 0.1 , and the OVX returns are in the other 19 quantiles, i.e., q 1 α 1 for α 1   ϵ   0.05 ,   0.10 , ,   0.95 . OVX is found to have a dependence of negative and positive signs on WTI. This outcome, though, is mainly significant in the lower quantiles, suggesting that the structure of dependence between the two series is mostly seen during the low OVX state. Next, we consider the scenario when WTI returns are in the median quantile, i.e., q 2 α 2 for α 2 = 0.5 , with positive (negative) and significant cross-quantilogram estimates, for some lags, at the upper (lower) quantiles of OVX returns (see Figure 8). When WTI returns are lower than the median, OVX returns have a positive (negative) dependence on WTI returns in the lower (upper) quantiles. Finally, we present the upper quantile scenario, when oil returns are in the higher quantile, q 2 α 2 = 0.9 , and OVX returns are in the other 19 quantiles, i.e., q 1 α 1 for α 1   ϵ   0.05 ,   0.10 , ,   0.95 . In Figure 9, we see evidence of spillovers from WTI returns to OVX returns. Overall, there is no consistency in either the sign or the pattern of the cross-quantilogram estimates. Furthermore, a significant negative dependence is only seen for some lower quantiles of OVX returns. This suggests that the relationship between the two return series is very weak. Interestingly, a comparison of Figure 7, Figure 8 and Figure 9 reveals that the predictive ability of WTI returns over OVX returns is usually more noticeable during bearish market conditions.
In the case where OVX returns are in the lower quantile, i.e., q 2 α 2 for α 2 = 0.1 and WTI returns are the other 19 quantiles, i.e., q 1 α 1 for α 1   ϵ   0.05 ,   0.10 , ,   0.95 , the results are given in Figure 10. Overall, the sample cross-quantilograms are statistically insignificant, suggesting that no significant dependence exists between WTI returns and OVX returns when the latter is in lower quantiles. Similarly, once the OVX returns are at the medium quantile (Figure 11), there is no sufficient indication justifying the quantile association between the return series. This advocates that, once OVX is in a normal condition, OVX returns do not have the ability to predict WTI returns.
Figure 12 shows that when OVX is in the higher quantile, q 2 α 2 = 0.9 , the cross-quantilogram estimates are negative (positive) and statistically significant for some lags in the lower (upper) quantiles. This outcome indicates that when OVX returns are very high, revealing rising uncertainty in the crude oil market, OVX returns exhibit a substantial impact on WTI returns. However, this impact is only found in the extreme quantiles (very bearish or very bullish conditions) of oil price returns. A comparison of Figure 10, Figure 11 and Figure 12 shows that the forecasting ability of OVX returns over WTI returns is usually more pronounced during soaring investor sentiment about the state of the crude oil market.
To sum up, the results indicate that when WTI and OVX return series are in different quantiles, WTI returns can be useful to predict extreme low/high variations in OVX returns during bearish market conditions. Indeed, they show that the lower quantile of WTI returns positively and significantly influences the quantiles of OVX returns, but the dependence is reversed during bear market conditions. Under normal market circumstances, the dependence can be either positive or negative. These results are quite comparable to Benedetto et al. [33], who present evidence of predictability from OVX to WTI, although our results are more nuanced and account for upper, middle, and lower quantiles. For the reverse relationship, OVX returns are found to be helpful in predicting extreme fluctuations in WTI returns. Indeed, evidence indicates that high OVX returns have significant and negative effects on the lower quantiles of WTI returns, but the dependence turns positive at upper quantiles. Furthermore, during low and middle OVX returns, the evidence does not support any directional predictability, which validates the results of Liu et al. [21], who show that oil price returns are reacting more to OVX increase than to OVX decrease. This also confirms the findings of Agbeyegbe [41] and Ji and Fan [27], who postulate that OVX is favored as a gauge of investor fear over risk preference.
Broadly, our findings prove that the relationships between OVX and oil prices exhibit heterogeneous behavior across quantiles and asymmetric dependence. These findings are partly in line with Agbeyegbe [41], Chen and Zou [57], Silva Junior [45], Lin and Tsai [40], Shaikh [44], and Fousekis [46], who show evidence of asymmetric relationships between OVX and oil prices. Furthermore, our findings are partly in accordance with Aboura and Chevallier [38], who establish the existence of an inverse leverage effect, which characterizes a larger increase in volatility following an upsurge in oil prices than following a decrease in oil prices of the same magnitude. However, the findings do not agree with Chen et al. [39], who find that a feeble negative association exists between OVX changes and future crude oil price returns, and that extreme levels of OVX cannot appropriately forecast future fluctuations of returns. However, it is worth noting that the studies mentioned above do not rely on quantile approaches as we do, particularly in relation to causal links between the two return series. Our comprehensive analysis allows us to cover the detailed interplay between oil prices and oil implied volatility.

6. Conclusions

In this study, we analyze the Granger causal interplay and directional predictability between crude oil returns and the returns of the Crude Oil Volatility Index (OVX), under various market conditions. To this end, we apply various quantile-based tests and models, specifically the quantile unit root test developed by Galvao [47], the parametric test of Granger causality in quantiles of Troster [48], and the nonparametric cross-quantilogram approach of Han et al. [49].
While the Granger causality in quantiles analysis indicates a significant bi-directional predictability between oil returns and OVX returns, the results of the cross-quantilogram approach show a complete and more nuanced picture of the relationships at various quantiles. Specifically, the relationship exhibits heterogeneous behavior, which points to asymmetry across quantiles. Furthermore, the quantile dependence from oil returns to OVX returns is different from that from OVX returns to WTI returns. The predictive power of oil returns is generally seen when both return series are in similar lower quantiles, which means that when the oil market is bearish, it is possible for OVX to be declining. In opposite quantiles, however, there is evidence of quantile predictability at the middle quantile and at the tails of oil returns, particularly the left tail.
For the reverse relationship, OVX returns can be used to predict extreme low/high variations in WTI prices. Precisely, when the level of OVX is high, OVX returns have a negative and significant effect on the lower quantiles of oil returns. However, this dependence turns positive at the upper quantiles of oil returns. At middle and low OVX return quantiles, there is no evidence of directional predictability. The existence of causal links in tails advocates the adequacy and significance of our modelling strategy, which casts doubt on the suitability of the conventional conditional mean-founded causality models that give a partial picture of the WTI-OVX relationship.
The findings of the study are useful for investors in the oil market, as they offer a new perspective on the relationship between return and volatility in the market and point to the asymmetric dependence across quantiles. Although several previous studies present evidence of a significant relationship between oil returns and OVX, the innovation of our study is that it reveals that the causality is quantile-dependent and asymmetric. As a result, investors aiming to control the evolution of prices should take the level of implied volatility into account, because the level of investor sentiment in the oil market, measured by OVX, may drive upcoming variations in spot crude oil returns at some quantiles. At the same time, spot oil price fluctuations can affect implied volatility levels under specific conditions.
The results are also relevant to portfolio risk management and the hedging strategies of participants in the crude oil market. For example, investors can better forecast the return and volatility of the oil market, conditional on market states, and thus refine their trading and hedging strategies and assess their exposure to oil market risk. Additionally, they can anticipate variations in their oil portfolio values due to oil volatility shocks, and can thus adjust their investment and trading positions, considering crude oil market states such as extremely unstable or stable markets. More precisely, investors should consider an asymmetric hedging strategy under upside and downside variations in the levels of OVX, particularly extreme movements of OVX upset by oil-related events that increase the uncertainty and complexity of the international oil market.
Given that the quantile dependence between crude oil prices and oil implied volatility may depend on multitude exogenous variables, future studies could extend this research by accounting for the effect of other variables such as US stock market conditions.

Author Contributions

B.R. conceptualized the paper, investigated and provided all resources, defined the research methodology, performed the estimation and analysis, and wrote the original draft preparation. E.B. supervised the project, contributed to the interpretation of the results, and revised the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

Data Availability Statement

The datasets used and analyzed during the study are available from public websites.

Conflicts of Interest

The authors declare no competing interest.

References

  1. Aguilera, R.F.; Eggert, R.G.; CC, G.L.; Tilton, J.E. Depletion and the future availability of petroleum resources. Energy J. 2009, 30, 141–174. [Google Scholar] [CrossRef]
  2. Demirbas, A.; Omar Al-Sasi, B.; Nizami, A.S. Recent volatility in the price of crude oil. Energy Sources Part B Econ. Plan. Policy 2017, 12, 408–414. [Google Scholar] [CrossRef]
  3. Joo, Y.C.; Park, S.Y. Oil prices and stock markets: Does the effect of uncertainty change over time? Energy Econ. 2017, 61, 42–51. [Google Scholar] [CrossRef]
  4. Bouri, E.; Demirer, R.; Gupta, R.; Pierdzioch, C. Infectious diseases, market uncertainty and oil market volatility. Energies 2020, 13, 4090. [Google Scholar] [CrossRef]
  5. Dutta, A.; Bouri, E.; Noor, M.H. Climate bond, stock, gold, and oil markets: Dynamic correlations and hedging analyses during the COVID-19 outbreak. Resour. Policy 2021, 74, 102265. [Google Scholar] [CrossRef] [PubMed]
  6. Diaz, E.M.; Molero, J.C.; Perez de Gracia, F. Oil price volatility and stock returns in the G7 economies. Energy Econ. 2016, 54, 417–430. [Google Scholar] [CrossRef] [Green Version]
  7. Luo, X.; Qin, S. Oil price uncertainty and Chinese stock returns: New evidence from the oil volatility index. Finance Res. Lett. 2017, 20, 29–34. [Google Scholar] [CrossRef]
  8. Kinateder, H.; Wagner, N. Oil and Stock Market Returns: Direction, Volatility, or Liquidity? SSRN Electron. J. 2017. [Google Scholar] [CrossRef]
  9. Batten, J.A.; Kinateder, H.; Szilagyi, P.G.; Wagner, N.F. Liquidity, surprise volume and return premia in the oil market. Energy Econ. 2019, 77, 93–104. [Google Scholar] [CrossRef]
  10. Szakmary, A.; Ors, E.; Kyoung Kim, J.; Davidson, W.N. The predictive power of implied volatility: Evidence from 35 futures markets. J. Bank. Finance 2003, 27, 2151–2175. [Google Scholar] [CrossRef]
  11. Becker, R.; Clements, A.E.; White, S.I. Does implied volatility provide any information beyond that captured in model-based volatility forecasts? J. Bank. Finance 2007, 31, 2535–2549. [Google Scholar] [CrossRef] [Green Version]
  12. Mencía, J.; Sentana, E. Valuation of VIX derivatives. J. Finance Econ. 2013, 108, 367–391. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, Z.; Tseng, H.K.; Wu, J.S.; Ding, Z. Implied volatility relationships between crude oil and the U.S. stock markets: Dynamic correlation and spillover effects. Resour. Policy 2020, 66, 101637. [Google Scholar] [CrossRef]
  14. Dutta, A.; Bouri, E.; Saeed, T. News-based equity market uncertainty and crude oil volatility. Energy 2021, 222, 119930. [Google Scholar] [CrossRef]
  15. Lanne, M.; Saikkonen, P. Modeling the U.S. Short-Term Interest Rate by Mixture Autoregressive Processes. SSRN Electron. J. 2005, 1, 96–125. [Google Scholar] [CrossRef] [Green Version]
  16. Ali, F.D.; Bin, L.; Leqin, W. Revisiting the risk-return relation in the South African stock market. Afr. J. Bus. Manag. 2012, 6, 11411–11415. [Google Scholar] [CrossRef] [Green Version]
  17. Chang, K.L. Does the return-state-varying relationship between risk and return matter in modeling the time series process of stock return? Int. Rev. Econ. Finance 2016, 42, 72–87. [Google Scholar] [CrossRef]
  18. Aslanidis, N.; Christiansen, C.; Savva, C.S. Risk-return trade-off for European stock markets. Int. Rev. Finance Anal. 2016, 46, 84–103. [Google Scholar] [CrossRef] [Green Version]
  19. Abakah, E.J.A.; Tiwari, A.K.; Alagidede, I.P.; Gil-Alana, L.A. Re-examination of risk-return dynamics in international equity markets and the role of policy uncertainty, geopolitical risk and VIX: Evidence using Markov-switching copulas. Finance Res. Lett. 2022, 47, 102535. [Google Scholar] [CrossRef]
  20. Lundblad, C. The risk return tradeoff in the long run: 1836–2003. J. Finance Econ. 2007, 85, 123–150. [Google Scholar] [CrossRef]
  21. Liu, B.; Ji, Q.; Fan, Y. Dynamic return-volatility dependence and risk measure of CoVaR in the oil market: A time-varying mixed copula model. Energy Econ. 2017, 68, 53–65. [Google Scholar] [CrossRef]
  22. Wu, S.-J.; Lee, W.-M. Intertemporal risk–return relationships in bull and bear markets. Int. Rev. Econ. Finance 2015, 38, 308–325. [Google Scholar] [CrossRef]
  23. Glosten, L.R.; Jagannathan, R.; Runkle, D.E.; Breen, W.; Hansen, L.; Hess, P.; Hsieh, D.; Judson, R.; Kocher-Lakota, N.; Mcdonald, R.; et al. On the Relation between the Expected Value and the Volatility of the Nominal Excess Return on Stocks. J. Finance 1993, 48, 1779–1801. [Google Scholar] [CrossRef]
  24. Bollerslev, T.; Osterrieder, D.; Sizova, N.; Tauchen, G. Risk and return: Long-run relations, fractional cointegration, and return predictability. J. Finance Econ. 2013, 108, 409–424. [Google Scholar] [CrossRef] [Green Version]
  25. Chen, M. Risk-return tradeoff in Chinese stock markets: Some recent evidence. Int. J. Emerg. Mark. 2015, 10, 448–473. [Google Scholar] [CrossRef]
  26. Basher, S.A.; Sadorsky, P. Hedging emerging market stock prices with oil, gold, VIX, and bonds: A comparison between DCC, ADCC and GO-GARCH. Energy Econ. 2016, 54, 235–247. [Google Scholar] [CrossRef] [Green Version]
  27. Ji, Q.; Fan, Y. Modelling the joint dynamics of oil prices and investor fear gauge. Res. Int. Bus. Finance 2016, 37, 242–251. [Google Scholar] [CrossRef]
  28. Haugom, E.; Langeland, H.; Molnár, P.; Westgaard, S. Forecasting volatility of the U.S. oil market. J. Bank. Finance 2014, 47, 1–14. [Google Scholar] [CrossRef] [Green Version]
  29. Lux, T.; Segnon, M.; Gupta, R. Forecasting crude oil price volatility and value-at-risk: Evidence from historical and recent data. Energy Econ. 2016, 56, 117–133. [Google Scholar] [CrossRef] [Green Version]
  30. Dutta, A. Modeling and forecasting oil price risk: The role of implied volatility index. J. Econ. Stud. 2017, 44, 1003–1016. [Google Scholar] [CrossRef]
  31. Lv, W. Does the OVX matter for volatility forecasting? Evidence from the crude oil market. Phys. A Stat. Mech. Its Appl. 2018, 492, 916–922. [Google Scholar] [CrossRef]
  32. Chen, H.; Liu, L.; Li, X. The predictive content of CBOE crude oil volatility index. Phys. A Stat. Mech. Its Appl. 2018, 492, 837–850. [Google Scholar] [CrossRef]
  33. Benedetto, F.; Mastroeni, L.; Quaresima, G.; Vellucci, P. Does OVX affect WTI and Brent oil spot variance? Evidence from an entropy analysis. Energy Econ. 2020, 89, 104815. [Google Scholar] [CrossRef]
  34. Maghyereh, A.I.; Awartani, B.; Bouri, E. The directional volatility connectedness between crude oil and equity markets: New evidence from implied volatility indexes. Energy Econ. 2016, 57, 78–93. [Google Scholar] [CrossRef] [Green Version]
  35. Bouri, E.; Jain, A.; Biswal, P.C.; Roubaud, D. Cointegration and nonlinear causality amongst gold, oil, and the Indian stock market: Evidence from implied volatility indices. Resour. Policy 2017, 52, 201–206. [Google Scholar] [CrossRef]
  36. Choi, S.Y.; Hong, C. Relationship between uncertainty in the oil and stock markets before and after the shale gas revolution: Evidence from the OVX, VIX, and VKOSPI volatility indices. PLoS ONE 2020, 15, e0232508. [Google Scholar] [CrossRef] [PubMed]
  37. Naeem, M.A.; Hasan, M.; Agyemang, A.; Chowdhury, M.I.H.; Balli, F. Time-frequency dynamics between fear connectedness of stocks and alternative assets. Int. J. Finance Econ. 2021, 1–14. [Google Scholar] [CrossRef]
  38. Aboura, S.; Chevallier, J. Leverage vs. feedback: Which effect drives the oil market? Finance Res. Lett. 2013, 10, 131–141. [Google Scholar] [CrossRef] [Green Version]
  39. Chen, Y.; He, K.; Yu, L. The Information Content of OVX for Crude Oil Returns Analysis and Risk Measurement: Evidence from the Kalman Filter Model. Ann. Data Sci. 2015, 2, 471–487. [Google Scholar] [CrossRef] [Green Version]
  40. Lin, J.B.; Tsai, W. The relations of oil price change with fear gauges in global political and economic environment. Energies 2019, 12, 2982. [Google Scholar] [CrossRef]
  41. Agbeyegbe, T.D. An inverted U-shaped crude oil price return-implied volatility relationship. Rev. Finance Econ. 2015, 27, 28–45. [Google Scholar] [CrossRef]
  42. Echaust, K.; Just, M. Tail dependence between crude oil volatility index and WTI oil price movements during the COVID-19 pandemic. Energies 2021, 14, 4147. [Google Scholar] [CrossRef]
  43. Li, S.; Li, J.; Lu, X.; Sun, Y. Exploring the dynamic nonlinear relationship between crude oil price and implied volatility indices: A new perspective from MMV-MFDFA. Phys. A Stat. Mech. Its Appl. 2022, 603, 127684. [Google Scholar] [CrossRef]
  44. Shaikh, I. The relation between implied volatility index and crude oil prices. Eng. Econ. 2019, 30, 556–566. [Google Scholar] [CrossRef] [Green Version]
  45. Silva, J.C.A., Jr. An S-Shaped Crude Oil Price Return-Implied Volatility Relation: Parametric and Nonparametric Estimations. Int. J. Econ. Finance 2017, 9, 54. [Google Scholar] [CrossRef] [Green Version]
  46. Fousekis, P. Crude oil price and implied volatility: Insights from non-parametric quantile regressions. Stud. Econ. Finance 2019, 36, 168–182. [Google Scholar] [CrossRef]
  47. Galvao, A.F. Unit root quantile autoregression testing using covariates. J. Econom. 2009, 152, 165–178. [Google Scholar] [CrossRef]
  48. Troster, V. Testing for Granger-causality in quantiles. Econom. Rev. 2018, 37, 850–866. [Google Scholar] [CrossRef]
  49. Han, H.; Linton, O.; Oka, T.; Whang, Y. The Cross-Quantilogram: Measuring Quantile Dependence and Testing Directional Predictability between Time Series. J. Econom. 2016, 193, 251–270. [Google Scholar] [CrossRef] [Green Version]
  50. Yang, M. Journal of Economic Dynamics & Control The risk return relationship: Evidence from index returns and realised variances. J. Econ. Dyn. Control 2019, 107, 103732. [Google Scholar] [CrossRef]
  51. Agnolucci, P. Volatility in crude oil futures: A comparison of the predictive ability of GARCH and implied volatility models. Energy Econ. 2009, 31, 316–321. [Google Scholar] [CrossRef]
  52. Herrera, A.M.; Karaki, M.B. The effects of oil price shocks on job reallocation. J. Econ. Dyn. Control. 2015, 61, 95–113. [Google Scholar] [CrossRef]
  53. Herrera, A.M.; Karaki, M.B.; Rangaraju, S.K. Where do jobs go when oil prices drop? Energy Econ. 2017, 64, 469–482. [Google Scholar] [CrossRef] [Green Version]
  54. Karaki, M.B. Nonlinearities in the response of real GDP to oil price shocks. Econ. Lett. 2017, 161, 146–148. [Google Scholar] [CrossRef]
  55. Bahel, E.; Marrouch, W.; Gaudet, G. The economics of oil, biofuel and food commodities. Resour. Energy Econ. 2013, 35, 599–617. [Google Scholar] [CrossRef] [Green Version]
  56. Ma, F.; Wahab, M.I.M.; Huang, D.; Xu, W. Forecasting the realized volatility of the oil futures market: A regime switching approach. Energy Econ. 2017, 67, 136–145. [Google Scholar] [CrossRef]
  57. Chen, Y.; Zou, Y. Examination on the Relationship between OVX and Crude Oil Price with Kalman Filter. Procedia Comput. Sci. 2015, 55, 1359–1365. [Google Scholar] [CrossRef] [Green Version]
  58. Raggad, B. Can implied volatility predict returns on oil market? Evidence from Cross-Quantilogram Approach. Resour. Policy 2023, 80, 103277. [Google Scholar] [CrossRef]
  59. Koenker, R.; Xiao, Z. Unit root quantile autoregression inference. J. Am. Stat. Assoc. 2004, 99, 775–787. [Google Scholar] [CrossRef] [Green Version]
  60. Troster, V.; Shahbaz, M.; Uddin, G.S. Renewable energy, oil prices, and economic activity: A Granger-causality in quantiles analysis. Energy Econ. 2018, 70, 440–452. [Google Scholar] [CrossRef]
  61. Ye, C.; Chen, Y.; Inglesi-Lotz, R.; Chang, T. CO2 emissions converge in China and G7 countries? Further evidence from Fourier quantile unit root test. Energy Environ. 2020, 31, 348–363. [Google Scholar] [CrossRef]
  62. Sakov, A.; Bickel, P.J. An Edgeworth expansion for the m out of n bootstrapped median. Stat. Probab. Lett. 2000, 49, 217–223. [Google Scholar] [CrossRef]
  63. Broock, W.A.; Scheinkman, J.A.; Dechert, W.D.; LeBaron, B. A test for independence based on the correlation dimension. Econom. Rev. 1996, 15, 197–235. [Google Scholar] [CrossRef]
  64. Bai, J.; Perron, P. Computation and analysis of multiple structural change models. J. Appl. Econom. 2003, 18, 1–22. [Google Scholar] [CrossRef] [Green Version]
  65. Dickey, D.A.; Fuller, W.A. Distribution of the Estimators for Autoregressive Time Series with a Unit Root. J. Am. Stat. Assoc. 1979, 74, 427–431. [Google Scholar] [CrossRef]
  66. Phillips, P.C.B.; Perron, P. Testing for a unit root in time series regression. Biometrika 1988, 75, 335–346. [Google Scholar] [CrossRef]
  67. Zivot, E.; Andrews, D.W.K. Further evidence on the Great Crash, the oil-price shock, and the and Unit-Root hypothesis. J. Bus. Econ. Stat. 1992, 10, 251–270. [Google Scholar]
  68. Johansen, S. Estimation and Hypothesis Testing of Cointegration Vectors in Gaussian Vector Autoregressive Models. Econometrica 1991, 59, 1551–1580. [Google Scholar] [CrossRef]
  69. Johansen, S. Likelihood-Based Inference in Cointegrated Vector Autoregressive Models; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  70. Xiao, Z. Quantile cointegrating regression. J. Econom. 2009, 150, 248–260. [Google Scholar] [CrossRef] [Green Version]
  71. Lin, J.B.; Liang, C.C.; Tsai, W. Nonlinear relationships between oil prices and implied volatilities: Providing more valuable information. Sustainability 2019, 11, 3906. [Google Scholar] [CrossRef]
Figure 1. Evolution of daily levels and returns of WTI crude oil and OVX.
Figure 1. Evolution of daily levels and returns of WTI crude oil and OVX.
Mathematics 11 00528 g001
Figure 2. Q-Q plots of the series.
Figure 2. Q-Q plots of the series.
Mathematics 11 00528 g002
Figure 3. Sample cross-quantilograms, to capture directional predictability from WTI to OVX, when both return series are in the same quantiles. Results with lag k = 1, 2, …, 30. The dashed lines denote the 95% bootstrapped confidence intervals.
Figure 3. Sample cross-quantilograms, to capture directional predictability from WTI to OVX, when both return series are in the same quantiles. Results with lag k = 1, 2, …, 30. The dashed lines denote the 95% bootstrapped confidence intervals.
Mathematics 11 00528 g003
Figure 4. Box–Ljung test statistic Q ^ α p for each lag p using ρ ^ α k , to capture directional predictability from WTI to OVX, when both return series are in the same quantiles. The dashed lines are the 95% bootstrap confidence intervals.
Figure 4. Box–Ljung test statistic Q ^ α p for each lag p using ρ ^ α k , to capture directional predictability from WTI to OVX, when both return series are in the same quantiles. The dashed lines are the 95% bootstrap confidence intervals.
Mathematics 11 00528 g004
Figure 5. Sample cross-quantilograms, to capture directional predictability from OVX to oil, when both return series are in the same quantiles. See notes to Figure 3.
Figure 5. Sample cross-quantilograms, to capture directional predictability from OVX to oil, when both return series are in the same quantiles. See notes to Figure 3.
Mathematics 11 00528 g005
Figure 6. Box–Ljung test statistic Q ^ α p for each lag p using ρ ^ α k , to capture directional predictability from OVX to WTI, when both return series are in the same quantiles. See notes to Figure 4.
Figure 6. Box–Ljung test statistic Q ^ α p for each lag p using ρ ^ α k , to capture directional predictability from OVX to WTI, when both return series are in the same quantiles. See notes to Figure 4.
Mathematics 11 00528 g006
Figure 7. Sample cross-quantilograms for WTI returns in the lower quantile, i.e., q 2 α 2 for α 2 = 0.1 , to capture directional predictability from WTI return to OVX. See notes to Figure 3.
Figure 7. Sample cross-quantilograms for WTI returns in the lower quantile, i.e., q 2 α 2 for α 2 = 0.1 , to capture directional predictability from WTI return to OVX. See notes to Figure 3.
Mathematics 11 00528 g007
Figure 8. Sample cross-quantilograms for WTI returns in the median quantile, i.e., q 2 α 2 for α 2 = 0.5 , to capture directional predictability from WTI to OVX returns. See notes to Figure 3.
Figure 8. Sample cross-quantilograms for WTI returns in the median quantile, i.e., q 2 α 2 for α 2 = 0.5 , to capture directional predictability from WTI to OVX returns. See notes to Figure 3.
Mathematics 11 00528 g008
Figure 9. Sample cross-quantilograms for WTI returns in the upper quantile, i.e., q 2 α 2 for α 2 = 0.9 and OVX in the other 19 quantiles, i.e., q 1 α 1 for α 1   ϵ   0.05 ,   0.10 , ,   0.95 . See notes to Figure 3.
Figure 9. Sample cross-quantilograms for WTI returns in the upper quantile, i.e., q 2 α 2 for α 2 = 0.9 and OVX in the other 19 quantiles, i.e., q 1 α 1 for α 1   ϵ   0.05 ,   0.10 , ,   0.95 . See notes to Figure 3.
Mathematics 11 00528 g009
Figure 10. Sample cross-quantilograms for OVX returns in the lower quantile, i.e., q 2 α 2 for α 2 = 0.1 , to capture directional predictability from OVX to WTI returns. See notes to Figure 3.
Figure 10. Sample cross-quantilograms for OVX returns in the lower quantile, i.e., q 2 α 2 for α 2 = 0.1 , to capture directional predictability from OVX to WTI returns. See notes to Figure 3.
Mathematics 11 00528 g010
Figure 11. Sample cross-quantilograms for OVX returns in the median quantile, i.e., q 2 α 2 for α 2 = 0.5 , to capture directional predictability from OVX to WTI returns. See notes to Figure 3.
Figure 11. Sample cross-quantilograms for OVX returns in the median quantile, i.e., q 2 α 2 for α 2 = 0.5 , to capture directional predictability from OVX to WTI returns. See notes to Figure 3.
Mathematics 11 00528 g011
Figure 12. Sample cross-quantilograms for OVX returns in the upper quantile, i.e., q 2 α 2 for α 2 = 0.9 , to capture directional predictability from OVX to WTI returns. See notes to Figure 3.
Figure 12. Sample cross-quantilograms for OVX returns in the upper quantile, i.e., q 2 α 2 for α 2 = 0.9 , to capture directional predictability from OVX to WTI returns. See notes to Figure 3.
Mathematics 11 00528 g012
Table 1. Summary statistics of returns series for WTI and OVX.
Table 1. Summary statistics of returns series for WTI and OVX.
WTIOVX
Mean−0.0360.037
Median0.053−0.27
Maximum37.47457.821
Minimum−39.087−43.991
Std. Dev.2.7474.938
Skewness-0.4021.057
Kurtosis35.97416.026
Jarque-Bera153,031.3 ***24,497.87 ***
Ljung–Box85.20 ***51.134 ***
Observations33763376
Correlation−0.325 ***
Notes: The sample period is 10 May 2007 to 17 April 2020. *** denotes statistical significance at the 1% level.
Table 2. Results of Brock–Dechert–Scheinkma (BDS) statistics from VAR residuals.
Table 2. Results of Brock–Dechert–Scheinkma (BDS) statistics from VAR residuals.
mWTIOVX
20.029 ***0.018 ***
30.056 ***0.033 ***
40.074 ***0.04 ***
50.083 ***0.045 ***
60.087 ***0.046 ***
20.029 ***0.018 ***
Note: The values denote the BDS test built on the residuals of WTI and OVX returns in a VAR. m represents the embedding dimension of the BDS test. *** signifies rejection of the null of residuals being i.i.d. at the 1% significance level.
Table 3. Results for conventional unit root tests for the two series.
Table 3. Results for conventional unit root tests for the two series.
OVXWTI
ADF Test−59.922 ***−62.137 ***
PP Test−60.180 ***−62.107 ***
ZA Test−36.084 ***−29.688 ***
Note: ***, **, and * indicate the rejection of the null hypothesis at the 1% 5% and 10% significance levels, respectively.
Table 4. Quantile autoregression unit root test results.
Table 4. Quantile autoregression unit root test results.
WTIOVX
τ β τ T-StatCV β τ T-StatCV
0.050.047−10.595−3.37−0.154−23.545−2.31
0.10.031−19.564−3.35−0.111−42.006−2.31
0.150.049−25.391−3.293−0.092−48.172−2.31
0.20.075−29.558−3.236−0.055−55.28−2.313
0.250.065−34.652−3.236−0.047−60.029−2.385
0.30.038−38.02−3.208−0.039−62.469−2.53
0.350.023−43.397−3.195−0.041−63.142−2.588
0.40.015−47.051−3.136−0.03−65.601−2.65
0.450.011−48.649−3.073−0.029−67.283−2.644
0.5−0.002−51.437−3.07−0.024−65.054−2.704
0.55−0.012−53.711−3.066−0.02−62.884−2.742
0.6−0.01−52.972−3.010.007−59.946−2.854
0.65−0.023−51.49−2.9640.016−54.166−2.895
0.7−0.044−49.292−2.8840.018−49.539−2.888
0.75−0.067−45.158−2.7640.038−45.41−2.902
0.8−0.092−39.279−2.5750.067−41.718−2.901
0.85−0.121−34.714−2.4910.052−38.062−2.914
0.9−0.168−28.373−2.310.088−22.758−3.036
0.95−0.149−13.036−2.310.075−7.644−3.154
Notes: This table displays the outcomes of the quantile unit root test for the grid of 19 quantiles, τ = 0.05, …, 0.95. The procedure consists of testing the null hypothesis of the unit root H0: β(τ) = 1. The null hypothesis of unit roots is rejected at quantile τ when the test-statistics are less than 5% the critical value (CV).
Table 5. Results of linear and quantile cointegration tests.
Table 5. Results of linear and quantile cointegration tests.
Panel A: Johansen Linear Cointegration Test
Trace statistic
H0: rank = 0 (15.41)
Max. eigenvalue statistic
H0: rank = 0 (14.07)
WTI-OVX20.202 (**)18.992 (**)
Panel B: Quantile Cointegration Test
modelcoefficient S u p τ V ^ n τ critical values
1%5%10%
WTI versus OVXbeta7406.825 (***)250.470141.558115.485
gamma935.6308 (***)20.81413.5429.543
Notes: Panel A shows the result for the test of Johansen [68,69] for oil price and OVX, taken in natural logarithm form. Numbers in parentheses next to H0: rank = 0 are the 5% critical values of the corresponding test statistic. Panel B shows the test statistic of the quantile cointegration model between the considered variables computed using an equally spaced grid of 19 quantiles [0.05; 0.95]. *** and ** denote the rejection of the null hypothesis at the 1% and 5% significance level, respectively.
Table 6. p-values of the standard test of Granger causality (in mean).
Table 6. p-values of the standard test of Granger causality (in mean).
Number of Lags
123
H 0 :  WTI does not Granger-cause OVX0.3410.350.441
H 0 :  OVX does not Granger-cause WTI0.00 **0.001 **0.002 **
Notes: ** denotes statistical significance at the 5% level.
Table 7. Granger causality in quantiles: subsample p-values.
Table 7. Granger causality in quantiles: subsample p-values.
WTI to OVXOVX to WTI
Number of LagsNumber of Lags
Quantiles123123
0.050.00030.00030.00030.00030.00030.0003
0.10.00030.00030.00030.00030.00030.0003
0.150.00030.00030.00030.00030.00030.0003
0.20.00030.00030.00030.00030.00030.0003
0.250.00030.00030.00030.00030.00030.0003
0.30.00030.00030.00030.00030.00030.0003
0.350.00030.00030.00030.01510.00310.0031
0.40.16240.14450.35520.34170.79670.7699
0.450.25970.31210.16730.00340.01110.0089
0.50.00030.00030.00030.00030.00030.0003
0.550.00030.00030.00030.00030.00030.0003
0.60.00030.00030.00030.00030.00030.0003
0.650.00030.00030.00030.00030.00030.0003
0.70.00030.00030.00030.00030.00030.0003
0.750.00030.00030.00030.00030.00030.0003
0.80.00030.00030.00030.00030.00030.0003
0.850.00030.00030.00030.00030.00030.0003
0.90.00030.00030.00030.00030.00030.0003
0.950.00860.00620.01230.935611
Note: The values in bold are the p-values indicating the rejection of the null hypothesis of no quantile Granger causality at the 1% significance level.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Raggad, B.; Bouri, E. Quantile Dependence between Crude Oil Returns and Implied Volatility: Evidence from Parametric and Nonparametric Tests. Mathematics 2023, 11, 528. https://doi.org/10.3390/math11030528

AMA Style

Raggad B, Bouri E. Quantile Dependence between Crude Oil Returns and Implied Volatility: Evidence from Parametric and Nonparametric Tests. Mathematics. 2023; 11(3):528. https://doi.org/10.3390/math11030528

Chicago/Turabian Style

Raggad, Bechir, and Elie Bouri. 2023. "Quantile Dependence between Crude Oil Returns and Implied Volatility: Evidence from Parametric and Nonparametric Tests" Mathematics 11, no. 3: 528. https://doi.org/10.3390/math11030528

APA Style

Raggad, B., & Bouri, E. (2023). Quantile Dependence between Crude Oil Returns and Implied Volatility: Evidence from Parametric and Nonparametric Tests. Mathematics, 11(3), 528. https://doi.org/10.3390/math11030528

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