. International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 87 International Journal of Energy Economics and Policy ISSN: 2146-4553 available at http: www.econjournals.com International Journal of Energy Economics and Policy, 2019, 9(3), 87-105. Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress Zouheir Ahmed Mighri*, Majid Ibrahim Alsaggaf Department of Finance and Insurance, College of Business, University of Jeddah, Saudi Arabia. *Email: zmighri@gmail.com Received: 10 December 2018 Accepted: 03 March 2019 DOI: https://doi.org/10.32479/ijeep.7643 ABSTRACT This paper attempts to estimate the relationship between oil prices and financial stress using weekly data for the period December 31, 1993 to July 15, 2016. The analysis is carried out using the cointegration framework. Both the linear and non-linear models for cointegration and related error correction models are estimated. The paper finds the threshold cointegration model more suitable than the linear cointegration models. It finds evidence of asymmetry in the adjustment process to equilibrium. It also finds that regimes with negative (below the threshold) changes of deviations adjust much faster than regimes with positive (above the threshold) changes of deviations, especially during a crisis period. Also, bi-direction causality is reported between the two variables. Keywords: Threshold Cointegration, Asymmetric Adjustment, Asymmetric Error Correction, Financial Stress, Oil Prices, Financial Crisis JEL Classifications: C22, C32, C58, G11 1. INTRODUCTION In the literature, it is well recognized that oil price shocks have detrimental effect on economic activity in developed and developing countries (Cunado and de Gracia, 2003; Cunado and de Gracia, 2005; Hamilton, 2011), especially for oil-importers. Although the effect of oil shocks on the macro-economy seems to have weakened through time, Kilian (2008) argued that this is partly due to increased demand for industrial output, which offsets the negative impact of an increase in oil price. Besides, Rafiq et al. (2009) argued that the negative impact of an increase in oil price is usually found to be higher than the positive impact of a fall in oil prices. In addition, Balke et al. (2002) claimed that monetary policy alone cannot account for this asymmetry. Transaction costs and financial stress are among the factors that lead to the asymmetric effect. The financial stress index literature is a rapidly developing one. Existing studies either focus on only constructing a financial stress index for a single country (Illing and Liu, 2006; Hakkio and Keeton, 2009; Morales and Estrada, 2010; Holló, 2012; Nazlioglu et al., 2015) or both on constructing an index for numerous countries and evaluating the link between financial stress and economic activity to examine how well financial stress index identifies known periods of financial distress (Slingenberg and de Haan, 2011; Cardarelli et al., 2011; Holló et al., 2012; Cevik et al., 2013; Mallick and Sousa, 2013; Chau and Deesomsak, 2014; Islami and Kurz-Kim, 2014). For one country, financial stress indices combine more indicators into one statistic than multi-country stress indices (Kliesen et al., 2012; Vermeulen et al., 2015). Most studies use market data, but some of them use both mixed market and balance sheet data (Holló et al., 2012), while others consider only balance sheet data (Morales and Estrada, 2010). Different ways are used by authors in order to combine indicators into an aggregate financial stress index. While most studies take the average of standardized variables, others use principal components analysis (Illing and Liu, 2006; Hakkio and Keeton, 2009). More recently, authors used portfolio theory based aggregation schemes This Journal is licensed under a Creative Commons Attribution 4.0 International License Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 201988 that take into account the correlation structure of stress indicators in order to quantify the level of systemic stress (Holló et al., 2012). The financial stress index is a relatively new concept, and the literature on impacts of oil prices on macroeconomic and financial variables is a lot broader than the financial stress index literature. The work developed by Chen et al. (2014) is probably the first study in the literature that examines the link between financial stress and oil prices. In their article, Chen et al. (2014) extended the Killian (2009) framework to identify an exogenous shock arising from changes in financial market conditions and examine the consequent macroeconomic impacts of oil price changes. Using the Kansas City financial stress index, global oil production, global real economic activity, and real oil prices, they found that financial stress index shocks trigger a significant negative response in real oil prices. Furthermore, a second study by Nazlioglu et al. (2015) examined whether there is a volatility transmission between oil prices and financial stress by means of the volatility spillover test. They used west texas intermediate (WTI) crude oil prices and Cleveland financial stress index (CFSI) for the period 1991-2014 and divide the sample into pre-crisis, in-crisis, and post-crisis periods due to the downward trend in oil price in 2008. According to these authors, the dynamic relationship between oil prices and financial stress can exist through two channels: Their impact on economic activity and on investor behavior. Their empirical results indicate that oil prices and financial stress index are dominated by long-run volatility. They also argued that a rise in oil prices depresses economic activity, may put pressure on credit markets, and negatively affect stock markets and the banking system. Besides, they claimed that, in times of high financial stress, economic activity slows down, leading to low energy demand and declining oil prices. In addition, crude oil markets are seen by investors as alternative investment areas to financial markets. With respect to oil price shocks, investors adjust their portfolios, which will have repercussions on financial asset prices. On the other hand, investors will be obliged to change their portfolios due to increased financial stress, which will have an effect on oil markets. Furthermore, Nazlioglu et al. (2015) claimed that financial stress influences economic activity through the bank lending channel via decreasing the amount of available credits and through financial leverage via changes in credit worthiness of borrowing businesses. Even though the relationship between financial stress and real economic activity or growth is well-studied (Demirguc-Kunt and Levine, 2001; Levine, 2005; Illing and Liu, 2006), the inter-temporal link between oil price and financial stress index is not yet well explored. This study examines whether there is asymmetric relationship between world oil prices and financial stress index. Considering the leading role of the U.S. financial system all over the world, the financial stress index for U.S. is taken as representative of the global financial stress. To the extent of our knowledge, this study is the first to explicitly examine asymmetric threshold cointegration between financial stress and world oil markets by employing the methodology developed by Enders and Siklos (2001). Most of studies adopt a linear cointegration framework (Engle and Granger, 1987), which assumes a linear long-run relationship among economic variables and a linear adjustment towards the equilibrium. However, the linear structure has been challenged as many economic variables display a nonlinear or asymmetric effect in their long-run relationship and short-term adjustment process (Granger and Lee, 1989; Enders and Granger, 1998). For cross-listings between crude oil prices and financial stress, the rationale of nonlinear modeling is more straightforward. Given the intricacies of the trading environment between these two markets, such as transaction costs, short-sell restrictions and exchange rate risks, arbitragers may only appear when price deviations from the equilibrium are large enough to cover their transaction costs and risk premia, implying an asymmetric adjustment process. There have been some extensions of linear cointegration models to capture this asymmetry. For instance, Balke and Fomby (1997) generalize the cointegration analysis to allow for a threshold effect in the adjustment process. Enders and Granger (1998) and Enders and Siklos (2001) expand the Engle-Granger two-step cointegration test by allowing for the possibility of asymmetric adjustment processes. In this study, we employ the Enders-Siklos threshold cointegration test to explore the long run asymmetric equilibrium relationship between oil prices and financial stress. The data set includes weekly observations from December 31, 1993 to July 15, 2016, and is divided into three sub-periods due to the downward trend in oil prices in 2008: The pre-oil crisis, the oil crisis, and the post oil crisis (pre-crisis, in-crisis, and post-crisis hereafter) periods. The remainder of the article is organized as follows. Section 2 describes the data. Section 3 outlines the econometric methodology. Section 4 presents the descriptive statistics and time series properties of data and discusses the empirical results. Section 5 is devoted to concluding remarks. 2. DATA SOURCES AND DESCRIPTION In this study, we use two variables, namely a measure of financial stress and oil prices at weekly frequency. The decision to carry out the analysis at weekly frequency is to better account for the dynamic relationships between oil and financial markets during the 2007-2009 global financial crisis. For world oil prices, we use the WTI spot crude oil prices, obtained from the FRED database of the St. Louis Federal Reserve Bank1. Given that our oil prices are weekly, we employ the St. Louis Fed’s Financial Stress Index (STLFSI)2 provided by the Federal Reserve Bank of St. Louis. While there are other financial stress index measures for the US, like the Chicago fed index, the Kansas city fed index (KCFSI), and the CFSI3, none of these indexes are available at weekly frequency. 1 https://fred.stlouisfed.org. 2 The St. Louis Fed’s Financial Stress Index (STLFSI) is based on 18 weekly data series. The actual index is constructed using a principal components analysis, which is a statistical method of extracting factors responsible for the comovements of the 18 variable groups. It is assumed that financial stress is the primary factor influencing this comovement, and by extracting this factor (the first principal component) financial stress index can be created. 3 Note that the CFSI, as a measure of stress in financial markets, has been unavailable since May 9, 2016, due to the discovery of errors that overestimated stress in the real estate and securitization markets. https://fred.stlouisfed.org Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 89 The STLFSI measures the degree of financial stress in the markets. This index is more comprehensive and overcomes the potential criticisms of focusing solely on one indicator. In combining several indicators, it has a broad coverage as it covers three important areas: (i) Interest rates (such as federal fund rate; 2 year, 10 year, and 30 year treasury; and corporate bond yield); (ii) yield curve (such as 10 year minus 3-month treasury; corporate bond minus 10 year treasury; 3 month TED spread); and (iii) other counterparty risk indicators (such as J.P. Morgan emerging markets bond index, Chicago board options exchange market volatility index, Merrill lynch bond market volatility index). Each of these variables captures some aspect of financial stress. The average value of the index, which begins in late 1993, is designed to be zero. Thus, zero is viewed as representing normal financial market conditions. Values below zero suggest below-average financial market stress, while values above zero suggest above-average financial market stress. Note that an increase in financial stress will be associated with higher funding costs and greater economic uncertainty, resulting in declining real economic activity. Moreover, an increased financial stress will render financial investors more risk averse, which will discourage investment in asset markets, resulting in falling asset prices, including oil prices (Hakkio and Keeton, 2009; Davig and Hakkio, 2010). The data set includes weekly observations from December 31, 1993 to July 15, 2016 and it is divided into three sub-periods: The pre-crisis period from December 31, 1993 to July 27, 2007, the crisis period from August 3, 2007 to March 27, 2009, and the post-crisis period from April 3, 2009 to July 15, 2016. Even though the WTI starts earlier, the starting date of the sample is constrained by availability of STLFSI, and all available data since the start of this study is included. The recent global financial crisis has some unique features, such as the length, breadth, and crisis sources. Numerous studies use major economic and financial events in order to determine the crisis length and source ad hoc (Baur, 2012; Dimitriou and Kenourgios, 2013; Dimitriou et al., 2013; Mighri and Mansouri, 2014). Besides, the choice of the sub-periods is based on the downward trend in oil prices within the crisis date (Mollick and Assefa, 2013; Turhan et al., 2013). These studies suggested that when oil prices are used, separate analyses are necessary before, at and after the crisis period. In this study, the length of the global financial crisis and its phases are specified following an economic approach. We define a relatively long crisis period based on all major international financial and economic news events representing the global financial crisis. We use the official timelines provided by Federal Reserve Board of St. Louis (2009) and the BIS (2009), among others, in order to choose the crisis period. According to these studies, the timeline of the global financial crisis is separated in four phases. Phase 1 described as “initial financial turmoil spans from 1 August 2007 to 15 September 2008. Phase 2 is defined as “sharp financial market deterioration” and spans from 16 September 2008 to 31 December 2008. Phase 3 described as “macroeconomic deterioration” spans from January 1, 2009 to March 31, 2009. Phase 4 described as a phase of “stabilization and tentative signs of recovery” (post-crisis period), including a financial market rally, spans from 1 April 2009 to the end of the sample period. For that reason, the crisis can be defined from August 2007 to March 2009 covering the first three phases. In the light of the literature, we therefore question the impact of the recent global financial crisis on the financial stress and oil price link and thus the data is divided into three sub-periods. 3. ECONOMETRIC METHODOLOGY Cointegration has been widely used to investigate relationship among price variables. The two major cointegration methods are Johansen and Engle-Granger two-step approaches. Both of them assume symmetric relationship between variables. In recent years, threshold cointegration has been increasingly used in price transmission studies. Balke and Fomby (1997) proposed a two-step approach for examining threshold cointegration on the basis of the approach developed by Engle and Granger (1987). Enders and Granger (1998) and Enders and Siklos (2001) further generalize the standard Dickey-Fuller test by allowing for the possibility of asymmetric movements in time-series data. This makes it possible to test for cointegration without maintaining the assumption of a symmetric adjustment to a long-term equilibrium. Thereafter, the method has been widely applied to analyze asymmetric price transmission. In this study, linear cointegration, threshold cointegration, and asymmetric error correction models are employed to examine the oil price and financial stress dynamics. These models will be able to assess asymmetric price dynamics in both the long term and short term. 3.1. Linear Cointegration Analysis In this study, the focus variables are weekly as well as monthly price series of crude oil and financial stress in the United States. As usual, their stochastic properties of non-stationarity and order of integration can be evaluated using the Augmented Dickey-Fuller (ADF) (Dickey and Fuller, 1979), and Phillips-Perron (Phillips and Perron, 1988) unit root tests. If both the price series appear to have a unit root, then it is appropriate to conduct cointegration analysis to assess their interaction. Econometric literature proposes different methodological alternatives to empirically analyze the long-run relationships and dynamics interactions between two or more time-series variables. Two cointegration methods widely used are the full information maximum likelihood-based Johansen approach and Engle-Granger two-step approach (Engle and Granger, 1987; Enders, 2004). The Johansen approach is a multivariate generalization of the Dickey-Fuller test (Johansen, 1988; Johansen and Juselius, 1990). It concentrates on the relationship between the rank of a matrix and its characteristic roots in a vector autoregression. The Johansen approach starts with a vector autoregressive model and then reformulates it into a vector error correction model as follows: Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 201990 Vt = π1Vt−1+...+πkVt−K+εt (1) ∆ Γ ∆ ΠV K V Vt i t i t K ti = − + + =∑ − −1 1  (2) Where Vt is a vector of the price at date (week) t for crude oil (yt) and financial stress (xt), K is the number of lags, and εt is the error term. The relationship among the coefficients for the two equations is given as: 1 p i jj I  = Γ = − +∑ (3) 1 K h h I  = Π = − +∑ (4) Where I is an identity matrix. Two types of tests, i.e., the trace and maximum eigenvalue statistics, can be used to detect the number of cointegrating vectors, r, among the variables in Vt. The Engle-Granger two-stage approach focuses on the time series property of the residuals from the long-term equilibrium relationship (Engle and Granger, 1987). The first step of the analysis consists in determining a break point into the relationship that defines the long run relationship between the crude oil prices and financial stress index: yt = ξ0+ξ1xt+εt (5) ∆ ∆ε ρε ϕ εt t i p i t i tz= + +− = −∑1 1 (6) Where yt and xt denote the oil prices and financial stress, respectively, ξ0, ξ1, ρ and φi are parameters to be estimated, εt is the disturbance term, which should be stationary if any long- run relationship exists between the two integrated price series, indicates the first difference, ˆt is the estimated residuals, ρ measures the speed of convergence of the system, zt is a white noise disturbance term, and p denotes the number of lags. The parameter ξ1 indicates the long-run elasticity of price transmission and gives the magnitude of adjustment of the crude oil price to variations of the financial stress index. If ξ1<1, changes in the financial stress index are not fully passed onto the crude oil price. In the first stage of estimating the long-term relationship among the variables yt and xt, the financial stress is chosen to be placed on the right side and assumed to be the driving force. In the second stage, the estimated residuals ˆ t are used to conduct a unit root test (Engle and Granger, 1987). The number of lags is chosen so there is no serial correlation in the regression residuals. It can be selected using the Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), or Ljung-Box Q test. If the null hypothesis of ρ=0 is rejected, then the residual series from the long-term equilibrium is stationary and the focal variables of yt and xt are cointegrated. Rejecting the null hypothesis of no cointegration ρ=0 in favor of the alternative hypothesis −2<ρ<0 implies that the {εt} sequence is stationary with mean zero. Any deviations from the long-run value of the disturbance term t are ultimately eliminated. Convergence is assured if −2<ρ<0. As such, Eq. (5) is an attractor such that εt can be written as an error correction model. The change in εt equals ρ multiplied by εt−1 regardless of whether εt−1≥0 or εt−1<0. 3.2. Cointegration Analysis with Structural Breaks Residual-based cointegration tests (Engle and Granger, 1987) assume that cointegrating vectors are constant over time. However, if there is a regime shift in the series, there will be a shift in the cointegrating vector as well. In such circumstances, these standard tests could lead to incorrect inferences about the long- run relationship of the price series. Furthermore, Phillips (1986) shows that if a structural break exists in the data, but is omitted from the cointegration relationship, this could lead to spurious rejections when the null of no cointegration is wrongly rejected. For the Engle and Granger (1987) test, such spurious rejections tend to occur for breaks that are located either too early in the sample or when the magnitude of the break increases. Thus, the power of the Engle and Granger (1987) test to find cointegration is severely affected by the presence of breaks in the level or the trend function in the cointegration relationship. Gregory and Hansen (1996a; 1996b) addressed this issue and proposed a residual-based cointegration test that allows for the possibility of regime shifts either in the intercept or the entire vector of coefficients. Gregory and Hansen (1996a; 1996b) analyzed four models and then tested the null hypothesis of no cointegration. Model 1 (see Eq.5) is the standard cointegration model where no changes in the intercept or a trend function are allowed under the null hypothesis. The other three models include shifts in either the intercept (Level shift model C) or trend (Level shift model with trend C/T) or shifts in the intercept and slope vector of coefficients (Regime shift model C/S). Model C/S is unique in the sense it allows the long-run equilibrium relationship to rotate as well as shift in parallel fashion. The break point in any model is determined endogenously within the data series. Level shift model C can be expressed as follows: y xt tt t t= + + +ξ ξ φ ξ0 0 10 ’ ’’  (7) In this parameterization, 0 ’ represents the intercept before the shift, and 0 ’’ denotes the change in the intercept at the time of the shift. Level shift model with trend C/T can be represented by: y t xt tt t t= + + + +ξ ξ φ β ξ υ0 0 10 ’ ’’ (8) Regime shift model C/S is given as: y x xt tt t t tt t= + + + +ξ ξ φ ξ ξ φ κ0 0 1 10 0 ’ ’’ ’ ’’ (9) In this case, 0 ’ and 0 ’’ are as in the level shift model C. 1 ’ denotes the cointegrating slope coefficients before the regime shift, and 1 ’’ denotes the change in the slope coefficients. A time trend into the regime shift model (C/S/T) could be also introduced: Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 91 y t x xt tt t t tt t= + + + + +ξ ξ φ β ξ ξ φ η0 0 1 10 0 ’ ’’ ’ ’’ (10) In these four models above, the structural break is modeled by the introduction of a dummy variable tt0 , which takes values (0,1) depending on the nature of the structural break. φ τt if t t if t t = > ≤    1 0 0 0 (11) Where t0 is the unknown parameter denoting the timing of the change point. In all four models postulated, the null hypothesis of no cointegration can be tested by examining whether the residuals of the ordinary least-squares (OLS) regression applied to Eqs. (5) and (7)-(10), respectively are stationary processes. The procedure for computing the test statistic for each possible regime shift t0 ∈T involves four steps. In essence, it involves the search for the smallest value of either the modified Phillips-Perron (PP) ( Z * and Zt * ) or ADF(ADF*) test statistic across all possible break points: Z Z t t T   * inf ( )= ∈0 0 (12) Z Z tt t T t * inf ( )= ∈0 0 (13) ADF ADF t t T * inf ( )= ∈0 0 (14) 3.3. Threshold Cointegration Analysis The implicit assumption of linear and symmetric adjustment (Engle and Granger, 1987) is problematic. Enders and Siklos (2001) proposed a two-regime threshold cointegration approach to entail asymmetric adjustment in cointegration analysis. They argued that the Engle-Granger cointegration test is likely to lead to misspecification errors when the adjustment of the error correction term is asymmetric. They remedy this error by expanding the Engle-Granger two-step cointegration test to incorporate an asymmetric error correction term. In the next step, we determine whether or not the disturbance term εt is stationary by considering an asymmetric test methodology in the form of Threshold Autoregressive (TAR) cointegration model as proposed by Enders and Granger (1998) and Enders and Siklos (2001). The alternative model modifies Eq. (6) such that: ∆εt=It ρ1(εt−1−τ)+(1−It) ρ2 (εt-1−τ)+μt (15) Where ρ1, ρ2 are coefficients, τ is the value of the threshold, μt is a white-noise disturbance and It is the Heaviside indicator such that. I if ift t t = ≥ <    − − 1 0 1 1 ε τ ε τ (16) In order for {εt} to be stationary, a necessary condition is −2<(ρ1,ρ2)<0. If the variance of μt is sufficiently large, it is also possible for one value of ρj to be in the range of −2 and 0 and for the other value to equal zero. Although there is no convergence in the regime with the unit-root (i.e., the regime in which ρj=0), large realizations of μt will switch the system into the convergent regime. In both cases, under the null assumption of no cointegration between the variables, the -statistic for the null hypothesis ρ1=ρ2=0 has a nonstandard distribution. Rejecting this assumption means that Eq. (15) is an attractor such that the equilibrium value of the {εt} is τ. The adjustment process is (ρ1εt−1−τ) if the lagged value of εt is above its long-run equilibrium value, while if the lagged value of εt is below its long-run equilibrium value, the adjustment is ρ2(εt−1−τ). If −1<|ρ1|<|ρ2|<0, negative discrepancies will be more persistent than positive discrepancies. Moreover, Tong (1983) showed that the OLS estimates of ρ1 and ρ2 have an asymptotic multivariate normal distribution if the sequence {εt} is stationary. Therefore, if the null assumption ρ1=ρ2=0 is rejected, it is possible to test for symmetric adjustment (i.e., ρ1=ρ2) using a standard F-test. Rejecting both the null assumptions ρ1=ρ2=0 and ρ1=ρ2 indicates the existence of threshold cointegration and asymmetric adjustment. Since the exact nature of the nonlinearity may not be known, Enders and Siklos (2001) consider another kind of asymmetric cointegration test methodology that allows the adjustment to be contingent on the change in εt−1 (i.e., ∆εt−1) instead of the level of εt−1. In this case, the Heaviside indicator of Eq. (16) becomes. I if ift t t = ≥ <    − − 1 0 1 1 ∆ ∆ ε τ ε τ (17) This specification is especially relevant when the adjustment is such that the series exhibits more “momentum” in one direction than in the other (Thompson, 2006; Kuo and Enders, 2004; Enders and Siklos, 2001; Enders and Granger, 1998). That is, the speed of adjustment depends on whether εt is increasing (i.e., widening) or decreasing (i.e., narrowing). According to Thompson (2006), among others, if |ρ1|<|ρ2|, then increase in εt tend to persist, whereas decreases revert back to the threshold quickly. The resulting model is called momentum-threshold autoregressive (M-TAR) cointegration model. The TAR model captures asymmetrically deep movements if, for instance, positive deviations are more prolonged than negative deviations. The M-TAR model allows the autoregressive decay to depend on ∆εt−1. As such, the M-TAR specification can capture asymmetrically “sharp” movements in {εt} sequence (Caner and Hansen, 2001). In both the TAR and M-TAR cointegration processes, the null assumption of ρ1=ρ2=0 could be tested, while the null hypothesis of symmetric adjustment may be tested by the restriction, ρ1=ρ2. Generally, there is no presumption to whether to use TAR or M-TAR specifications. Thus, it is recommended to select the adjustment mechanism by a model selection criterion such as AIC or BIC. Furthermore, if the errors in Eq. (15) are serially correlated, it is possible to use the augmented form of the test: ∆ ∆ε ρ ε τ ρ ε τ ϕ εt t t t t i p i t i tI I v= −( )+ −( ) −( )+ +− − = −∑1 1 2 1 11 (18) To use the tests, we first regress εt on a constant and call the residuals , { ˆ }t which are the estimates of (εt−1−τ). In a second step, we set the indicator according to Eq. (16) or Eq. (17) and estimate the following regression: Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 201992 ( ) ( ) ( )1 1 2 1 1ˆ ˆ ˆ1 ˆ p t t t t t i t i ti I I v        − − −=∆ = − + − − + ∆ +∑ (19) The number of lags p is specified to account for serially correlated residuals and it can be selected using AIC, BIC, or Ljung-Box Q test. In several applications, there is no reason to expect the threshold to correspond with the attractor (i.e., τ=0). In such circumstances, it is necessary to estimate the value of along with the values of ρ1 and ρ2. A consistent estimate of the threshold t can be obtained by adopting the methodology of Chan (1993). A super consistent estimate of the threshold value can be attained with several steps. First, the process involves sorting in ascending order the threshold variable, i.e., 1ˆt − for the TAR model or the 1 ˆ t −∆ for the M-TAR model. Second, the potential threshold values are determined. If the threshold value is to be meaningful, the threshold variable must actually cross the threshold value (Enders, 2004). Thus, the threshold value τ should lie between the maximum and minimum values of the threshold variable. In practice, the highest and lowest 15% of the values were removed from the search to ensure an adequate number of observations on each side. The middle 70% values of the sorted threshold variable are used as potential threshold values. Third, the TAR or M-TAR model is estimated with each potential threshold value. The sum of squared errors (SSE) for each trial can be calculated and the relationship between the SSE and the threshold value can be examined. Finally, the threshold value yielding the lowest SSE is deemed to be the consistent estimate of the threshold. Given these considerations, a total of four models are used in this study. They are TAR- Eq. (16) with τ=0; consistent TAR-Eq. (16) with τ estimated; MTAR- Eq. (17) with τ=0; and consistent MTAR- Eq. (17) with τ estimated. Since there is generally no presumption on which specification is used, it is recommended to choose the appropriate adjustment mechanism via model selection criteria of AIC and BIC (Enders and Siklos, 2001). A model with the lowest AIC and BIC will be used for further analysis. Insights into the asymmetric adjustments in the context of a long term cointegration relationship can be obtained with two tests. First, an F-test is used to examine the null assumption of no cointegration (H0:ρ1=ρ2=0) 4 against the alternative of cointegration with either TAR or M-TAR threshold adjustment. Let Φ and Φ* denote the F-statistics for testing the null assumption of ρ1=ρ2=0 under the TAR and the M-TAR specifications, respectively. The distributions of Φ and Φ* are determined by the form of the attractor. The second one is a standard F-test to assess the null assumption of symmetric adjustment in the long-term equilibrium (H0: ρ1=ρ2). Rejection of the null hypothesis indicates the existence of an asymmetric adjustment process. 4 The null hypothesis of non stationarity is rejected if the sample value of F-test statistic exceeds the Enders-Granger critical value. The critical values of the -statistics for the null hypothesis ρ1=ρ2=0 using the TAR and M-TAR specifications are reported in the first and second panels of Table 1 in Kuo and Enders (2004). 3.4. Asymmetric Error Correction Model with Threshold Cointegration According to Engle and Granger (1987), if all considered variables are cointegrated, then there will be a corresponding error correction model (ECM). The finding could be extended to threshold cointegration. This means that, if yt and xt are threshold cointegrated, then the ECM could be constructed as follows: ∆ ∆ ∆ x Z Z x y t x x t x t xj t j xj t j x t j p j p= + + + + + + − + − − − − − = =∑ ∑θ δ δ α β ν 1 1 1 1 , (20) And ∆ ∆ ∆ y Z Z x y t y y t y t j p yj t j yj t j y t j p= + + + + + + − + − − − = − − ∑ ∑ =θ δ δ α β ν 1 1 1 1 , (21) Where 1 1ˆt t tZ I  + − −= and 1 1ˆ(1 ) t t tZ I  − − −= − ; the parameters δ+ and δ− represent the adjustment speed of the coefficients of different sized deviations; θ is a constant; αj and βj are the coefficients of the lagged difference terms; p is the number of lags and νt is a white noise. The subscripts x and y are used in order to differentiate between the coefficients of variables xt and yt, respectively. t denotes time, and j represents lags. The equilibrium correction specification (ECM) of Engle and Granger (1987) assumes that the adjustment process due to disequilibrium among the variables is symmetric. In order to incorporate asymmetries, two extensions on the ECM model have been made. Error correction terms and first differences on the variables are decomposed into positive and negative values, as proposed by Granger and Lee (1989). The second extension adds the threshold cointegration mechanism to the Granger and Lee (1989) approach. The resulting asymmetric error correction model with threshold cointegration has the following form: ∆ ∆ ∆ x Z Z x x t x x t x t xj t j xj t j j j j p p= + + + + + + − + − − − + − + − − − = =∑ ∑θ δ δ α α 1 1 1 1 == + − + − − − =∑ ∑+ +1 1 p p xj t j xj t j x tj y yβ β ν∆ ∆ , (22) And ∆ ∆ ∆ y Z Z x x t y y t y t yj t j yj t j j j j p p= + + + + + + − + − − − + − + − − − = =∑ ∑θ δ δ α α 1 1 1 1 == =∑ ∑ + − + − − −+ + 1 1 p p jyj t j yj t j y t y yβ β ν∆ ∆ , (23) The Heaviside indicator function is constructed from Eq. (16) or Eq. (18). The superscripts “+” and “−” indicate that the variables are split into positive and negative components. The first differences are defined as follows: Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 93 ∆ ∆x x x x x x x xt j t j t j t j t j t j t j t j − + − − − − − − − − − − = { } = − ≥ <  max , , , 0 0 1 1 1   ∆ ∆x x x x x x x xt j t j t j t j t j t j t j t j − − − − − − − − − − − − = { } = − < ≥  min , , , 0 0 1 1 1   ∆ ∆y y y y y y y yt j t j t j t j t j t j t j t j − + − − − − − − − − − − = { } = − ≥ <  max , , , 0 0 1 1 1   ∆ ∆y y y y y y y yt j t j t j t j t j t j t j t j − − − − − − − − − − − − = { } = − < ≥  min , , , 0 0 1 1 1   The lag p is specified to account serially correlated residuals and is selected using AIC statistic and Ljung-Box Q test. The above specifications are able to distinguish between long-run and short- run adjustments of xt and yt. The long-run adjustment is determined by the parameters δ+ and δ−, whereas, the short-run adjustment is governed by the parameters α α βj j j + − +, , and  j − for j=1,…, p. If  x x + −≠ and  y y + −≠ , then both xt and yt exhibit asymmetry in long-run adjustment. If either  xj xj + −≠ or. or both, xt displays asymmetry in short-run adjustment. Besides, if either  yj yj + ≠ − or  yj yj + −≠ or both, yt displays asymmetry in short-run adjustment. In this paper, four types of single or joint null hypotheses and F-tests are examined (Meyer and Von Cramon-Taubadel, 2004; Frey and Manera, 2007; Sun, 2011; Chen and Zhu, 2015; Mighri and Mansouri, 2016). The first type is the Granger causality test to examine the lead-lag relationship between xt and yt. The null hypothesis that xt does not lead yt can be tested by restricting H yj yj01 0:  + −= = for all lags j simultaneously and then employing an F-test. Similarly, the null hypothesis that yt does not lead xt can be tested by restricting H xj xj02 0:   + −= = for all lags j simultaneously and then employing an F-test. In our empirical analysis, we expect to see one of the following: If one variable Granger-causes the other, then the former variable leads the latter; if there is no causal relationship between the two variables, then there is no obvious connection between the two variables; or, if the two variables mutually Granger-cause each other, then the two variables are closely linked to each other. The second type of hypothesis is concerned with the distributed lag asymmetric effect o n i t s o w n v a r i a b l e ; t h a t i s , H xj xj03 0:  + −= = a n d H yj yj04 0:   + −= = . The third type of the null hypothesis is the cumulative symmetric effect which can be expressed as H j p xj j p xj05 1 1 : = + = −∑ ∑=  f o r x t a n d H j p yj j p yj06 1 1 : = + = −∑ ∑=  , for yt. Finally, the equilibrium adjustment path asymmetry can be examined with the null hypotheses of H07 :   + −= for each equation estimated (i.e.,  x x + −= for xt and  y y + −= for yt) to examine whether it is possible to get back to equilibrium after a shock, and if it is the case, how long it will take. 4. EMPIRICAL RESULTS 4.1. Descriptive Statistics and Unit Root Test Figure 1 displays the time series plots for the oil prices and St. Louis Fed Financial Stress Index (STLFSI). Three observations can be made (i) Oil prices and STLFSI have an evident comovement in general, which reveals a high possibility of cointegration between these two series. (ii) Although oil prices and STLFSI move together most of the time during our sample period, they also display divergent movement indicating possible nonlinear cointegration. (iii) The two series tend to move more closely during and after the crisis relative to the pre-crisis period. It seems that the link between STLFSI and world oil prices have changed through time, which motivated us to concentrate on the sub-sample analysis. Table 1 reports summary statistics of oil prices and STLFSI for different samples in order to examine to what extent the descriptive statistics of the oil prices and the STLFSI differ across these sub-periods. The highest mean and standard deviation are observed for oil prices during all the sub-periods. As a simple measurement for volatility, the standard deviations of oil prices are higher in the crisis period compared to those of the pre- crisis and post-crisis periods. Skewness is a simple measure of asymmetry and kurtosis is a measurement for peaked or flatted distribution relative to a Gaussian distribution. We observe that oil price has negative skewness and is left tailed in both the crisis and post-crisis periods; although it has positive skewness and hence right tailed in the pre-crisis period. This stylized fact is Figure 1: Dynamics of weekly oil prices and St. Louis Fed Financial Stress Index Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 201994 also supported by the Jarque-Bera test statistic which rejects the null assumption of normality. The STLFSI has its positive and the greatest mean in the in-crisis period, which means that there is a significant or higher stress. It has positive mean in the pre-crisis and negative mean in the post- crisis. As illustrated in Figure 1, these features provide that the pre-crisis period is characterized as a normal stress period while moderate stress levels are observed in the post-crisis period. The STLFSI has smallest standard deviation during the crisis-period compared to the pre-and post-crisis periods. The skewness and kurtosis measures indicate deviation from normality. Skewness shows that STLFSI is right tailed in the crisis and post-crisis periods; it has left tail in the pre-crisis period. Kurtosis indicates that financial stress is less peaked during the pre-crisis and in-crisis periods compared to the post-crisis period. Moreover, the Jarque- Bera statistic shows non-normal behavior of financial stress. The different data characteristics apparent in the summary statistics for different considered periods lead to the question of whether the correlations between oil prices and STLFSI vary across these sub- periods as well. At first glance, correlation appears stronger in the pre-crisis and in-crisis periods than that in the post-crisis period. However, the positive correlation in the in-crisis period turns back to negative following the global financial crisis. This means that as financial stress goes from significant level to moderate level, oil prices go up. This reversed relationship becomes more apparent after 2008 according to Figure 1. The non-stationarity properties of the oil prices and financial stress indices are investigated by applying the ADF (Dickey and Fuller, 1979), PP (Phillips and Perron, 1988), and Kwiatkowski-Phillips- Schmidt-Shin (KPSS) (Kwiatkowski et al., 1992) unit root tests. For both ADF and PP tests, we consider three different cases: The case without constant and trend, the case with constant, as well as the case with both constant and trend. The results for the ADF and PP unit root tests are summarized in both Tables 2 and 3. These tests indicate that oil prices are characterized by a unit root process, implying that the shocks are permanent and not corrected over time. For the STLFSI, even though both the unit root tests support evidence on stationary process in the post-crisis period, this mean reverting process is not supported in the full-sample, pre-crisis and in-crisis periods. In addition, the KPSS test also reveals that both WTI and STLFSI series fail to reject the alternative assumption of a unit root at the 1% level of significance or better. Conversely, these series accept the null hypothesis of a stationary at 1% level of significance or better when tested for a unit root in first differences. Therefore, we conclude that both oil prices and STLFSI are integrated processes of order one, or unit root processes. 4.2. Results of the Linear Cointegration Analysis To conduct the linear cointegration analyses, we use both the Johansen and Engle-Granger approach. The application of the Johansen approach requires the determination of a lag length for the model, which is based on the lowest AIC and BIC. Without prior information, three model specifications with trend, constant, or no intercept are entailed (Table 4). For instance, with only a trend included, the Johansen maximum eigenvalue statistic (λmax) is 21.152 for the null hypothesis of no cointegrating vector between the prices of WTI and STLFSI. These are significant at the 5% level, indicating that the null hypothesis is rejected. However, for the null hypothesis of one cointegrating vector, the λmax tatistic decreases to 5.472, which is not statistically significant at all. Therefore, the maximum eigenvalue statistic conclude that there is one cointegrating vector. Similarly, the Johansen trace statistic (λtrace) also supports the conclusion that WTI and STLFSI are cointegrated. The Engle-Granger cointegration test is implemented through two steps. In the first step, the long-term relationship between the price series is estimated, as specified in Eq. (17). The estimates for the coefficients on the financial stress index (i.e., ξ1) are highly statistically significant in all cases. In the second step, the residual series is used to conduct a unit root test with the specification in Eq. (18). As reported in Table 5, the sufficiently optimal lag, used for addressing the problem of serial correlation, is chosen based on the AIC and Ljung-Box Q statistics. Besides, the statistics from the unit root test (i.e., ρ) are statistically significant. The results of the Engle-Granger cointegration test provide evidence for the alternative hypothesis of linear cointegration in all cases except for the post crisis period. 4.3. Results of the Cointegration Analysis in the Presence of Structural Breaks Using the Gregory and Hansen (1996a; 1996b) tests, we tested for a bivariate cointegration relationship between oil prices and Table 1: Descriptive statistics Variable Whole sample Pre-crisis In-crisis Post-crisis WTI STLFSI WTI STLFSI WTI STLFSI WTI STLFSI Mean 51.139 0.000 31.351 0.215 87.666 1.798 79.623 −0.811 Median 42.520 0.091 26.340 0.418 91.180 1.099 85.660 −0.934 Maximum 142.520 5.701 75.630 1.470 142.520 5.701 112.300 2.640 Minimum 11.000 −1.657 11.000 −0.906 32.980 −0.077 28.140 −1.657 SD 31.306 1.000 16.670 0.570 29.969 1.483 21.507 0.712 Skewness 0.557 1.299 1.121 −0.479 −0.178 1.039 −0.671 2.124 Kurtosis 2.040 7.514 3.056 1.980 2.059 2.802 2.265 9.015 Jarque−bera 106.004*** (0.000) 1330.457*** (0.000) 148.664*** (0.000) 57.889*** (0.000) 3.670 (0.160) 15.800*** (0.000) 37.155*** (0.000) 860.726*** (0.000) Observations 1177 1177 709 709 87 87 381 381 Correlation −0.417 − −0.758 − −0.619 − −0.309 − ***, **, and *indicate respectively statistical significance at the 1%, 5%, and 10% levels. WTI: West texas intermediate, STLFSI: St. Louis fed’s financial stress index Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 95 financial stress in the US markets. Results of the residual-based tests for cointegration in models with regime shifts are reported in Table 6. In almost the cases, the hypothesis of cointegration with a structural break is not supported at better than 1% significance level. Taking into account the structural break in the cointegrating relationship between oil prices and financial stress, the adjustment of the oil price to changes in the financial stress in the long-run is explored allowing for the possibility of asymmetries in the error correction process. To deal with the issue of statistical inference in a cointegrated system with structural breaks, both TAR and MTAR models (Enders and Siklos, 2001) are used. 4.4. Results of the Threshold Cointegration Analysis From the linear cointegration tests, the price transmission mechanism between oil prices and financial stress index may be asymmetric. To investigate this possibility, it is necessary to go further than the usual concept of cointegration in order to allow for asymmetric cointegration and thus asymmetric price Table 3: Unit root tests for STLFSI Test Whole sample Pre-crisis Crisis Post-crisis Levels 1st diff. Levels 1st diff. Levels 1st diff. Levels 1st diff. ADF unit root test Exogenousa C CLT C CLT CLT C CLT NONE Lag lengthc 7 6 1 0 1 0 1 0 T-Stat. −2.377 −13.275 −1.827 −22.308 −2.266 −6.318 −5.272 −14.758 Test crit. valuese: 1% level −3.436 −3.966 −3.439 −3.971 −4.070 −4.070 −3.982 −2.571 Test crit. valuese: 5% level −2.864 −3.414 −2.865 −3.416 −3.464 −3.464 −3.422 −1.942 PP unit root test Exogenousb CLT CLT C CLT CLT CLT CLT NONE Bandwidthd 6 3 2 4 5 4 5 8 Adj. T-Stat. −3.115 −25.206 −1.681 −22.229 −2.095 −6.397 −5.414 −15.226 Asymptotic crit. valuesee: 1% level −3.966 −3.966 −3.439 −3.971 −4.068 −4.070 −3.982 −2.571 Asymptotic crit. valuesee: 5% level −3.414 −3.414 −2.865 −3.416 −3.463 −3.464 −3.422 −1.942 KPSS unit root test Exogenousb CLT CLT CLT CLT C C CLT CLT Bandwidthf 26 5 21 0 7 4 15 8 LM Stat 1.3792 0.0268 0.6042 0.0620 0.8146 0.0679 0.3727 0.1688 Asymptotic crit. valuesee: 1% level 0.7390 0.2160 0.2160 0.2160 0.7390 0.7390 0.2160 0.2160 Asymptotic crit. valuesee: 5% level 0.4630 0.1460 0.1460 0.1460 0.4630 0.4630 0.1460 0.1460 CLT: Constant, Linear Trend. C indicates Constant, aModel selection is based on Schwarz information criterion, bModel selection is based on Newey-West Bandwith, cLag length selection is based on Schwarz information criterion, maxlag=22, dThe bandwidth selection is defined by using Bartlett kernel, eMacKinnon (1996), fModel selection is based on Newey-West bandwith using Bartlett kernel. STLFSI: St. Louis Fed Financial Stress Index, ADF: Augmented dickey fuller, PP: Phillips-perron, KPSS: Kwiatkowski−Phillips−Schmidt−Shin Table 2: Unit root tests for oil prices (WTI) Test Whole sample Pre-crisis Crisis Post-crisis Levels 1st diff. Levels 1st diff. Levels 1st diff. Levels 1st diff. ADF unit root test Exogenousa CLT CLT C CLT C C C C Lag lengthc 1 0 3 2 0 0 1 0 T−Stat. −1.792 −28.098 0.061 −14.208 −0.359 −7.757 −1.422 −15.028 Test crit. valuese: 1% level −3.966 −3.966 −3.439 −3.971 −3.508 −3.509 −3.447 −3.447 Test crit. valuese: 5% level −3.414 −3.414 −2.865 −3.416 −2.896 −2.896 −2.869 −2.869 PP unit root test Exogenousb CLT CLT C CLT C C C C Bandwidthd 14 12 13 17 4 4 6 2 Adj. t-Stat. −2.296 −29.186 0.443 −22.696 −0.660 −7.939 −1.344 −15.011 Asymptotic crit. valuese: 1% level −3.966 −3.966 −3.439 −3.971 −3.508 −3.509 −3.447 −3.447 Asymptotic crit. valuese: 5% level −3.414 −3.414 −2.865 −3.416 −2.896 −2.896 −2.869 −2.869 KPSS unit root test Exogenousb CLT CLT CLT CLT CLT CLT CLT CLT Bandwidthf 26 14 22 14 7 4 15 4 LM Stat. 0.3054 0.0527 0.5704 0.0405 0.2633 0.1408 0.5033 0.0440 Asymptotic crit. valuese: 1% level 0.2160 0.2160 0.2160 0.2160 0.2160 0.2160 0.2160 0.2160 Asymptotic crit. valuese: 5% level 0.1460 0.1460 0.1460 0.1460 0.1460 0.1460 0.1460 0.1460 CLT: Constant, Linear Trend, C indicates Constant, aModel selection is based on Schwarz information criterion, bModel selection is based on Newey–West bandwith, cLag length selection is based on Schwarz information criterion, maxlag=22, dThe bandwidth selection is defined by using Bartlett kernel, eMacKinnon (1996), fModel selection is based on Newey–West bandwith using bartlett kernel. WTI: West Texas intermediate, ADF: Augmented dickey fuller, PP: Phillips-perron, KPSS: Kwiatkowski−Phillips−Schmidt−Shin Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 201996 transmission. We conduct a nonlinear cointegration analysis by using the threshold auto-regression models. A total of four models are considered in this study. They are TAR with τ=0, consistent TAR (C.TAR) with τ estimated, M-TAR with τ=0 and consistent M-TAR (C.M-TAR) with τ estimated. To address possible serial correlation in the residual series, we select an appropriate lag by specifying a maximum lag of 16. We use AIC, BIC and Ljung-Box Q statistics for diagnostic analyses on the residuals. In most cases, the value of the threshold τ is unknown and has to be estimated along the values of ρ1 and ρ1. Table 4: Results of the Johansen cointegration tests Test Specification Lag Statistic Critical value (%) 10 5 1 Johansen λmax r=1 Trend 5 5.472 10.49 12.25 16.26 r=0 Trend 5 21.152** 16.85 18.96 23.65 r=1 Constant 5 4.199 7.52 9.24 12.97 r=0 Constant 5 9.511 13.75 15.67 20.2 r=1 None 5 4.098 6.5 8.18 11.65 r=0 None 5 9.511 12.91 14.9 19.19 Johansen λtrace r≤1 Trend 5 5.472 10.49 12.25 16.26 r=0 Trend 5 26.624** 22.76 25.32 30.45 r≤1 Constant 5 4.199 7.52 9.24 12.97 r=0 Constant 5 13.71 17.85 19.96 24.6 r≤1 None 5 4.098 6.5 8.18 11.65 r=0 None 5 13.609 15.66 17.95 23.52 r is the number of cointegrating vectors. The critical values are from Enders (2004). ***, ** and *denote statistical significance at the 1%, 5% and 10% levels, respectively Table 5: Results of the linear (Engle-Granger) cointegration tests Variable Whole sample Pre-crisis In-crisis Post-crisis First step ξ0 51.139*** [61.64] 36.118*** [82.690] 110.160*** [27.518] 72.039*** [45.268] ξ1 −13.067*** [−15.74] −22.169*** [−30.900] −12.511*** [−7.266] −9.355*** [−6.338] Second step Lags (ρ) 9 3 1 1 ρ −0.005* [−1.975] −0.017** [−2.435] −0.047* [−1.763] −0.006 [−0.939] φ1 0.253*** [8.617] 0.185*** [4.887] 0.235** [2.228] 0.262*** [5.245] φ2 −0.104*** [−3.439] −0.070* [−1.837] – – φ3 0.147*** [4.820] 0.089** [2.377] – – φ4 −0.098** [−3.195] – – – φ5 0.043 [1.402] – – – φ6 −0.041 [0.186] – – – φ7 −0.012 [−0.411] – – – φ8 0.009 [0.329] – – – φ9 0.058* [1.980] – – – AIC 5332.006 2966.987 540.394 1723.438 BIC 5387.784 2989.806 547.792 1735.267 QLB (4) 0.998 0.793 0.124 0.87 QLB (8) 1.000 0.171 0.429 0.815 QLB (12) 0.127 0.115 0.263 0.518 The numbers in the brackets are t-values. ***, **and *denote statistical significance at the 1%, 5% and 10% levels, respectively Table 6: Results of Gregory-Hansen cointegration test Period Level shift model C Level shift model with trend C/T t-Stat. t0 Lag Date t-Stat. t0 Lag Date Whole sample −3.883 565 9 2004:10:22 −3.853 565 9 2004:10:22 Pre-crisis −4.108 568 3 2004:11:12 −4.741 568 3 2004:11:12 In-crisis −3.133 60 0 2008:09:19 −5.043 71 3 2008:12:05 Post-crisis −4.739 293 1 2014:11:07 −4.817 293 1 2014:11:07 Period Regime shift model C/S Regime shift model with trend C/S/T t-Stat. t0 Lag Date t-Stat t0 Lag Date Whole sample −5.170 795 12 2009:03:20 −3.888 565 9 2004:10:22 Pre-crisis −4.917 546 1 2004:06:11 −4.48 565 1 2004:10:22 In-crisis −5.728 71 3 2008:12:05 −3.288 60 0 2008:09:19 Post-crisis −5.069 293 1 2014:11:07 −4.822 293 1 2014:11:07 t-Stat. indicate smallest t-statistics using Gregory–Hansen cointegration test among possible break points. Three asterisks *** (resp. **, *) denote rejection of the null hypothesis at the 1% (resp. 5%, 10%) significance level. t0 denotes the break point corresponding to the smallest t-statistic. Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 97 We follow the Chan’s (1993) method to estimate the threshold values for consistent TAR and M-TAR models. The empirical results of the threshold cointegration tests with an unknown threshold value using TAR, M-TAR and their consistent counterparts are reported in Tables 7 and 8. The value is the optimal threshold for the indicator function. Under these conditions, we can reject the null hypothesis of threshold cointegration (ρ1=ρ2=0) for all considered cases. This means that there exists a cointegrating relationship between oil prices and financial stress index. Since WTI is cointegrated with STLFSI using the consistent TAR or M-TAR models, we examine whether their adjustment coefficients are different across positive and negative errors. This procedure is achieved by verifying the existence of an asymmetric cointegration, i.e., testing the null assumption of ρ1=ρ2. Notice that the asymmetry test only makes sense when the two previous tests reject the null hypothesis. That is, if the ρi coefficients estimated for the threshold are significantly different from zero, then the regression is nontrivial and testing for symmetry makes all the sense. Based on the “Principle of Parsimony,” AIC, BIC and Ljung-Box Q statistics suggest that the most applicable model for variables’ adjustment to long-run equilibrium is the M-TAR model with consistent threshold value for the whole sample, while the consistent TAR model is the best one for the in-crisis and post- crisis periods. It turns out that different lag specifications in the models have little impact of the final threshold values selected. The variation of the SSE by threshold value for consistent M-TAR model with a lag of twelve is presented in Figure 2. The lowest SSE for the consistent M-TAR model is 6851.566 at the threshold value of −1.093. Similarly, the best threshold value with the lowest SSE is estimated to be 24.604 for the consistent TAR model. Finally, while the four nonlinear threshold cointegration models have similar results (Table 7), the consistent M-TAR model has the lowest AIC statistic of 5301.544 and BIC statistic of 5377.438, and therefore, is deemed to be the best. As shown in Tables 7 and 8, we found limited evidence of asymmetric price transmission between oil prices and financial stress. Therefore, oil prices became cointegrated with the STLFSI, the adjustment mechanism is asymmetric and the speed of adjustment to the equilibrium is different when the last equilibrium error has different signs. This means that the change in the equilibrium error has a different impact on the adjustment speed to the new equilibrium. Focusing on the results from the consistent M-TAR model, the F-test for the null hypothesis of no cointegration has a statistic of 6.165 and it is highly significant at the 1% level. Thus, the oil prices and financial stress index are cointegrated with threshold adjustment. Furthermore, the F- statistic for the null hypothesis of symmetric price transmission has a value of 8.746 and it is also significant at the 1% level. Therefore, the adjustment process is asymmetric when WTI and STLFSI adjust to achieve the long- term equilibrium. When considering the in-crisis period and focusing on the results from the consistent TAR model, the F-test for the null assumption of no cointegration has a statistic of 3.235 and it is statistically significant at the 5% level. Thus, WTI and STLFSI are cointegrated with threshold adjustment. In addition, the F statistic for the null assumption of symmetric price transmission has a value of 3.277 and it is statistically significant at the 10% level. Therefore, during the global financial crisis period, the adjustment process is asymmetric when WTI and STLFSI adjust to achieve the long- term equilibrium. Focusing on the findings from the consistent TAR model for the post-crisis period, the F-test for the null hypothesis of no cointegration has a statistic of 2.32 and it is significant at the 10% level. Therefore, the oil prices and financial stress index are cointegrated with threshold adjustment during the post-crisis period. Furthermore, the F statistic for the null hypothesis of symmetric price transmission has a value of 3.753 and it is also significant at the 10% level. Thus, during the post-crisis period, the adjustment process is asymmetric when WTI and STLFSI adjust to achieve the long-term equilibrium. 4.5. Results of the Error Correction Model Given the nonlinear threshold cointegration results, the final step in our analysis is to proceed with the asymmetric error correction model in order to investigate the movement of the oil price and STLFSI index series in a long-run equilibrium relationship. For the whole sample period, the results of our estimations of the consistent M-TAR error correction models are illustrated in Table 9. Diagnostic analyses on the residuals with AIC, BIC and Ljung-Box Q statistics select a lag of eight for the model. The consistent M-TAR model is the best from the threshold cointegration analyses and the error correction terms are constructed using Eq. (17) and Eq. (19). Results show that WTI is cointegrated with STLFSI index and it also exhibits asymmetric adjustments. Besides, the short-term equilibrium adjustment process mainly occurs with STLFSI-index since δ+ = δ−. Moreover, there are three situations to reduce the price deviations between the two variables if they are cointegrated (Chen et al., 2013). Given the case STLFSI-index price is larger than WTI price, there are three situations to reduce the price deviations: (i) STLFSI-index price goes down and WTI price goes up; (ii) STLFSI-index price goes down and WTI price goes down as well, but STLFSI-index price drops more; (iii) STLFSI-index price goes up and WTI price goes up, but STLFSI-index price increases less. In our empirical results, for regimes with positive shocks (STLFSI- index price is higher than WTI price), the adjustment coefficient for STLFSI-index is 0.0001 and -0.003 for WTI, which means that, in the next period, WTI price will go up and STLFSI-index price will go down, and thus, the price deviation will decrease. For regimes with negative shocks (STLFSI-index price is lower than WTI price), the adjustment coefficient for STLFSI-index is −0.0002 and −0.018 for WTI, which means that, in the next period, WTI price will go down and STLFSI-index price will go down as well, Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 201998 Ta bl e 7: R es ul ts o f t he n on lin ea r (t hr es ho ld ) c oi nt eg ra ti on te st s V ar ia bl e W ho le s am pl e P re -c ri si s TA R C .T A R M TA R C .M TA R TA R C .T A R M TA R C .M TA R L ag s (ρ ) 12 12 12 12 7 7 7 7 T hr es ho ld (τ ) 0. 00 0 24 .6 04 0. 00 0 −1 .0 93 0. 00 0 −8 .2 37 0. 00 0 1. 22 ρ 1 −0 .0 05 [− 1. 63 7] −0 .0 08 ** * [− 2. 64 5] 0. 00 1 [0 .1 17 ] 0. 00 1 [0 .1 26 ] −0 .0 09 [− 0. 95 3] −0 .0 08 [− 0. 94 8] −0 .0 14 [− 1. 47 5] −0 .0 06 [− 0. 40 5] ρ 2 −0 .0 04 [− 0. 96 5] 0. 00 1 [0 .2 17 ] −0 .0 1* ** [− 2. 75 8] −0 .0 17 ** * [− 3. 51 0] −0 .0 26 ** [− 2. 32 3] −0 .0 3* * [− 2. 50 4] −0 .0 18 * [− 1. 66 5] −0 .0 19 ** [− 2. 34 0] φ 1 0. 25 8* ** [8 .7 80 ] 0. 26 0* ** [8 .8 65 ] 0. 23 9* ** [7 .7 64 ] 0. 22 6* ** [7 .2 60 ] 0. 17 2* ** [4 .5 25 ] 0. 17 2* ** [4 .5 31 ] 0. 17 3* ** [4 .5 19 ] 0. 16 8* ** [4 .3 49 ] φ 2 −0 .0 94 ** [− 3. 11 9] −0 .0 92 ** [− 3. 05 5] −0 .0 98 ** [− 3. 24 0] −0 .0 99 ** [− 3. 27 7] −0 .0 64 * [− 1. 67 9] −0 .0 64 * [− 1. 66 8] −0 .0 64 * [− 1. 66 2] −0 .0 64 * [− 1. 66 8] φ 3 0. 13 9* ** [4 .6 02 ] 0. 14 1* ** [4 .6 45 ] 0. 13 9* ** [4 .5 87 ] 0. 14 0* ** [4 .6 50 ] 0. 08 9* * [2 .3 23 ] 0. 08 9* * [2 .3 21 ] 0. 08 9* * [2 .3 26 ] 0. 08 9* * [2 .3 37 ] φ 4 −0 .1 02 ** * [− 3. 34 0] −0 .1 00 ** [− 3. 28 7] −0 .1 00 ** [− 3. 28 4] −0 .0 99 ** [− 3. 27 4] 0. 04 8 [1 .2 49 ] 0. 04 9 [1 .2 59 ] 0. 04 8 [1 .2 45 ] 0. 04 6 [1 .1 98 ] φ 5 0. 04 4 [1 .4 29 ] 0. 04 5 [1 .4 68 ] 0. 04 9 [1 .5 89 ] 0. 05 3* [1 .7 20 ] 0. 03 2 [0 .8 25 ] 0. 03 3 [0 .8 54 ] 0. 03 3 [0 .8 49 ] 0. 03 2 [0 .8 43 ] φ 6 −0 .0 41 [− 1. 32 0] −0 .0 39 [− 1. 25 8] −0 .0 39 [− 1. 29 0] −0 .0 38 [− 1. 24 8] −0 .0 75 * [− 1. 94 9] −0 .0 74 * [− 1. 93 6] −0 .0 74 * [− 1. 93 7] −0 .0 76 ** [− 1. 98 4] φ 7 −0 .0 13 [− 0. 41 1] −0 .0 11 [− 0. 37 3] −0 .0 12 [− 0. 38 3] −0 .0 14 [− 0. 45 1] −0 .0 70 * [− 1. 83 4] −0 .0 68 * [− 1. 79 0] −0 .0 69 * [− 1. 82 6] −0 .0 69 * [− 1. 82 3] φ 8 0. 02 1 [0 .6 69 ] 0. 02 1 [0 .7 01 ] 0. 01 6 [0 .5 25 ] 0. 01 8 [0 .5 82 ] – – – – φ 9 0. 05 6* [1 .8 32 ] 0. 05 8* [1 .8 87 ] 0. 05 6* [1 .8 42 ] 0. 05 6* [1 .8 42 ] – – – – φ 1 0 −0 .0 43 [− 1. 42 8] −0 .0 41 [− 1. 35 6] −0 .0 45 [− 1. 49 3] −0 .0 44 [− 1. 45 6] – – – – φ 1 1 −0 .0 83 ** [− 2. 73 2] −0 .0 81 ** [− 2. 67 4] −0 .0 79 ** [− 2. 62 3] −0 .0 80 ** [− 2. 67 1] – – – – φ 1 2 0. 08 9* * [3 .0 37 ] 0. 09 1* * [3 .1 12 ] 0. 08 8* * [2 .9 95 ] 0. 08 9* * [3 .0 61 ] – – – – To ta l o bs . 11 77 11 77 11 77 11 77 70 9 70 9 70 9 70 9 C oi nt o bs . 11 64 11 64 11 64 11 64 70 1 70 1 70 1 70 1 A IC 53 10 .3 32 53 06 .8 51 53 06 .2 80 53 01 .5 44 29 50 .7 20 29 49 .8 70 29 52 .0 54 29 51 .2 76 B IC 53 86 .2 26 53 82 .7 45 53 82 .1 74 53 77 .4 38 29 96 .2 46 29 95 .3 95 29 97 .5 79 29 96 .8 01 Q LB (4 ) 0. 99 9 0. 99 9 0. 99 9 0. 99 8 1. 00 0 1. 00 0 1. 00 0 1. 00 0 Q LB (8 ) 1. 00 0 1. 00 0 1. 00 0 1. 00 0 1. 00 0 1. 00 0 1. 00 0 1. 00 0 Q LB (1 2) 1. 00 0 1. 00 0 1. 00 0 1. 00 0 0. 88 0 0. 87 7 0. 88 0 0. 84 2 N o C I: Φ (H 0 : ρ 1= ρ 2 =0 ) 1. 79 4 [0 .1 67 ] 3. 52 1* * [0 .0 30 ] 3. 80 5* * [0 .0 23 ] 6. 16 5* ** [0 .0 02 ] 3. 06 1* * [0 .0 47 ] 3. 48 5* * [0 .0 31 ] 2. 39 7* [0 .0 92 ] 2. 78 4* [0 .0 62 ] N o A PT : F (H 0 : ρ 1= ρ 2 ) 0. 03 0 [0 .8 61 ] 3. 47 5* [0 .0 63 ] 4. 04 1* * [0 .0 45 ] 8. 74 6* ** [0 .0 03 ] 1. 38 2 [0 .2 40 ] 2. 22 4 [0 .1 36 ] 0. 06 4 [0 .8 00 ] 0. 83 3 [0 .3 62 ] ** *, * * an d *d en ot e st at is tic al s ig ni fic an ce a t t he 1 % , 5 % a nd 1 0% le ve ls , r es pe ct iv el y. N um be rs in p ar en th es es a re P v al ue s. N um be rs in b ra ck et s ar e t- va lu es Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 99 but WTI drops more and thus the price deviation will decrease. The adjusted R-squared value is 0.243 for the STLFSI-index and 0.131 for WTI. Moreover, the AIC and BIC statistics for WTI are both larger than those for the STLFSI-index. This means that the model specification is better fitted on the WTI price. Using the estimation results of the asymmetric ECM with nonlinear threshold cointegration, we also conduct the hypothesis testing described in Section 3 (paragraph 3.4). The hypotheses of Granger causality between the series are assessed with F-tests. The F-statistic of 5.724 reveals that STLFSI does Granger cause WTI. Besides, the F-statistic of 4.941 indicates that WTI does Granger cause STLFSI. This indicates that, in the short-term, both variables affect each other. Similarly, the F-statistic of 16.899 for STLFSI discloses that the lagged index series have significant impacts on its own index. Furthermore, the F-statistic of 5.134, for WTI, reveals that the lagged price series have significant impacts on its own price. Thus, in the short term, STLFSI and WTI have been evolving more dependently. A number of hypotheses are examined for asymmetric price transmission. The first one is the distributed lag asymmetric effect. In each price equation, the equality of the corresponding positive and negative coefficients for each of the eight lags is tested; in total, there are sixteen F-tests for this hypothesis. It turns out that an important number of them is statistically significant and distributed lag asymmetric effect does exist. Furthermore, the cumulative asymmetric effects are also examined. The largest F-statistic is 9.963 but only one of the four statistics is statistically significant at the 1% level. Thus, cumulative effects are asymmetric. The fourth examined asymmetry is the momentum equilibrium adjustment path asymmetries (H07:δ +=δ−). For STLFSI, the F-statistic is 1.706 with a P = 0.192. The point estimates of the coefficients for the error correction terms are 0.0001 with a t-value of 0.668 for positive error correction term and −0.0002 with a t-value of –0.82 for the negative one. In contrast, for WTI price, the F-statistic is 7.952 with a P = 0.005. Thus, there is momentum equilibrium adjustment asymmetry. The point estimates are −0.003 with a t-value of −0.746 for positive deviations and −0.018 with a t-value of −3.405 for negative deviations. The magnitude suggests that in the short term the WTI responds to the positive deviations by 0.3% in a week but by 1.8% to negative deviations. Measured in response time, positive and negative deviations take, respectively, 333.333 and 55.556 weeks to be fully digested. Therefore, in the short-term, WTI has a much faster reaction to negative deviations from long-term equilibrium than positive deviations. Further findings for the in-crisis period are reported in Table 9. Diagnostic analyses on the residuals with AIC, BIC and Ljung- Box Q statistics select a lag of four for the model. The consistent TAR model is the best from the threshold cointegration analyses and the error correction terms are constructed using Eq. (16) and Eq. (19). During the global financial crisis period, results show that WTI is cointegrated with STLFSI index and it also exhibits asymmetric adjustments. Besides, for regimes with positive shocks, the adjustment coefficient for STLFSI-index is 0.004 and −0.036 for WTI, which means that, in the next period, WTI price will go up and STLFSI-index price will go down, and thus, the T ab le 8 : R es ul ts o f t he n on lin ea r (t hr es ho ld ) c oi nt eg ra ti on te st s (c on ti nu ed ) V ar ia bl e In -c ri si s P os t- cr is is TA R C .T A R M TA R C .M TA R TA R C .T A R M TA R C .M TA R L ag s (ρ ) 4 1 1 1 1 1 1 1 T hr es ho ld (τ ) 0. 00 0 −2 7. 00 0 0. 00 0 –6 .1 5 0. 00 0 17 .0 19 0. 00 0 –1 .3 95 ρ 1 –0 .0 43 [– 1. 16 1] –0 .0 19 [– 0. 62 7] –0 .0 13 [– 0. 35 1] –0 .0 29 [– 0. 99 9] –0 .0 2* [– 1. 96 1] –0 .0 27 ** [– 2. 15 3] –0 .0 06 [– 0. 67 4] –0 .0 1 [– 1. 39 8] ρ 2 –0 .0 39 [– 0. 91 5] –0 .1 27 ** [– 2. 46 9] –0 .0 87 ** [– 2. 23 8] –0 .1 37 ** [– 2. 11 6] 0. 00 1 [0 .2 05 ] 0. 00 04 [0 .0 61 ] –0 .0 06 [– 0. 65 4] 0. 00 6 [0 .5 29 ] φ 1 0. 30 9* * [2 .7 94 ] 0. 24 5* * [2 .3 50 ] 0. 22 8* * [2 .1 69 ] 0. 22 7* * [2 .1 65 ] 0. 26 4* ** [5 .3 07 ] 0. 26 5* ** [5 .3 14 ] 0. 26 2* ** [5 .2 33 ] 0. 26 2* ** [5 .2 58 ] φ 2 –0 .1 29 –1 .1 46 – – – – – – – φ 3 0. 27 2* *2 .4 38 – – – – – – – φ 4 –0 .1 80 [– 1. 61 5] – – – – – – – To ta l o bs . 87 87 87 87 38 1 38 1 38 1 38 1 C oi nt o bs . 82 85 85 85 37 9 37 9 37 9 37 9 A IC 52 3. 59 7 53 9. 06 3 54 0. 36 3 54 0. 02 1 17 22 .4 19 17 21 .6 74 17 25 .4 38 17 24 .0 74 B IC 54 0. 44 4 54 8. 83 4 55 0. 13 3 54 9. 79 2 17 38 .1 69 17 37 .4 24 17 41 .1 88 17 39 .8 24 Q LB (4 ) 0. 99 6 0. 19 1 0. 22 4 0. 37 8 0. 90 6 0. 93 3 0. 87 0. 83 1 Q LB (8 ) 0. 99 3 0. 55 1 0. 60 1 0. 78 5 0. 81 6 0. 84 5 0. 81 5 0. 82 6 Q LB (1 2) 0. 94 2 0. 41 2 0. 32 9 0. 57 4 0. 53 0. 58 8 0. 51 8 0. 52 5 N o C I: Φ (H 0 : ρ 1= ρ 2 =0 ) 1. 05 8 [0 .3 52 ] 3. 23 5* * [0 .0 44 ] 2. 56 4* [0 .0 83 ] 2. 73 9* [0 .0 71 ] 1. 94 6 [0 .1 44 2] 2. 32 * [0 .0 99 ] 0. 43 9 [0 .6 45 ] 1. 11 9 [0 .3 28 ] N o A PT : F (H 0 : ρ 1= ρ 2 ) 0. 00 4 [0 .9 50 ] 3. 27 7* [0 .0 74 ] 1. 98 3 [0 .1 63 ] 2. 32 1 [0 .1 31 ] 3. 00 7* [0 .0 84 ] 3. 75 3* [0 .0 53 ] 0. 00 02 [0 .9 95 ] 1. 35 6 [0 .2 45 ] ** *, * * an d *d en ot e st at is tic al s ig ni fic an ce a t t he 1 % , 5 % a nd 1 0% le ve ls , r es pe ct iv el y. N um be rs in p ar en th es es a re P v al ue s. N um be rs in b ra ck et s ar e t- va lu es Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019100 price deviation will decrease. For regimes with negative shocks, the adjustment coefficient for STLFSI-index is −0.002 and −0.076 for WTI, which means that, in the next period, WTI price will go down and STLFSI-index price will go down as well, but WTI drops more and thus the price deviation will decrease. The adjusted R-squared value is 0.337 for the STLFSI-index and 0.208 for WTI price. Likewise, the AIC and BIC statistics for WTI are both larger than those for the STLFSI-index, which means that the model specification is better fitted on the WTI price during the in-crisis period. Furthermore, the F-statistic of 1.811 reveals that STLFSI does Granger cause WTI. Moreover, the F-statistic of 3.098 indicates that WTI does Granger cause STLFSI. This indicates that, in the short-term, both variables affect each other during the crisis period. In the same way, the F-statistic of 2.34 for STLFSI discloses that the lagged index series have significant impacts on its own index. Additionally, the F-statistic of 0.781, for WTI, reveals that the lagged price series have insignificant impacts on its own price. Accordingly, in the short term, STLFSI and WTI have been evolving more dependently during the in-crisis period. The results for hypothesis description are also reported in Table 9. In each price equation, the equality of the corresponding positive and negative coefficients for each of the four lags is tested; in total, there are eight F-tests for this assumption. It turns out that a few numbers of them is statistically significant and the distributed lag asymmetric effect does exist. Furthermore, the cumulative asymmetric effects are not statistically significant and thus, cumulative effects are symmetric. In addition, the momentum equilibrium adjustment path asymmetry is examined. For STLFSI, the F-statistic is 1.893 with a P = 0.174. The point estimates of the coefficients for the error correction terms are 0.004 with a t-value of 1.889 for positive error correction term and −0.002 with a t-value of 0.502 for the negative one. On the contrary, for WTI price, the F-statistic is 0.377 with a P = 0.174. Thus, there is no momentum equilibrium adjustment asymmetry. The estimation and diagnostic results for the post-crisis period are reported in Table 9. Diagnostic analyses on the residuals with AIC, BIC and Ljung-Box Q statistics select a lag of five for the model. The consistent TAR model is the best from the threshold cointegration analyses and the error correction terms are constructed using Eq. (16) and Eq. (19). During the post- crisis period, results show that WTI is cointegrated with STLFSI index and it also exhibits asymmetric adjustments. Besides, for regimes with positive shocks, the adjustment coefficient for STLFSI-index is −0.0003 and −0.031 for WTI, which means that, in the next period, WTI price will go up and STLFSI-index price will go down, and thus, the price deviation will decrease. For regimes with negative shocks, the adjustment coefficient for STLFSI-index is −0.0001 and 0.002 for WTI, which means that, in the next period, WTI price will go down and STLFSI-index price will go down as well, but WTI drops more and thus the price deviation will decrease. The AIC and BIC statistics for WTI are both larger than those for the STLFSI-index, which means that the model specification is better fitted on the WTI price during the post-crisis period. Furthermore, the F-statistic of 0.953 reveals that STLFSI does not Granger cause WTI. Moreover, the F-statistic of 0.736 indicates that WTI does not Granger cause STLFSI. This indicates that, in the short-term, both variables do not affect each other during the post-crisis period. Also, the F-statistic of 2.938 for STLFSI discloses that the lagged index series have significant impacts on its own index. Additionally, the F-statistic of 3.869, for WTI, reveals that the lagged price series have significant impacts on its own price. Accordingly, in the short term, STLFSI and WTI have been evolving more independently during the post-crisis period. Finally, the results for hypothesis description are summarized in Table 9. First, they indicate the existence of some distributed lag asymmetric effects. Second, the cumulative effects are asymmetric. Third, there is a momentum equilibrium adjustment path asymmetry. For STLFSI, the F-statistic is 0.189 with a P = 0.664. The point estimates of the coefficients for the error correction terms are −0.0003 with a t-value of −0.619 for positive error correction term and −0.0001 with a t-value of −0.266 for the negative one. In contrast, for WTI price, the F-statistic is 3.274 with a P = 0.071. Furthermore, the point estimates are −0.031 with a t-value of −1.947 for positive deviations and 0.002 with a t-value of 0.216 for negative deviations. The magnitude suggests that in the short term the WTI responds to the positive deviations by 3.1% in a week but by 0.2% to negative deviations. Measured in response time, positive and negative deviations take, respectively, 32.26 and Figure 2: Threshold value for M-TAR (whole sample) Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 101 Ta bl e 9: R es ul ts o f t he a sy m m et ri c er ro r co rr ec ti on m od el s w it h th re sh ol d co in te gr at io n V ar ia bl e W ho le s am pl e (C .M − TA R ; l ag =8 ) In –c ri si s (C . T A R ; l ag =4 ) P os t– cr is is (C . T A R ; l ag =5 ) ST L F SI W T I ST L F SI W T I ST L F SI W T I C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c θ 0. 00 1 0. 14 4 −0 .1 12 −0 .6 84 0. 06 5 0. 71 3 2. 83 7* 1. 91 8 0. 01 5 1. 31 8 −0 .6 5* −1 .6 75 1  + 0. 43 3* ** 9. 41 2. 64 9* ** 2. 75 8 0. 56 5* ** 3. 09 8 4. 71 5 1. 60 4 −0 .0 33 −0 .3 41 5. 56 5* 1. 67 5 2  + 0. 05 8 1. 19 8 −6 .4 75 ** * −6 .3 55 0. 09 5 0. 47 −7 .1 25 ** −2 .2 0. 04 2 0. 43 1 0. 30 1 0. 08 9 3  + 0. 07 6 1. 52 8 1. 18 8 1. 13 8 0. 19 2 0. 95 3 2. 82 4 0. 87 −0 .0 89 −0 .9 05 1. 86 5 0. 54 6 4  + −0 .1 61 ** * −3 .1 91 −2 .9 01 ** * −2 .7 51 −0 .3 45 * −1 .7 44 −3 .2 59 −1 .0 23 0. 05 1 0. 52 3 −2 .0 82 −0 .6 09 5  + −0 .3 09 ** * −6 .0 15 3. 67 3* ** 3. 42 – – – – 0. 05 5 0. 56 3 −0 .2 16 −0 .0 64 6  + 0. 27 2* ** 5. 23 4 −0 .7 5 −0 .6 91 – – – – – – – – 7  + −0 .2 82 ** * −5 .4 1 0. 98 8 0. 90 8 – – – – – – – – 8  + 0. 00 4 0. 07 −3 .8 08 ** * −3 .5 66 – – – – – – – – 1  − 0. 13 5* * 2. 21 1 −1 .6 12 −1 .2 65 −0 .2 21 −0 .8 73 0. 46 3 0. 11 4 0. 41 9* ** 4. 11 7 1. 27 4 0. 35 9 2  − −0 .0 93 −1 .5 53 −0 .5 2 −0 .4 15 0. 01 7 0. 06 9 −3 .4 15 −0 .8 64 0. 01 6 0. 15 1 0. 25 0. 06 8 3  − 0. 10 5* 1. 74 3 −0 .9 12 −0 .7 27 0. 26 1 1. 17 0. 97 2 0. 27 0. 14 1 1. 35 2 −7 .3 28 ** −2 .0 14 4  − −0 .0 18 −0 .3 11 2. 02 8* 1. 65 2 −0 .3 04 −1 .3 42 6. 58 6* 1. 80 4 0. 08 9 0. 86 1 −1 .3 83 −0 .3 85 5  − 0. 16 2* ** 2. 78 −2 .4 23 ** −1 .9 93 – – – – 0. 01 1 0. 11 1 3. 44 1 0. 98 6 6  − −0 .0 12 −0 .2 13 −2 .4 11 ** −2 .0 05 – – – – – – – – 7  − 0. 06 5 1. 17 1 0. 21 5 0. 18 5 – – – – – – – – 8  − 0. 16 3* ** 2. 95 7 −0 .5 17 −0 .4 48 – – – – – – – – 1 + 0. 00 8* ** 2. 69 7 0. 10 7* 1. 84 2 0. 02 5. 1. 59 −0 .2 4 −0 .9 6 0. 00 1 0. 37 5 0. 19 5* 1. 93 4 2  + −0 .0 03 −1 .1 19 −0 .0 05 −0 .0 95 −0 .0 05 −0 .3 −0 .3 07 −1 .1 99 −0 .0 06 ** −2 .0 11 0. 17 * 1. 68 9 3  + −0 .0 05 −1 .6 28 0. 11 2* 1. 94 5 −0 .0 4* * −2 .4 13 0. 16 9 0. 63 6 0. 00 2 0. 53 2 0. 05 7 0. 56 6 4  + −0 .0 04 −1 .2 76 0. 04 3 0. 75 2 −0 .0 28 * −1 .7 01 −0 .0 07 −0 .0 28 0. 00 2 0. 59 3 0. 11 9 1. 18 4 0. 00 03 0. 12 0. 03 3 0. 57 4 – – – – −0 .0 01 −0 .3 45 0. 05 6 0. 56 5 6  + −0 .0 04 −1 .5 49 −0 .0 96 * −1 .6 91 – – – – – – – – 7  + 0. 00 3 0. 99 2 0. 04 8 0. 84 5 – – – – – – – – (C on td ... ) Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019102 (C on td ... ) V ar ia bl e W ho le s am pl e (C .M − TA R ; l ag =8 ) In –c ri si s (C . T A R ; l ag =4 ) P os t– cr is is (C . T A R ; l ag =5 ) ST L F SI W T I ST L F SI W T I ST L F SI W T I C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c 8  + 0. 00 7* ** 2. 62 2 0. 20 1* ** 3. 57 8 – – – – – – – – 1 − −0 .0 07 ** −2 .5 39 0. 27 ** * 4. 69 2 −0 .0 15 −0 .9 18 0. 51 9* 1. 90 7 −0 .0 04 . −1 .4 97 0. 39 7* ** 3. 92 4 2  − 0. 00 5* 1. 66 4 −0 .1 5* ** −2 .6 27 −0 .0 04 −0 .2 53 0. 15 8 0. 58 7 0. 00 2 0. 61 3 −0 .1 86 * −1 .7 92 3  − 0. 01 2* ** 4. 49 4 0. 01 1 0. 19 1 0. 05 6* ** 3. 58 5 −0 .0 65 −0 .2 56 0. 00 1 0. 37 8 −0 .0 62 −0 .6 03 4  − −0 .0 07 ** * −2 .5 79 −0 .0 54 −0 .9 47 −0 .0 06 −0 .3 07 0. 13 5 0. 46 7 −0 .0 01 −0 .4 75 −0 .1 39 −1 .3 52 −0 .0 01 −0 .2 49 0. 13 9* * 2. 40 9 – – – – 0. 00 1 0. 26 3 0. 08 9 0. 87 3 6  − 0. 00 1 0. 26 2 0. 06 3 1. 10 3 – – – – – – – – 7  − −0 .0 03 −0 .9 43 −0 .0 47 −0 .8 16 – – – – – – – – 8  − −0 .0 14 ** * −5 .1 07 0. 04 4 0. 77 5 – – – – – – – – δ+ 0. 00 01 0. 66 8 –0 .0 03 −0 .7 46 0. 00 4* 1. 88 9 −0 .0 36 −1 .0 95 −0 .0 00 3 −0 .6 19 −0 .0 31 * −1 .9 47 δ− −0 .0 00 2 −0 .8 2 −0 .0 18 ** * −3 .4 05 −0 .0 02 −0 .5 02 −0 .0 76 −1 .3 6 −0 .0 00 1 −0 .2 66 0. 00 2 0. 21 6 D ia gn os tic s R –s qu ar ed 0. 26 5 – 0. 15 6 – 0. 48 4 – 0. 38 4 – 0. 11 8 – 0. 12 5 – A dj us te d R –s qu ar ed 0. 24 3 – 0. 13 1 – 0. 33 7 – 0. 20 8 – 0. 06 3 – 0. 07 0 – F– st at . 11 .9 94 ** * (0 .0 00 ) 6. 15 9* ** (0 .0 00 ) 3. 28 2* ** (0 .0 00 2) 2. 18 5* * (0 .0 11 9) 2. 15 0* ** (0 .0 02 2) 2. 28 3* ** (0 .0 01 ) St at . D W 1. 98 7 (0 .8 26 ) 1. 98 0 (0 .6 74 ) 2. 02 3 (0 .9 18 ) 2. 01 0 (0 .8 34 ) 1. 99 6 (0 .8 56 ) 2. 00 1 (0 .8 78 ) A IC −2 01 8. 22 5 – 50 81 .8 21 – 40 .8 15 – 49 6. 73 4 – −8 83 .4 38 – 17 80 .6 11 – B IC −1 83 5. 95 5 – 52 64 .0 90 – 88 .9 49 – 54 4. 86 9 – −7 89 .1 92 – 18 74 .8 57 – Q L B (4 ) 0. 99 5 – 0. 99 8 – 0. 92 2 – 0. 90 2 – 0. 99 9 – 0. 99 6 – Q L B (8 ) 0. 95 1 – 1. 00 0 – 0. 68 4 – 0. 47 8 – 0. 11 2 – 0. 90 9 – Q L B (1 2) 0. 75 9 – 0. 98 7 – 0. 78 5 – 0. 62 5 – 0. 22 4 – 0. 75 7 – H yp ot he si s de sc ri pt io n F –s ta ti st ic P va lu e F –s ta ti st ic P va lu e F –s ta ti st ic P va lu e F –s ta ti st ic P va lu e F –s ta ti st ic P va lu e F –s ta ti st ic P va lu e G ra ng er c au sa lit y te st 01 : 0 yi yj H   + − = = 16 .8 99 ** * 0. 00 0 5. 72 4* ** 0. 00 0 2. 34 ** 0. 02 9 1. 81 1* 0. 09 1 2. 93 8* ** 0. 00 1 0. 95 3 0. 48 5 02 : 0 xi xj H   + − = = 4. 94 1* ** 0. 00 0 5. 13 4* ** 0. 00 0 3. 09 8* ** 0. 00 5 0. 78 1 0. 62 1 0. 73 6 0. 69 0 3. 86 9* ** 0. 00 0 D is tr ib ut ed la g as ym m et ri c ef fe ct 03 1 1 : 0 x x H   + − = = 11 .5 05 ** * 0. 00 1 5. 40 3* * 0. 02 0 4. 86 2* * 0. 03 1 0. 54 8 0. 46 2 7. 98 2* ** 0. 00 5 0. 59 3 0. 44 2 03 2 2 : 0 x x H   + − = = 2. 95 4* 0. 08 6 10 .4 44 ** * 0. 00 1 0. 04 2 0. 83 9 0. 36 6 0. 54 8 0. 02 6 0. 87 3 0. 00 0 0. 99 3 03 3 3 : 0 x x H   + − = = 0. 10 0 0. 75 2 1. 26 1 0. 26 2 0. 03 8 0. 84 6 0. 10 7 0. 74 5 1. 99 3 0. 15 9 2. 62 5 0. 10 6 03 4 4 : 0 x x H   + − = = 2. 55 8 0. 11 0 6. 98 2* ** 0. 00 8 0. 01 3 0. 90 9 2. 94 2* 0. 09 1 0. 05 3 0. 81 8 0. 01 5 0. 90 2 Ta bl e 9: (C on tin ue d) Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 103 V ar ia bl e W ho le s am pl e (C .M − TA R ; l ag =8 ) In –c ri si s (C . T A R ; l ag =4 ) P os t– cr is is (C . T A R ; l ag =5 ) ST L F SI W T I ST L F SI W T I ST L F SI W T I C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c C oe ffi ci en t t– st at is ti c 03 5 5 : 0 x x H   + − = = 27 .5 26 ** * 0. 00 0 10 .5 62 ** * 0. 00 1 – – – – 0. 07 4 0. 78 6 0. 43 1 0. 51 2 03 6 6 : 0 x x H   + − = = 9. 92 7* ** 0. 00 2 0. 77 6 0. 37 9 – – – – – – – – 03 7 7 : 0 x x H   + − = = 15 .5 28 ** * 0. 00 0 0. 17 7 0. 67 4 – – – – – – – – 03 8 8 : 0 x x H   + − = = 3. 30 0* 0. 07 0 3. 21 3* 0. 07 3 – – – – – – – – 04 1 1 : 0 y y H   + − = = 9. 80 8* ** 0. 00 2 2. 84 0* 0. 09 2 1. 98 5 0. 16 4 2. 72 4 0. 10 4 1. 30 4 0. 25 4 1. 49 3 0. 22 3 04 2 2 : 0 y y H   + − = = 2. 69 9 0. 10 1 2. 20 4 0. 13 8 0. 00 0 0. 98 5 1. 01 5 0. 31 8 2. 58 1 0. 10 9 4. 62 2* * 0. 03 2 04 3 3 : 0 y y H   + − = = 13 .0 22 ** * 0. 00 0 1. 10 4 0. 29 4 11 .2 26 ** * 0. 00 1 0. 25 6 0. 61 5 0. 00 8 0. 93 0 0. 51 9 0. 47 2 04 4 4 : 0 y y H   + − = = 0. 57 7 0. 44 8 1. 01 4 0. 31 4 0. 52 0. 47 4 0. 08 1 0. 77 7 0. 43 3 0. 51 1 2. 45 1 0. 11 8 04 5 5 : 0 y y H   + − = = 0. 04 9 0. 82 5 1. 22 8 0. 26 8 – – – – 0. 13 8 0. 71 1 0. 04 0 0. 84 2 04 6 6 : 0 y y H   + − = = 1. 15 7 0. 28 2 2. 77 3* 0. 09 6 – – – – – – – – 04 7 7 : 0 y y H   + − = = 1. 32 1 0. 25 1 0. 97 4 0. 32 4 – – – – – – – – 04 8 8 : 0 y y H   + − = = 21 .0 07 ** * 0. 00 0 2. 69 0 0. 10 1 – – – – – – – – C um ul at iv e as ym m et ri c ef fe ct 05 1 1 : p p xj xj j j H   + − = = = ∑ ∑ 9. 96 3* ** 0. 00 2 0. 06 8 0. 79 4 1. 20 4 0. 27 7 0. 45 3 0. 50 3 6. 14 9* * 0. 01 4 1. 01 2 0. 31 5 06 1 1 : p p yj yj j j H   + − = = = ∑ ∑ 2. 20 6 0. 13 8 0. 54 6 0. 46 0 1. 44 7 0. 23 3 1. 14 9 0. 28 8 0. 00 2 0. 96 7 2. 04 6 0. 15 3 Eq ui lib riu m ad ju stm en t p at h as ym m et ry H 07 : δ + = δ− 1. 70 6 0. 19 2 7. 95 2* ** 0. 00 5 1. 89 3 0. 17 4 0. 37 7 0. 54 1 0. 18 9 0. 66 4 3. 27 4* 0. 07 1 Q LB (P ) d en ot es th e si gn ifi ca nc e le ve l f or th e L ju ng -B ox Q s ta tis tic , w hi ch te st s se ri al c or re la tio n ba se d on P a ut oc or re la tio n co ef fic ie nt s. N um be rs in p ar en th es es a re P v al ue s. * ** , * *, a nd * de no te s ta tis tic al s ig ni fic an ce a t t he 1 % , 5 % a nd 1 0% le ve ls , r es pe ct iv el y Ta bl e 9: (C on tin ue d) Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019104 500 weeks to be fully digested. Therefore, in the short-term, WTI has a much faster reaction to positive deviations from long-term equilibrium than negative deviations. 5. CONCLUDING REMARKS In this article, we investigated the dynamic relationship between oil prices and financial stress over the period from December 31, 1993 to July 15, 2016. In particular, we focused on the linkages between those variables in both the long-run and short-run horizons under both the linear and nonlinear (threshold) cointegration framework. As an extension of preceding studies, we make use of the methodology developed by Enders and Siklos (2001), based on a nonlinear (threshold) cointegration model allowing for nonlinear adjustment to long-run equilibrium. The results from the conventional linear cointegration approaches suggested that we can reject the null hypothesis of no cointegration. According to our empirical results, the null hypothesis of linear cointegration between oil prices and financial stress is rejected in favor of a threshold cointegration model in the sense that the short- term adjustments to the equilibrium are asymmetric depending on the deviation from the equilibrium. Thereafter, using the consistent TAR and MTAR specifications, we found evidence of asymmetry in the adjustment process to equilibrium. These findings would suggest the presence of a significant nonlinear behavior in the oil price-financial stress relationship. We find an asymmetric effect in the short-term adjustment process. Regimes with negative (below the threshold) changes of deviations adjust much quicker than regimes with positive (above the threshold) changes of deviations, especially during the crisis period. After incorporating the asymmetric adjustment and using a Granger causality test, we find a bi-directional causality between oil prices and financial stress index, indicating that these variables affect have been evolving more dependently in the short term, in particular during the crisis period. Our results also reveal the existence of three types of statistical significant asymmetries, namely the distributed lag asymmetric effect, the cumulative asymmetric effects, and the momentum equilibrium adjustment path asymmetries. Future research could extend the two-regime threshold cointegration model to three or more regimes. REFERENCES Balke, N.S., Brown, S.P.A., Yücel, M.K. (2002), Oil price shocks and the US economy: Where does the asymmetry originate? Energy Journal, 23(3), 27-52. Balke, N.S., Fomby, T.B. (1997), Threshold cointegration. International Economic Review, 38(3), 627-645. Baur, D.G. (2012), Financial contagion and the real economy. Journal of Banking and Finance, 36(10), 2680-2692. BIS. (2009), The International Financial Crisis: Timeline, Impact and Policy Responses in Asia and the Pacific. Basel: Bank for International Settlements. Caner, M., Hansen, B. (2001), Threshold autoregression with a unit root. Econometrica, 69(6), 1555-1596. Cardarelli, R., Elekdag, S., Lall, S. (2011), Financial stress and economic contractions. Journal of Financial Stability, 7(2), 78-97. Cevik, E.I., Dibooglu, S., Kenc, T. (2013), Measuring financial stress in Turkey. Journal of Policy Modeling, 35(2), 370-383. Chan, K.S. (1993), Consistency and limiting distribution of the least squares estimator of a threshold autoregressive model. Annals of Statistics, 21(1), 520-533. Chau, F., Deesomsak, R. (2014), Does linkage fuel the fire? The transmission of financial stress across the markets. International Review of Financial Analysis, 36, 57-70. Chen, H., Choi, P.M.S., Hong, Y. (2013), How smooth is price discovery? Evidence from cross-listed stock trading. Journal of International Money and Finance, 32, 668-699. Chen, H., Zhu, Y. (2015), An empirical study on the threshold cointegration of Chinese A and H cross-listed shares. Journal of Applied Statistics, 42(11), 2406-2419. Chen, W., Hamori, S., Kinkyo, T. (2014), Macroeconomic impacts of oil prices and underlying financial shocks. Journal of International Financial Markets, Institutions and Money, 29, 1-12. Cunado, J., de Gracia, F.P. (2003), Do oil price shocks matter? Evidence for some European countries. Energy Economics, 25(2), 137-154. Cunado, J., de Gracia, F.P. (2005), Oil prices, economic activity, and inflation: Evidence for some Asian countries. The Quarterly Review of Economics and Finance, 45(1), 65-83. Davig, T., Hakkio, C. (2010), What is the effect of financial stress on economic activity? Federal reserve bank of Kansas city. Economic Review, 2010, 35-62. Demirguc-Kunt, A., Levine, R. (2001), Financial Structures and Economic Growth: A Cross Country Comparison of Banks, Markets, and Development. Cambridge, MA: MIT Press. Dickey, D.A., Fuller, W.A. (1979), Distribution of the estimators for autoregressive time series with a unit root. Journal of the American Statistical Association, 74(366), 427-431. Dimitriou, D., Kenourgios, D. (2013), Financial crises and dynamic linkages among international currencies. Journal of International Financial Markets, Institutions and Money, 26, 319-332. Dimitriou, D., Kenourgios, D., Simos, T. (2013), Global financial crisis and emerging stock market contagion: A multivariate FIAPARCH- DCC approach. International Review of Financial Analysis, 30, 46-56. Enders, W. (2004), Applied Econometric Time Series. New York: John Wiley and Sons, Inc. p480. Enders, W., Granger, C.W.F. (1998), Unit-root tests and asymmetric adjustment with an example using the term structure of interest rates. Journal of Business and Economic Statistics, 16(3), 304-311. Enders, W., Siklos, P.L. (2001), Cointegration and threshold adjustment. Journal of Business and Economic Statistics, 19(2), 166-176. Engle, R., Granger, C.W.J. (1987), Cointegration and error correction: Representation, estimation, and testing. Econometrica, 55(2), 251-276. Federal Reserve Board of St. Louis. (2009), The Financial Crisis: A Timeline of Events and Policy Actions. St. Louis: Federal Reserve Board of St. Louis. Frey, G., Manera, M. (2007), Econometric models of asymmetric price transmission. Journal of Economic Surveys, 21(2), 349-415. Granger, C.W.J., Lee, T.H. (1989), Investigation of production, sales, and inventory relationships using multicointegration and non- symmetric error correction models. Journal of Applied Econometrics, 4, 145-159. Gregory, A.W., Hansen, B.E. (1996a), Residual-based tests for cointegration in models with regime shifts. Journal of Econometrics, 70(1), 99-126. Gregory, A.W., Hansen, B.E. (1996b), Tests for cointegration in models with regime and trend shifts. Oxford Bulletin of Economics and Statistics, 58(3), 555-560. Hakkio, C., Keeton, W. (2009), Financial stress: What is it, how can it be Mighri and Alsaggaf: Asymmetric Threshold Cointegration and Nonlinear Adjustment between Oil Prices and Financial Stress International Journal of Energy Economics and Policy | Vol 9 • Issue 3 • 2019 105 measured, and why does it matter? Federal reserve bank of Kansas city. Economic Review, 2009, 5-50. Hamilton, J. (2011), Nonlinearities and the macroeconomic effects of oil prices. Macroeconomic Dynamics, 15(3), 364-378. Holló, D. (2012), A System-Wide Financial Stress Indicator for the Hungarian Financial System. MNB Occasional Papers No. 105. Holló, D., Kremer, M., Lo Duca, M. (2012), CISS-A B Composite Indicator of Systemic Stress in the Financial System. European Central Bank Working Paper No. 1426. Illing, M., Liu, Y. (2006), Measuring financial stress in a developed country: An application to Canada. Journal of Financial Stability, 2, 243-265. Islami, M., Kurz-Kim, J.R. (2014), A single composite financial stress indicator and its real impact in the Euro area. International Journal of Finance and Economics, 19, 204-211. Johansen, S. (1988), Statistical analysis of cointegration vectors. Journal of Economic Dynamics and Control, 12(2-3), 231-254. Johansen, S., Juselius, K. (1990), Maximum likelihood estimation and inference on cointegration with applications to the demand for money. Oxford Bulletin of Economics and Statistics, 52(2), 169-210. Kilian, L. (2008), The economic effects of energy price shocks. Journal of Economic Literature, 46(4), 871-909. Killian, L. (2009), Not all oil price shocks are alike: Disentangling demand and supply shocks in the crude oil market. American Economic Review, 99(3), 1053-1069. Kliesen, K.L., Owyang, M.T., Vermann, E.K. (2012), Disentangling diverse measures: A survey of financial stress indexes. Federal reserve bank of St. Louis. Economic Review, 94(5), 369-397. Kuo, S.H., Enders, W. (2004), The term structure of Japanese interest rates: The equilibrium spread with asymmetric dynamics. Journal of the Japanese and International Economies, 18(1), 84-98. Kwiatkowski, D., Phillips, P.C.B., Schmidt, P., Shin, Y. (1992), Testing the null hypothesis of stationarity against the alternative of a unit root. Journal of Econometrics, 54(1-3), 159-178. Levine R. (2005), Finance and growth: Theory and evidence. In: Aghion, P., Durlauf, S., editors. Handbook of Economic Growth. Vol. 1. Amsterdam: North-Holland Elsevier. p865-934. MacKinnon, J.G. (1996), Numerical distribution functions for unit root and cointegration tests. Journal of Applied Econometrics, 11(6), 601-618. Mallick, S.K., Sousa, R.M. (2013), The real effects of financial stress in the Eurozone. International Review of Financial Analysis, 30, 1-17. Meyer, J., Von Cramon-Taubade, S. (2004), Asymmetric price transmission: A survey. Journal of Agricultural Economics, 55(3), 581-611. Mighri, Z., Mansouri, F. (2014), Modeling international stock market contagion using multivariate fractionally integrated APARCH approach. Cogent Economics and Finance, 2(1), 963632. Mighri, Z., Mansouri, F. (2016), Asymmetric price transmission within the Argentinean stock market: An asymmetric threshold cointegration approach. Empirical Economics, 51(3), 1115-1149. Mollick, A.V., Assefa, T.A. (2013), U.S. stock returns and oil prices: The tale from daily data and the 2008-2009 financial crisis. Energy Economics, 36, 1-18. Morales, M.A., Estrada, D. (2010), A financial stability index for Columbia. Annals of Finance, 6(4), 555-581. Nazlioglu, S., Soytas, U., Gupta, R. (2015), Oil prices and financial stress: A volatility spillover analysis. Energy Policy, 82, 278-288. Phillips, P.C.B., Perron, P. (1988), Testing for a unit root in time series regression. Biometrika, 75(2), 335-346. Rafiq, J., Salim, R., Bloch, H. (2009), Impact of crude oil price volatility on economic activities: An empirical investigation in the Thai economy. Resources Policy, 34(3), 121-132. Slingenberg, J.W., de Haan, J. (2011), Forecasting Financial Stress. De Nederlandsche Bank Working Paper No. 292. Sun, C. (2011), Price dynamics in the import wooden bed market of the United States. Forest Policy and Economics, 13(6), 479-487. Thompson, M.A. (2006), Asymmetric adjustment in the prime lending- deposit rate spread. Review of Financial Economics, 15(4), 323-329. Tong, H. (1983), Threshold Models in Non-Linear Time Series Analysis. New York: Springer-Verlag. Turhan, I., Hacıhasanoglu, E., Soytas, U. (2013), Oil prices and emerging market exchange rates. Emerging Markets Finance and Trade, 49(1), 21-36. Vermeulen, R., Hoeberichts, M., Vašíček, B., Žigraiová, D., Šmídková, K., de Haan, J. (2015), Financial stress indices and financial crises. Open Economic Review, 26(3), 383-406.