Next Article in Journal
Introduction of a Corporate Security Risk Management System: The Experience of Poland
Next Article in Special Issue
Pricing Cat Bonds for Cloud Service Failures
Previous Article in Journal
Health Insurance Patterns of Older Veterans: Evidence from the Health and Retirement Study
Previous Article in Special Issue
Cap Rates as a Function of Real Economic Growth
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Co-Jumps, Co-Jump Tests, and Volatility Forecasting: Monte Carlo and Empirical Evidence

1
Department of Finance, Jack Welch College of Business & Technology, Sacred Heart University, 5151 Park Avenue, Fairfield, CT 06825, USA
2
US Economics Research, Barclays, 745 Seventh Avenue, New York, NY 10019, USA
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2022, 15(8), 334; https://doi.org/10.3390/jrfm15080334
Submission received: 29 June 2022 / Revised: 24 July 2022 / Accepted: 25 July 2022 / Published: 28 July 2022
(This article belongs to the Special Issue Mathematical and Empirical Finance)

Abstract

:
This study classifies jumps into idiosyncratic jumps and co-jumps to quantitatively identify systematic risk and idiosyncratic risk by utilizing high-frequency data. We found that systematic risk occurs more frequently and has larger magnitudes than the idiosyncratic risk in an individual asset, which indicates that volatilities from one sector are largely derived from the contagious effect of other sectors. We further investigated the importance of idiosyncratic jumps and co-jumps to predict the sector-level S&P500 exchange-traded fund (ETF) volatility. It was found that the predictive content of co-jumps is higher than that of idiosyncratic jumps, suggesting that systematic risk is more informative than idiosyncratic risk in volatility forecasting. Additionally, we carried out Monte Carlo experiments designed to examine the relative performances of the four co-jump tests. The findings indicate that the BLT test and the co-exceedance rule of the LM test outperform other tests, while the co-exceedance rule of the LM test has larger power and a smaller empirical size than that of the BLT test.

1. Introduction

Effective stock volatility prediction remains a key goal for empirical finance researchers. For example, in portfolio management, agents utilize such predictions when allocating assets and managing risk. Various discrete-time models (e.g., ARCH-GARCH models) and continuous-time models (e.g., stochastic volatility and related models) have been developed to address this issue. Key papers in the discrete-time framework include Harvey et al. (1994), Andersen and Lund (1997), Hansen and Lunde (2005), Brandt and Jones (2006), and the references therein. With the availability of high-frequency intraday data, Andersen et al. (2001) developed a nonparametric volatility measure, which fostered the development of a number of time series volatility forecasting models, including mixed data sampling (MIDAS) models (see Ghysels et al. 2006; Ghysels and Sinko 2011) and long-memory ARMA models (e.g., Andersen et al. 2003; Koopman et al. 2005). In particular, this has led to a new class of parsimonious models, called heterogeneous autoregressive models of realized volatility (HAR-RV), including lag terms of RV over different time horizons. The original HAR forecasting model, introduced by Müller et al. (1997) and Corsi (2009), has been augmented by separating jump and diffusion components, including a variety of jump variation and intensity-related variables (e.g., Andersen et al. 2007; Corsi et al. 2010; Patton and Sheppard 2015; Lahaye and Neely 2020). These models have been successfully applied in various empirical applications.
In this study, we classified jumps into idiosyncratic jumps and co-jumps to decompose jump risk into idiosyncratic risk and systematic risk, and further investigated the importance of idiosyncratic jumps and co-jumps to predict sector-level S&P500 exchange-traded fund (ETF) volatility. Previous literature mostly focused on studying jumps in individual assets or a portfolio index and examining the predictability of jumps in forecasting volatilities of different target assets without decomposition. Our motivation for decomposing jumps into idiosyncratic jumps and co-jumps is to quantitatively identify idiosyncratic risk and systematic risk in the individual asset by utilizing 5-min price data and, thus, evaluate the predictability of idiosyncratic risk and systematic risk in forecasting volatilities of the target asset. Our impetus for this decomposition is the hypothesis that systematic jumps (risk) may be more informative than idiosyncratic jumps (risk) in volatility prediction. In addition, the identification of systematic jumps from idiosyncratic jumps is important on its own. In theory, idiosyncratic jumps can be alleviated in a weighted portfolio, whereas systematic jumps cannot. Hence, the failure to distinguish idiosyncratic jumps (idiosyncratic risk) from co-jumps (systematic risk) can lead to dire consequences when managing portfolios.
To evaluate this hypothesis, co-jumps and idiosyncratic jumps are included in predictive HAR-type regressions, which enables us to assess their relative marginal predictive content. Specifically, we compare the following alternative models: (i) the basic HAR model, which includes the continuous component of realized volatility; (ii) the HAR model that includes the continuous component and the co-jump component; (iii) the augmented HAR model that includes the continuous component and the idiosyncratic jump component; and (iv) the augmented HAR model that includes all three components. We predict the volatility of a variety of target assets, including the major S&P500 sector-level ETFs: financial (XLF), technology (XLK), consumer staple (XLP), utility (XLU), and consumer discretionary (XLY). Co-jumps and idiosyncratic jumps among sector-level ETFs are detected using the co-exceedance rule of the LM jump test in Lee and Mykland (2007) and Gilder et al. (2014).
Our experiment focused on one-day-ahead volatility predictions using rolling windows. Our sample period was from November 2008 to December 2019. For model evaluation, we used in-sample adjusted R 2 , out-of-sample R 2 , and mean absolute forecasting error (MAFE). In addition, we implement the model confidence set (MCS) tests of Hansen et al. (2011) to compare the predictive abilities of all candidate models. We summarize our main findings below.
First, among the four alternative HAR models in our forecasting experiment, models augmented with co-jumps and with both co-jumps and idiosyncratic jumps showed the largest predictive accuracy gains, in terms of in-sample adjusted R 2 , out-of-sample R 2 , and the mean absolute forecasting error (MAFE). This result was robust across different forecasting targets, i.e., S&P 500 sector ETF volatilities. For example, models augmented with co-jumps showed smaller MAFEs than the baseline HAR model and HAR model augmented with idiosyncratic jumps for all five target sectors, as listed in Table 10. Moreover, models combining both co-jumps and idiosyncratic jumps showed even smaller MAFEs, except for the financial sector. Moreover, the model confidence set (MCS) test, which assigns the value of 1 to the best model and a larger number for superior models, also suggests that models augmented with co-jumps or models augmented with both co-jumps and idiosyncratic jumps are superior to other models.
Second, models with co-jumps resulted in higher explanatory power than models with idiosyncratic jumps in both the in-sample and ex-ante forecasting experiments, in line with the hypothesis that systematic risks (as represented by co-jumps) potentially play a larger role than idiosyncratic risks in volatility forecasting.
Third, co-jumps occurred generally two to three times more frequently than idiosyncratic jumps during our sample period from 2007 to 2019, and the frequency of occurrences of both the co-jumps and idiosyncratic jumps increased especially during significant market turmoil, such as the 2008 financial crisis and the 2011 European debt crisis. This finding is consistent with the prevalence of systematic risk compared with the idiosyncratic risk in the financial market.
Fourth, we also found that co-jumps have a larger magnitude and variation compared to idiosyncratic jumps, suggesting systematic risk may pose a larger impact on the financial market than idiosyncratic risk.
Finally, given the importance of co-jump tests in our empirical analysis, we additionally carried out a series of Monte Carlo experiments to examine the relative performances of a number of widely used co-jump tests, including: (i) the BLT co-jump test in Bollerslev et al. (2008); (ii) the JT co-jump test in Jacod and Todorov (2009); (iii) the co-exceedance rule in Gilder et al. (2014) using the LM test in Lee and Mykland (2007); and (iv) the co-exceedance rule using the intersection of the LM test and BNS test in Barndorff-Nielsen and Shephard (2006). While univariate jump tests have been researched extensively1, the co-jump testing literature used in our empirical analysis is relatively nascent. One strand of the literature identifies co-jumps based on the identification of jumps in a portfolio (e.g., Bollerslev et al. 2008). Another strand uses univariate jump tests to identify co-jumps in multivariate processes (e.g., Gilder et al. 2014). A third strand develops co-jump tests by direct examination of multiple price processes (see Jacod and Todorov 2009; Mancini and Gobbi 2012; Bandi and Reno 2016; Bibinger and Winkelmann 2015; Caporin et al. 2017).
Therefore, based on our Monte Carlo experiments, our key finding is that the BLT test and the co-exceedance rule of the LM test outperform other tests, while the co-exceedance rule of the LM test has larger power and a smaller empirical size than that of the BLT test. The JT test, the BLT test, and the co-exceedance rule of the LM test all have good powers, except for the co-exceedance rule of the intersection of the LM and BNS tests, and this result is robust across all the data-generating processes (DGPs) considered in our experiment. Additionally, the BLT test and the co-exceedance rule of the LM test both have empirical sizes close to the nominal size in all DGPs. The empirical size of the JT test is mixed, under- or over-sized in different DGPs, and the co-exceedance rule of the LM and BND test is significantly undersized. Interestingly, adding idiosyncratic jumps into the DGPs increases the size of the JT test and lowers the size of the other three tests.
The rest of this paper is organized as follows. Section 2 outlines the setup and summarizes the co-jump tests examined in the paper. Section 3 discusses the Monte Carlo experiments. Section 4 outlines all volatility forecasting models used in the sequel. Section 5 summarizes the setup of our forecasting experiments. Section 6 briefly discusses the data used in our analysis, and Section 7 summarizes our key empirical findings. Concluding remarks are contained in Section 9.

2. Co-Jump Identification

2.1. Setup

Let P t be the log price of an asset at time t. We assume P t follows an Itô semimartingale:
P t = P 0 + 0 t b s d s + 0 t σ s d W s + J t
In the above equation, both b s and σ s are adapted, càdlàg, and locally bounded. W s is a standard Brownian motion and J t is a compound Poisson process (CPP):
J t = i = 1 N t Y i
where N t is a Poisson process, representing the number of jumps in the interval [ 0 , t ] . The jump magnitudes Y i ’s are i i d random variables.
Consider a finite time horizon [ 0 , t ] that contains n high-frequency observations of the log price process. A typical time horizon is one day. Let Δ n = t / n be the sampling frequency. The intraday return at the equidistant interval is [ ( i 1 ) Δ n , i Δ n ] is r i = P i Δ n P ( i 1 ) Δ n .
In the high-frequency literature, it is widely recognized that realized volatility (RV)2 is a consistent estimator of quadratic variation, namely:
R V t = i = 1 n r i 2 u . c . p . 0 t σ s 2 d s + s T ( Δ P s ) 2 = Q V t = I V t + J V t ,
where u.c.p. denotes convergence in probability, uniformly in time. The integrated volatility ( I V t = 0 t σ s 2 d s ) measures the contribution of the continuous component to the total quadratic variation. The jump variation J V t = s T ( Δ P s ) 2 measures the contribution from the discontinuous (jump) component.
There are many estimators of IV. For example, Barndorff-Nielsen and Shephard (2004) propose a bipower variation estimation B P V t to estimate integrated volatility. Namely, define the standardized realized bipower variation as:
B P V t = ( μ 1 ) 2 i = 2 n | r i | | r i 1 | ,
where μ 1 = E ( | Z | ) = 2 Γ ( 1 ) / Γ ( 1 / 2 ) = 2 / π , where Z is a standard normal random variable and Γ ( · ) denotes the gamma function.

2.2. Co-Jump Tests

In this section, we briefly summarize the co-jump testing methods used in the sequel.

2.2.1. BLT Co-Jump Test

Bollerslev et al. (2008) proposed a BLT test to detect co-jumps in a large ensemble of stocks. They developed a theoretical foundation that showed how only co-jumps (not idiosyncratic jumps) could be detected in a large equal-weighted index. Let M denote the total number of assets under consideration. The BLT mean cross-product test statistic is defined as:
m c p t , i = 2 M ( M 1 ) j = 1 M 1 l = j + 1 M r i j r i l , i = 1 , , n and t = 1 , , T ,
where r i j = P i Δ n j P ( i 1 ) Δ n j for j = 1 , , M .
The mcp statistic has a nonzero means; the corresponding studentized test statistic is:
z m c p , t , i = m c p t , i m c p t ¯ s m c p , t , for i = 1 , , n a n d t = 1 , , T .
where
m c p t ¯ = 1 n m c p t = 1 n i = 1 n m c p t , i ,
s m c p , t = 1 n i = 1 n ( m c p t , i m c p t ¯ ) 2 .
The critical value under the null is calculated by bootstrapping the test statistics z m c p , t , i using Monte Carlo simulations.

2.2.2. JT Co-Jump Test

Jacod and Todorov (2009) constructed two test statistics to identify co-jumps under two different null hypotheses of joint and disjoint jumps. The test statistics were proposed for detecting co-jumps in a bivariate process. Co-jumps among multivariate processes can be detected from the combination of bivariate processes. The test statistics of the common jump Φ n ( j ) and disjoint jump Φ n ( d ) are defined as:
Φ n ( j ) = V ( f , k Δ n ) t V ( f , Δ n ) t and Φ n ( d ) = V ( f , Δ n ) t V ( g 1 , Δ n ) t V ( g 2 , Δ n ) t ,
where k is an integer greater than 1 and V ( f , k Δ n ) t is defined as:
V ( f , k Δ n ) t = i = 1 [ t / k Δ n ] f ( P i Δ n P ( i 1 ) Δ n ) .
The functions f ( x ) , g 1 ( x ) , and g 2 ( x ) used in the cited paper are given as follows:
f ( x ) = ( x 1 x 2 ) 2 , g 1 ( x ) = ( x 1 ) 4 , g 2 ( x ) = ( x 2 ) 4 .
They show that the test statistics for the null hypothesis of disjoint jumps Φ n ( d ) converge in probability to 0 on Ω T ( d ) and the null hypothesis with common jumps Φ n ( j ) converges in probability to 1 on Ω T ( j ) (refer to Theorem 3.1 therein). Here, Ω T ( j ) and Ω T ( d ) are defined as:
Ω T ( j ) = { ω : on [ 0 , t ] the process r i 1 r i 2 is not identically 0 } ,
Ω T ( d ) = { ω : on [ 0 , t ] the processes r i 1 and r i 2 are not identically 0 , but the process r i 1 r i 2 is } .
The authors constructed critical regions of the two statistics as:
C n ( j ) = { | Φ n ( j ) 1 | c n ( j ) } and C n ( d ) = { Φ n ( d ) c n ( d ) } ,
where c n ( j ) = V ^ n ( j ) / α and c n ( d ) = V ^ n ( d ) / α with significance level α . The term V ^ n ( j ) is defined as:
V ^ n ( j ) = Δ n ( k 1 ) F ^ t n V ( f , Δ n ) t
with F ^ t n given by:
F ^ t n = 2 k n Δ n i = 1 + k n [ t / Δ n ] k n 1 m I n ( i ) ( r i 1 ) 2 ( r i 2 ) 2 × ( r i 1 r m 2 + r m 1 r i 2 ) 2 1 { r i > α Δ n ω , r m α Δ n ω } ,
where I n ( i ) = I n , ( i ) I n , + ( i ) and I n , ( i ) = { i k n , i k n + 1 , , i 1 } if i > k n and I n , + ( i ) = { i + 2 , i + 3 , , i + k n + 1 } . The term V ^ n ( d ) is defined as:
V ^ n ( d ) = Δ n F ^ t n + A ^ T n V ( g 1 , Δ n ) t V ( g 2 , Δ n ) t
where
F ^ t n = 2 k n Δ n i = 1 + k n [ t / Δ n ] k n 1 m I n ( i ) ( ( r i 1 ) 2 ( r m 2 ) 2 + ( r m 1 ) 2 ( r i 2 ) 2 ) × 1 { r i > α Δ n ω , r m α Δ n ω }
A ^ t n = 1 Δ n i = 1 t / Δ n f ( r i ) 1 { r i α Δ n ω } .
The local window k n = 1 / Δ n and the truncation level of α Δ n ω = 0.03 × Δ n 0.49 3.

2.2.3. Univariate Jump Co-Exceedance Rule

Gilder et al. (2014) proposed a co-exceedance-based co-jump detection method by applying univariate jump tests to individual stocks to identify co-jumps. They selected three univariate jump tests developed by Barndorff-Nielsen and Shephard (2006); Lee and Mykland (2007) and Andersen et al. (2010), respectively. The idea was to find the intersection between ABD jump test results and BNS jump test results (ABD∩BNS), the intersection between ABD jump test results and LM jump test results (ABD∩LM), the intersection between BNS jump test results and LM jump test results (BNS∩LM), and the intersection among three jump tests results (ABD∩LM∩BNS), and retain only the intersection of the identified jumps to avoid spurious test results. Co-jumps are detected if more than one univariate process has jumps at the same time interval.
In the Lee and Mykland (2007) jump test, they used the ratio of realized returns to the estimated instantaneous volatility and constructed a nonparametric jump test to identify the exact timings of the jumps at the intraday level. The test statistics to identify the jump during interval ( t + l / n , t + ( l + 1 ) / n ) is:
L ( t + ( l + 1 ) / n ) = X t + ( l + 1 ) / n X t + l / n σ ^ t + ( l + 1 ) / n ,
where
σ ^ t + ( l + 1 ) / n 2 1 K 2 i = l K + 1 l 2 | X t + ( i + 1 ) / n X t + i / n X t + i / n X t + ( i 1 ) / n | .
Here, K is the window size of a local movement of the process. It is also shown that:
max l A ¯ n | L ( t + ( l + 1 ) / n ) | C n S n ε , a s Δ n 0 ,
where ε has a cumulative distribution function P ( ε x ) = exp ( e x ) ,
C n = ( 2 log n ) 1 / 2 c log π + log ( log n ) 2 c ( 2 log n ) 1 / 2 and S n = 1 c ( 2 log n ) 1 / 2
c ≈ 0.7979 and A ¯ n is the set of l { 0 , 1 , , n } , so that there are no jumps in ( t + l / n , t + ( l + 1 ) / n ] .
Barndorff-Nielsen and Shephard (2006) developed a non-parametric (Hausman 1978) type test using the difference between the realized quadratic variation and realized bipower variation. Realized quadratic variation was estimated using realized volatility. The adjusted jump ratio test statistics can be defined as,
B N S = M 1 / 2 ϑ max ( 1 , Q P V ( μ 1 2 B P V ) 2 ) ( 1 B P V R V ) d N ( 0 , 1 )
where R V and B P V are the same as in (3) and (4), ϑ = π 2 / 4 + π 5 0.6090 . The realized quad-power variation Q P V is used to estimate integrated quarticity ( 0 t σ s 4 d s ) and can be given as:
Q P V = 1 Δ n j = 4 n | Δ j X | | Δ j 1 X | | Δ j 2 X | | Δ j 3 X | P μ 1 4 0 t σ s 4 d s .
The authors show that the null hypothesis of no jumps is rejected if the test statistic B N S is in excess of the critical value Φ α , where α is the test significance level, leading to a conclusion that there are jumps during the tested time interval.

3. Monte Carlo Experiments

The co-jump tests discussed above are widely used in the empirical literature to identify co-jumps among asset prices. In this section, we carry out a series of Monte Carlo experiments to evaluate the finite sample properties of the above co-jump tests: (i) BLT co-jump test; (ii) JT co-jump test; (iii) co-exceedance rule using the LM jump test; and (iv) the co-exceedance rule using the intersection of the LM test and BNS test. It is possible that both co-jumps and idiosyncratic jumps are present in the same period. Hence, it is important to check which test has the best performance in identifying co-jumps from idiosyncratic jumps.
We generate log prices for M = 10 assets using a Milstein discretion scheme from the following data-generating process (DGP):
d P t ( i ) = μ d t + V t ( i ) d W 1 , t ( i ) + d C J t ( i ) + d I J t ( i ) ,
d V t ( i ) = κ v ( θ v V t ( i ) ) d t + ζ V t ( i ) d W 2 , t ( i ) ,
where W 1 , t ( i ) and W 2 , t ( i ) are two standard Brownian motions with correlation ρ , where ρ = { 0 , 0.5 } . The co-jump and idiosyncratic jump are denoted by C J t ( i ) and I J t ( i ) , respectively. The jump component J t ( i ) (the sum of C J t ( i ) and I J t ( i ) ) are simulated as compound Poisson processes N t with intensity λ i :
J t ( i ) = j = 1 N t Y j ,
where Y j ’s are i i d normal, which characterizes the jump size. Under the null hypothesis, we first simulate DGPs without jumps as the first scenario. Then, we add idiosyncratic jumps I J t ( i ) onto the sample paths as the second scenario of the null hypothesis. Under the alternative, we add co-jumps C J t ( i ) to the null DGPs. When adding co-jumps to the null DGPs, we set the number of co-jumps C J t ( i ) among 10 assets equal to m = 3 , which is in line with what we find in real data4. Table 1 shows the parameter settings in the Monte Carlo experiments. We follow the parameter settings from (Huang and Tauchen 2005; Aït-Sahalia et al. 2009; Corradi et al. 2018). The jump intensity λ is set to be either 0.1 or 0.2. The sampling frequency is Δ n = 1 78 , which is consistent with 5-min sampling frequencies. Table 2 shows details of the data generating processes. In all experiments, we performed 1000 Monte Carlo replications.
The co-exceedance rule using the LM jump test and the intersection of the LM test and BNS test are implemented using the following approach5,
i = 1 M 1 { J u m p s ( i ) > 0 } 2 , co - jump 1 , no co - jump
where 1 { J u m p s ( i ) > 0 } is an indicator function equals 1 when a jump is detected in asset i during the intraday interval that contains time s, and takes value 0 when no jump is detected.
Since the JT test is a bivariate co-jump test, we use the following rule to detect co-jumps:
1 i < j M 1 { C o j u m p s ( i , j ) > 0 } 1 , co - jump = 0 , no co - jump
where 1 { C o j u m p s ( i , j ) > 0 } is an indicator function equals 1 when a co-jump is identified between asset i and asset j at time s. We use the test statistics Φ n ( d ) in the JT co-jump test to identify co-jumps.
The results are reported in Table 3, Table 4, Table 5 and Table 6, where the significance level is 5%. Table 3 contains the empirical sizes of the co-jump tests when the underlying process is continuous (i.e., DGPs 1–2). At the 5% significance level, the range of the empirical size for the BLT test is (0.0687, 0.0718); for the JT test, it is (0.0143, 0.0220); for the co-exceedance rule of the LM test, it is (0.0426, 0.0435); and for the co-exceedance rule of the intersection of the LM and BNS tests, it is (0.0014,0.0010). Hence, the BLT test and co-exceedance rule of the JT test are close to the nominal sizes of all DGPs. The JT test is undersized and the co-exceedance rule of the LM–BNS test is far below the nominal size.
Table 4 summarizes the empirical sizes of the co-jump tests when there are only idiosyncratic jumps (i.e., DGPs 3–10). Adding idiosyncratic jumps significantly increases the size of the JT test. For example, the size increases from 0.0143 to 0.1771 at a 5% significance level in DGPs 2 and 7, and these results with other DGPs are quite similar. For the BLT test, the co-exceedance rule of the LM test and the co-exceedance rule of the LM–BNS test, adding idiosyncratic jumps, slightly lower its size. Similarly, the BLT test and co-exceedance rule of the JT test are close to the nominal size of all DGPs. The JT test is significantly oversized and the co-exceedance rule of the LM–BNS test is far below the nominal size.
Table 5 shows the empirical power of the co-jump tests when there are no idiosyncratic jumps. The JT test, the BLT test, and the co-exceedance rule of the LM test have significantly larger powers compared to the co-exceedance rule of the intersection of the LM–BNS test. As noted in the table, at the 5 % significance level, the range of the empirical power of the BLT test is (0.9167, 0.9362), and the JT test is (0.9451, 0.9891), while the co-exceedance rule of the LM test is (0.9362, 0.9722), and the intersection of the LM–BNS test is (0.1702, 0.3409). Jump intensity and jump size have positive effects on power in most cases. For example, in DGPs 11 and 13, as the co-jump size σ C J increases from 0.3 to 0.4, the empirical power of the BLT co-jump test increases from 0.9255 to 0.9362, the co-exceedance rule of the LM test increases from 0.9362 to 0.9521, the co-exceedance rule of the intersection of the LM–BNS test increases from 0.1702 to 0.2606, and the JT test increases from 0.9451 to 0.9670.
Table 6 summarizes the empirical power of the co-jump tests when there are also idiosyncratic jumps. Similarly, the JT co-jump test, the BLT test, and the co-exceedance rule of the LM test have significantly larger powers than the co-exceedance rule of the intersection of the LM and BNS tests. It is also worth noting that adding idiosyncratic jumps lowers the empirical power of the co-jump tests, particularly for the BLT test and the co-exceedance rule of the LM test. The jump size impacts power positively for all co-jump tests. The jump intensity does not impact the power unanimously. For example, the jump intensity has negative effects on power in the BLT test, LM test, and LM–BNS, but it has positive effects on the power of the JT test. We have different results on the effects of jump intensity because when jump intensity increases, it increases both the co-jump intensity and idiosyncratic jump intensity.
Overall, the BLT test, the JT test, and the co-exceedance rule of the LM test all have good testing power, except for the co-exceedance rule of the LM–BNS test. Additionally, the BLT test and the co-exceedance rule of the LM test have empirical sizes close to the nominal sizes in all DGPs, while the JT test is mixed, under- or over-sized in different DGPs, and the co-exceedance rule of the LM–BNS test is significantly undersized in all DGPs. The co-exceedance rule of the LM test has good testing power and is more conservative under the null. Hence, this test is very unlikely to yield the false discovery of co-jumps. Therefore, we used this test in the following.

4. Forecasting Methods

4.1. Forecasting Models

In this section, we summarize all models we used in predicting sector-level equity volatility. With the recent availability of high-frequency financial data, the HAR-type models are now widely used in the literature for forecasting the time-varying return variance. The HAR-RV model, first proposed by Corsi (2009), is based on the idea of the HARCH model and the heterogeneous market hypothesis developed by Müller et al. (1997), in which the authors believe that different types of traders react to the market on different horizons and, consequently, affect the market differently. So the HAR-RV model incorporates different kinds of market volatilities over different time horizons to forecast volatility. One simple version of the HAR-RV model is:
R V t + h = β 0 + β d R V t + β w R V t 5 , t + β m R V t 22 , t + ϵ t + h ,
where R V t + h is the realized volatility defined in Equation (3) at day t + h . The multi-period normalized realized measures for realized volatility R V t h , t are defined as the average of the corresponding one-period measures. Namely,
R V t h , t = h 1 [ R V t h + 1 + R V t h + 2 + + R V t ] .
Built on the HAR-RV model, Andersen et al. (2007) further separated jumps from the continuous part and developed a new HAR-RV-CJ prediction model that contained the jump component, the continuous component, and the lag terms over different time horizons. More specifically, the HAR-RV-CJ model is given by:
R V t , t + h = β 0 + β c d C t + β c w C t 5 , t + β c m C t 22 , t + β j d J t + β j w J t 5 , t + β j m J t 22 , t + ϵ t + h ,
where C t is the variation of the continuous component on day t estimated by BPV defined in Equation (4) and J t is the corresponding jump variation estimated by
J t = R V t B P V t
Since the sampling interval in the empirical analysis is finite ( Δ n > 0 ), the above-defined J t can be negative in some cases. Andersen et al. (2007) suggested excluding the negative values of the estimates and truncating the jump measurements at 0. The jump component variation is calculated as:
J t = max [ R V t B P V t , 0 ]
The regressors J t h , t and C t h , t are defined in the same way as R V t h , t above.
In the sequel, we estimate various HAR-RV-J models that incorporate co-jumps ( C J t ) and idiosyncratic jumps ( I J t ) and analyze the predictive content for co-jumps and idiosyncratic jumps in the forecasting volatility.
Specification 1: HAR-RV-C model (benchmark model):
log R V t + 1 = β 0 + β c d log C t + β c w log C t 5 , t + β c m log C t 22 , t + ϵ t + 1 .
Specification 2: HAR-RV-C-CJ model with co-jumps ( C J t ):
log R V t + 1 = β 0 + β c d log C t + β c w log C t 5 , t + β c m log C t 22 , t + β c j d log C J t + β c j w log C J t 5 , t + β c j m log C J t 22 , t + ϵ t + 1 .
Specification 3: HAR-RV-C-IJ model with idiosyncratic jumps ( I J t ):
log R V t + 1 = β 0 + β c d log C t + β c w log C t 5 , t + β c m log C t 22 , t + β i j d log I J t + β i j w log I J t 5 , t + β i j m log I J t 22 , t + ϵ t + 1 .
Specification 4: HAR-RV-C-CJ-IJ model with co-jumps ( C J t ) and idiosyncratic jumps ( I J t ):
log R V t + 1 = β 0 + β c d log C t + β c w log C t 5 , t + β c m log C t 22 , t + β c j d log C J t + β c j w log C J t 5 , t + β c j m log C J t 22 , t + β i j d log I J t + β i j w log I J t 5 , t + β i j m log I J t 22 , t + ϵ t + 1
The regressors include the daily, weekly, and monthly lag terms. All of the model specifications are summarized in Table 7.

4.2. Co-Jump and Idiosyncratic Jump Variations

The co-exceedance rule of the LM test is utilized to identify co-jumps and idiosyncratic jumps among M assets, as defined in (29). The co-jump variation is measured using the power transformation of the instantaneous return:6
C J t = s = 1 n r s 2 * I s C J
where I s C J is a co-jump indicator for the target asset at intraday interval s on day t and r s is the return over the same interval. The number of observations per day is n = 78 . The idiosyncratic jump variation is
I J t = s = 1 n r s 2 * I s I J
where I s I J is an idiosyncratic jump indicator for the target asset at the intraday interval s on day t.
The ratios of co-jump and idiosyncratic jump contributions to the total variations are calculated as:
r a t i o C J = C J t R V t and r a t i o I J = I J t R V t .
The proportions of co-jumps and idiosyncratic jumps to the total number of jumps are defined as:
P o r p C J = N C J N C J + N I J and P o r p I J = N I J N C J + N I J ,
where N C J is the number of co-jumps and N I J is the number of idiosyncratic jumps7. The frequencies of co-jumps and idiosyncratic jumps are defined as the percentage of days identified as having co-jumps and idiosyncratic jumps, respectively:
F r e q u e n c y C J = 1 T i = 1 T 1 { s = 1 M I s C J } > 0 and F r e q u e n c y I J = 1 T i = 1 T 1 { s = 1 M I s I J } > 0 ,
where T is the total number of days in the sample.

5. Setup of the Forecasting Experiment

For the forecasting experiment, we carried out both an in-sample regression analysis and an out-of-sample forecast analysis. We used a rolling window estimation scheme, with window size T = 450 (i.e., approximately two years). We obtained sequences of daily out-of-sample volatility forecasts for the sample period from 13 November 2008 to 31 December 2019, which contained 2801 trading days. The in-sample regression analysis was also conducted using the same sample periods.
We used in-sample adjusted R 2 , the mean absolute forecasting error (MAFE), and out-of-sample R 2 to evaluate forecasting results among candidate models.
The out-of-sample R 2 is defined as:
out-of-sample R 2 = 1 t = 1 T ( y t y t ^ ) 2 t = 1 T ( y t y t ¯ ) 2
where y t ¯ is the historical average of volatility.
Moreover, we implemented the model confidence set (MCS) test of Hansen et al. (2011) to compare the forecasting performance amongst candidate models, in which a large p-value indicates a superior performance. The MCS test has a sequence of statistic tests, which are utilized to determine the “Superior Set Models” (SSM)8 from a collection of models, where the null hypothesis of equal predictive ability among candidate models is not rejected at a certain confidence level. Let d j h , i denote the loss difference between the forecasting model j and h:
d j h , i = l ( ϵ i j ) l ( ϵ i h ) , j , h = 1 , , Q , i = 1 , , n ,
where ϵ i j is the prediction error in model j at time span i, ϵ i h is the prediction error in the other model h, and l ( · ) represents the quadratic loss function. The null hypothesis of the equal predictive ability for a given set of models Q can be given as:
H 0 , Q : c j , h = 0 , for all j , h = 1 , 2 , , Q ,
H 1 , Q : c j , h 0 , for all j , h = 1 , 2 , , Q ,
where c j , h = E ( d j , h ) is finite- and time-independent. The first test statistic is constructed as:
T j , h = d ¯ j , h v a r ^ ( d ¯ j , h )
where d ¯ j , h = n 1 i = 1 n d j h , i and v a r ^ ( d ¯ j , h ) is the bootstrapped estimate of v a r ( d ¯ j , h ) . We performed a stationary bootstrap procedure of 2000 resamples. In Hansen et al. (2011), the first test statistics maps naturally into the following test statistic:
T R , Q = max j , h Q T j , h
where T j , h is calculated using Equation (49). Since the asymptotic distribution of the test statistic is nonstandard, the relevant distribution can be estimated with bootstrap methods. We calculate the p-values for the candidate forecasting models. The MCS test used to obtain the SSM can be summarized as the following step:
  • Set Q = Q 0 .
  • Test for the null hypothesis of equal predictive ability among candidate models in Q with significance level α . If the null hypothesis cannot be rejected, then the set SSM is the same as the set of all models. Otherwise, determine the worst model in the set Q.
  • Remove the worst model from Q, and proceed to step 2.
The choice of the worst model to be eliminated in step 2 is determined using an elimination rule that is consistent with the test statistics defined in Equation (49), which are
e R , Q = arg max j sup j Q d ¯ j , h v a r ^ ( d ¯ j , h )

6. Data Description

We collected data for nine SPDR exchange-traded funds (ETFs), which were designed to track nine sectors of the S&P 500 index, including the material sector (XLB), energy sector (XLE), financial sector (XLF), industrial sector (XLI), technology sector (XLK), consumer staples sector (XLP), utility sector (XLU), healthcare sector (XLV), and consumer discretionary sector (XLY). The sources of nine sector ETF intraday observations were from the trade and quote database (TAQ) through the Wharton research data service (WRDS). Table 8 shows the details of the forecasting targets. The full sample period was from 3 January 2006 to 31 December 2019 for a total of 3523 trading days. We excluded overnight returns from our sample. The previous tick method derived from Gençay et al. (2001) was used to filter out price data. Raw data were cleaned and outliers were removed based on the procedures in Brownlees and Gallo (2006). We used the 5-min sampling frequency to reduce the effect of microstructure noise, yielding 78 observations per day in most cases.

7. Empirical Findings

In this section, we discuss the main empirical findings. We begin by discussing the results from identifying co-jumps and idiosyncratic jumps among nine sector ETFs in the S&P 500 market. The first sense of the prevalence of co-jumps and idiosyncratic jumps can be formulated by inspecting Table 9 and Figure 1, Figure 2 and Figure 3. Then we turn to the one-day-ahead predictive performance of the forecasting models outlined in Section 4. We report the model confidence set results (Table 10), in-sample adjusted R 2 results (Table 11), out-of-sample R 2 results (Table 12), and additionally, the mean absolute forecasting error (MAFE) results (Table 13). Our main empirical findings are summarized as follows.
First, both co-jumps and idiosyncratic jumps occurred frequently from January 2007 to December 2019, especially during significant market turmoil, such as the 2008 financial crisis and the 2011 European debt crisis. Table 9 shows the summary statistics for co-jumps and idiosyncratic jumps. While both types of jumps occurred frequently from January 2008 to December 2019, the frequency of co-jumps was much higher. Moreover, the descriptive statistics also suggested that the magnitude (including both the upside jump and downside jump) of the co-jumps was typically larger than that of the idiosyncratic jumps. All these results are in line with the fact that there were two major crises in the sample period, during which the systematic jump risk would be more prevalent than the idiosyncratic jump risk.
Second, co-jumps occurred more frequently than idiosyncratic jumps. Figure 1 shows the number of co-jumps under different scenarios, i.e., across a different number of sectors from January 2007 to December 2019. Not surprisingly, co-jumps detected under scenarios with fewer sectors occurred more frequently than co-jumps identified under scenarios with a larger number of sectors. Figure 2 and Figure 3 depict the numbers of co-jumps and idiosyncratic jumps, respectively. The results in Figure 2 and Figure 3 indicate that co-jumps and idiosyncratic jumps were more densely populated during significant market turmoils, including the 2008 financial crisis, the 2011 European debt crisis, the 2015–2016 stock market selloff, and the stock market dip in December 2018–January 2019. This suggested higher levels of systematic jump risk and idiosyncratic jump risk during those market downturns. Moreover, as expected, the financial sector had more idiosyncratic jumps during the 2008 financial crisis compared with other sectors, which is consistent with the fact that this crisis was mainly triggered by the collapse of the financial sector.
Third, co-jumps and idiosyncratic jumps demonstrate marginal predictive content for volatility prediction, and this result was robust across different forecasting targets. The results in the model confidence set test (Table 10), in-sample adjusted R 2 (Table 11), out-of-sample R 2 (Table 12), and mean absolute forecasting error (Table 13) indicate that the HAR-RV-C-CJ-IJ forecasting model, which contains co-jumps and idiosyncratic jumps variations, outperforms other forecasting models (HAR-RV-C, HAR-RV-C-CJ, and HAR-RV-C-IJ) in most scenarios. Additionally, when comparing the HAR-RV-C benchmark model with HAR-RV-C-CJ, it is clear that adding co-jump variations into the model reduces the forecasting error and improves in-sample adjusted R 2 and out-of-sample R 2 in most cases. This suggests both systematic risk and idiosyncratic risk are useful in predicting volatility.
Fourth, the predictive content of co-jumps is higher than that of idiosyncratic jumps. Evidently, the model with co-jump variations (HAR-RV-C-CJ) represents the majority of entries with smaller MAFEs and larger out-of-sample R 2 and in-sample R 2 than the model with idiosyncratic jump variations (HAR-RV-C-IJ). For example, in Table 13, the HAR-RV-C-CJ model yields the lowest MAFEs for all sectors, and these MAFEs are comparatively lower than models using idiosyncratic jump information9. The same finding holds with the model confidence set testing results in Table 10. Comparing the HAR-RV-C-IJ model with idiosyncratic jump variations, the HAR-RV-C-CJ model with co-jump variations has larger p-values in all sectors. In Table 11 and Table 12, the HAR-RV-C-CJ model has higher in-sample R 2 and out-of-sample R 2 compared with the HAR-RV-C-IJ model. All these results suggest that systematic risk is more useful than idiosyncratic risk in volatility forecasting.

8. Policy Implications

In this paper, we distinguished jumps in sector-level ETF prices into idiosyncratic jumps, which represent large-scale price movements that occur mostly within one sector, and co-jumps, which are defined as large-scale price movements that take place across more than one sector around the same time. Specifically, idiosyncratic jumps reflect the influence of events that are mostly confined within one sector, while co-jumps, by the nature of price co-movements among more than one sector at the same time, capture the impact of events that can spill over to other sectors, leading to financial risk contagion.
In our empirical experiment, the augmentation of co-jumps leads to improvements in the forecasting of realized volatility for each of the five S&P 500 sectors. Co-jumps also demonstrate higher predictive power compared to idiosyncratic jumps and achieve even better predictive accuracy when combined with idiosyncratic jumps. These results put co-jumps as a key source of volatility for S&P 500 sectors, aside from within-sector risks or events that are captured by idiosyncratic jumps. Spillovers from other sectors account for a significant portion of volatility for each individual market sector, which should be adequately noticed by researchers and policymakers. Moreover, we also found that co-jumps occurred two to three times more frequently than idiosyncratic jumps in our sample period from 2007 to 2019, suggesting within-sector risks or events are more likely to spill over to other sectors.
For policymakers, it is crucial to be vigilant to within-sector risks or events to prevent systemic risks to the financial market, given the outsized effects on other sectors and the frequency of such scenarios. Significant risks or events of one sector could spill over to other sectors and cause increased volatility, resulting in contagion among different sectors.

9. Concluding Remarks

In this paper, we examined the usefulness of co-jump and idiosyncratic jump variation measures in forecasting the volatility of sector ETFs. The co-jumps are meant to represent the systematic risks across multiple assets, sectors, and markets, while the idiosyncratic jumps capture the idiosyncratic risks in an individual sector. In particular, using 5-min return data, we utilized extant jump and co-jump tests to disentangle individual sector idiosyncratic jumps and multiple sector co-jumps from January 2007 to December 2019. Perhaps not surprisingly, we found that co-jumps and idiosyncratic jumps occurred more frequently during the financial market turmoils, such as the 2008 financial crisis and the 2011 European debt crisis, affirming the severity of financial stress during the market turmoil. Moreover, we found that co-jumps have larger magnitudes and variations compared with idiosyncratic jumps, which suggests that systematic risk has a larger impact than idiosyncratic risk on the financial market. Additionally, we showed (via a series of forecasting experiments that utilized heterogeneous autoregressive (HAR) models) that volatility forecasting accuracy can be improved when the original HAR forecasting models are augmented by including co-jump and idiosyncratic jump variation measures. We also find that the predictive content of co-jumps is higher than that of idiosyncratic jumps, suggesting that systematic risk is more useful than idiosyncratic risk in volatility forecasting.
In addition, we carried out Monte Carlo simulations to examine the relative performances of four existing co-jump tests, including the BLT co-jump test by Bollerslev et al. (2008), the JT co-jump test by Jacod and Todorov (2009), the co-exceedance rule by Gilder et al. (2014) using the LM test by Lee and Mykland (2007), and the co-exceedance rule using the intersection of the LM test and BNS tests by Barndorff-Nielsen and Shephard (2006). We find that the BLT test and the co-exceedance rule of the LM test outperform the JT co-jump test and co-exceedance rule of the intersection of the LM and BNS tests. The BLT test, the JT test, and the co-exceedance rule of the LM test all have good testing power, except for the co-exceedance rule of the LM–BNS test. The BLT test and the co-exceedance rule of the LM test have empirical sizes close to the nominal sizes in all DGPs, while the JT test is mixed, under- or over-sized in different DGPs, and the co-exceedance rule of the LM–BNS test is significantly undersized.

Author Contributions

Conceptualization, W.P. and C.Y.; methodology, W.P. and C.Y.; software, W.P. and C.Y.; validation, W.P. and C.Y.; formal analysis, W.P. and C.Y.; investigation, W.P. and C.Y.; resources, W.P. and C.Y.; data curation, W.P. and C.Y.; writing—original draft preparation, W.P. and C.Y.; writing—review and editing, W.P. and C.Y.; visualization, W.P. and C.Y.; supervision, W.P. and C.Y.; project administration, W.P. and C.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors are grateful to Rodger Klein, Yuan Liao, Norman Swanson, and Xiye Yang for the useful comments and suggestions on the topics explored in this paper. We also thank the participants of the Southwestern Finance Association Annual Meeting, the Midwest Economic Association Annual Meeting, and workshops at Rutgers University for the helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

Notes

1
2
Bandi and Russell (2008) show in the presence of market microstructure noise, realized variance does not identify daily integrated variance of the frictionless equilibrium. They propose a more general treatment of the effect of market microstructure noise on realized variance estimates.
3
We follow the parameter settings from Jacod and Todorov (2009).
4
Among 10 assets, co-jumps occur more frequently among a small portion of all assets (e.g., m = 2 , 3 , 4 ) than a large portion of all assets (e.g., m = 5 , 6 , 7 , 8 , 9 , 10 ) at each time.
5
Gilder et al. (2014) also implement the same approach.
6
One approach to measure jump variation is to use the difference between realized volatility and bipower variation, as shown in Equations (34) and (35). Another approach to examine the jump power variation is formed using power transformation of the instantaneous return, i.e., r s q . Key papers in this area include (Ding et al. 1993; Ding and Granger 1996; Todorov and Tauchen 2010; Barndorff-Nielsen et al. 2008, and the references cited there in).
7
The co-exceedance rule of the LM test is utilized to identify co-jumps and idiosyncratic jumps among assets.
8
The set of models that consists of the best model(s).
9
Overall, the mean absolute forecasting errors in Table 13 are comparatively smaller with respect to results in the relevant literature (e.g., see Qiu et al. 2019).

References

  1. Aït-Sahalia, Yacine, and Jean Jacod. 2009. Testing for jumps in a discretely observed process. The Annals of Statistics 37: 184–222. [Google Scholar] [CrossRef]
  2. Andersen, Torben G., and Jesper Lund. 1997. Estimating continuous-time stochastic volatility models of the short-term interest rate. Journal of Econometrics 77: 343–77. [Google Scholar] [CrossRef]
  3. Andersen, Torben G., Tim Bollerslev, and Francis X. Diebold. 2007. Roughing it up: Including jump components in the measurement, modeling, and forecasting of return volatility. The Review of Economics and Statistics 89: 701–20. [Google Scholar] [CrossRef]
  4. Andersen, Torben G., Tim Bollerslev, Francis X. Diebold, and Paul Labys. 2001. The distribution of realized exchange rate volatility. Journal of the American Statistical Association 96: 42–55. [Google Scholar] [CrossRef]
  5. Andersen, Torben G., Tim Bollerslev, Francis X. Diebold, and Paul Labys. 2003. Modeling and forecasting realized volatility. Econometrica 71: 579–625. [Google Scholar] [CrossRef] [Green Version]
  6. Andersen, Torben G., Tim Bollerslev, Per Frederiksen, and Morten Ørregaard Nielsen. 2010. Continuous-time models, realized volatilities, and testable distributional implications for daily stock returns. Journal of Applied Econometrics 25: 233–61. [Google Scholar] [CrossRef] [Green Version]
  7. Bandi, Federico M., and Jeffrey R. Russell. 2008. Microstructure noise, realized variance, and optimal sampling. The Review of Economic Studies 75: 339–69. [Google Scholar] [CrossRef]
  8. Bandi, Federico M., and Roberto Reno. 2016. Price and volatility co-jumps. Journal of Financial Economics 119: 107–46. [Google Scholar] [CrossRef]
  9. Barndorff-Nielsen, Ole E., and Neil Shephard. 2004. Power and bipower variation with stochastic volatility and jumps. Journal of Financial Econometrics 2: 1–37. [Google Scholar] [CrossRef] [Green Version]
  10. Barndorff-Nielsen, Ole E., and Neil Shephard. 2006. Econometrics of testing for jumps in financial economics using bipower variation. Journal of Financial Econometrics 4: 1–30. [Google Scholar] [CrossRef]
  11. Barndorff-Nielsen, Ole E., Silja Kinnebrock, and Neil Shephard. 2008. Measuring Downside Risk-Realised Semivariance. Available online: https://ssrn.com/abstract=1262194 (accessed on 29 June 2022).
  12. Bibinger, Markus, and Lars Winkelmann. 2015. Econometrics of co-jumps in high-frequency data with noise. Journal of Econometrics 184: 361–78. [Google Scholar] [CrossRef]
  13. Bollerslev, Tim, Tzuo Hann Law, and George Tauchen. 2008. Risk, jumps, and diversification. Journal of Econometrics 144: 234–56. [Google Scholar] [CrossRef] [Green Version]
  14. Boswijk, H. Peter, Roger J. A. Laeven, and Xiye Yang. 2018. Testing for self-excitation in jumps. Journal of Econometrics 203: 256–66. [Google Scholar] [CrossRef]
  15. Brandt, Michael W., and Christopher S. Jones. 2006. Volatility forecasting with range-based EGARCH models. Journal of Business & Economic Statistics 24: 470–86. [Google Scholar]
  16. Brownlees, Christian T., and Giampiero M. Gallo. 2006. Financial econometric analysis at ultra-high frequency: Data handling concerns. Computational Statistics & Data Analysis 51: 2232–45. [Google Scholar]
  17. Caporin, Massimiliano, Aleksey Kolokolov, and Roberto Renò. 2017. Systemic co-jumps. Journal of Financial Economics 126: 563–91. [Google Scholar] [CrossRef] [Green Version]
  18. Corsi, Fulvio. 2009. A simple approximate long-memory model of realized volatility. Journal of Financial Econometrics 7: 563–91. [Google Scholar] [CrossRef]
  19. Corsi, Fulvio, Davide Pirino, and Roberto Reno. 2010. Threshold bipower variation and the impact of jumps on volatility forecasting. Journal of Econometrics 159: 276–88. [Google Scholar] [CrossRef] [Green Version]
  20. Corradi, Valentina, Mervyn J. Silvapulle, and Norman R. Swanson. 2018. Testing for jumps and jump intensity path dependence. Journal of Econometrics 204: 248–67. [Google Scholar] [CrossRef]
  21. Ding, Zhuanxin, and Clive W. J. Granger. 1996. Modeling volatility persistence of speculative returns: A new approach. Journal of Econometrics 73: 185–215. [Google Scholar] [CrossRef]
  22. Ding, Zhuanxin, Clive W. J. Granger, and Robert F. Engle. 1993. A long memory property of stock market returns and a new model. Journal of Empirical Finance 1: 83–106. [Google Scholar] [CrossRef]
  23. Gençay, Ramazan, Michel Dacorogna, Ulrich A. Muller, Olivier Pictet, and Richard Olsen. 2001. An Introduction to High-Frequency Finance. Cambridge: Academic Press. [Google Scholar]
  24. Ghysels, Eric, and Arthur Sinko. 2011. Volatility forecasting and microstructure noise. Journal of Econometrics 160: 257–71. [Google Scholar] [CrossRef]
  25. Ghysels, Eric, Pedro Santa-Clara, and Rossen Valkanov. 2006. Predicting volatility: Getting the most out of return data sampled at different frequencies. Journal of Econometrics 131: 59–95. [Google Scholar] [CrossRef] [Green Version]
  26. Gilder, Dudley, Mark B. Shackleton, and Stephen J. Taylor. 2014. Cojumps in stock prices: Empirical evidence. Journal of Banking & Finance 40: 443–59. [Google Scholar]
  27. Hansen, Peter R., Asger Lunde, and James M. Nason. 2011. The model confidence set. Econometrica 79: 453–97. [Google Scholar] [CrossRef] [Green Version]
  28. Hansen, Peter R., and Asger Lunde. 2005. A forecast comparison of volatility models: Does anything beat a GARCH (1, 1)? Journal of Applied Econometrics 20: 873–89. [Google Scholar] [CrossRef] [Green Version]
  29. Harvey, Andrew, Esther Ruiz, and Neil Shephard. 1994. Multivariate stochastic variance models. The Review of Economic Studies 61: 247–64. [Google Scholar] [CrossRef] [Green Version]
  30. Hausman, Jerry A. 1978. Specification tests in econometrics. Econometrica: Journal of the Econometric Society 46: 1251–71. [Google Scholar] [CrossRef] [Green Version]
  31. Huang, Xin, and George Tauchen. 2005. The relative contribution of jumps to total price variance. Journal of Financial Econometrics 3: 456–99. [Google Scholar] [CrossRef]
  32. Jacod, Jean, and Viktor Todorov. 2009. Testing for common arrivals of jumps for discretely observed multidimensional processes. The Annals of Statistics 37: 1792–838. [Google Scholar] [CrossRef] [Green Version]
  33. Jiang, George J., and Roel CA Oomen. 2008. Testing for jumps when asset prices are observed with noise—A “swap variance” approach. Journal of Econometrics 144: 352–70. [Google Scholar] [CrossRef]
  34. Koopman, Siem Jan, Borus Jungbacker, and Eugenie Hol. 2005. Forecasting daily variability of the S&P 100 stock index using historical, realised and implied volatility measurements. Journal of Empirical Finance 12: 445–75. [Google Scholar]
  35. Lahaye, Jérôme, and Christopher Neely. 2020. The role of jumps in volatility spillovers in foreign exchange markets: Meteor shower and heat waves revisited. Journal of Business & Economic Statistics 38: 410–27. [Google Scholar]
  36. Lee, Suzanne S., and Per A. Mykland. 2007. Jumps in financial markets: A new nonparametric test and jump dynamics. The Review of Financial Studies 21: 2535–63. [Google Scholar] [CrossRef]
  37. Mancini, Cecilia. 2009. Non-parametric threshold estimation for models with stochastic diffusion coefficient and jumps. Scandinavian Journal of Statistics 36: 270–96. [Google Scholar] [CrossRef]
  38. Mancini, Cecilia, and Fabio Gobbi. 2012. Identifying the brownian covariation from the co-jumps given discrete observations. Econometric Theory 28: 249–73. [Google Scholar] [CrossRef] [Green Version]
  39. Mukherjee, Arpita, Weijia Peng, Norman R. Swanson, and Xiye Yang. 2020. Financial econometrics and big data: A survey of volatility estimators and tests for the presence of jumps and co-jumps. In Handbook of Statistics. Amsterdam: Elsevier, vol. 42, pp. 3–59. [Google Scholar]
  40. Müller, Ulrich A., Michel M. Dacorogna, Rakhal D. Davé, Richard B. Olsen, Olivier V. Pictet, and Jacob E. Von Weizsäcker. 1997. Volatilities of different time resolutions—Analyzing the dynamics of market components. Journal of Empirical Finance 4: 213–39. [Google Scholar] [CrossRef]
  41. Patton, Andrew J., and Kevin Sheppard. 2015. Good volatility, bad volatility: Signed jumps and the persistence of volatility. Review of Economics and Statistics 97: 683–97. [Google Scholar] [CrossRef] [Green Version]
  42. Podolskij, Mark, and Daniel Ziggel. 2010. New tests for jumps in semimartingale models. Statistical Inference for Stochastic Processes 13: 15–41. [Google Scholar] [CrossRef]
  43. Qiu, Yue, Xinyu Zhang, Tian Xie, and Shangwei Zhao. 2019. Versatile HAR model for realized volatility: A least square model averaging perspective. Journal of Management Science and Engineering 4: 55–73. [Google Scholar] [CrossRef]
  44. Todorov, Viktor, and George Tauchen. 2010. Activity signature functions for high-frequency data analysis. Journal of Econometrics 154: 125–38. [Google Scholar] [CrossRef]
Figure 1. Number of co-jumps under different scenarios: Figure 1 shows the number of co-jumps identified under different scenarios, i.e., across nine sectors in the S&P 500 market from January 2007 to December 2019. The nine sectors are materials, energy, financial, industrial, technology, consumer staples, utility, healthcare, and consumer discretionary. The co-exceedance rule of the LM test is utilized to detect co-jumps among sectors. The shaded area in the figure denotes the financial crisis period from December 2007 to June 2009. See Section 4 and Section 7 for further details.
Figure 1. Number of co-jumps under different scenarios: Figure 1 shows the number of co-jumps identified under different scenarios, i.e., across nine sectors in the S&P 500 market from January 2007 to December 2019. The nine sectors are materials, energy, financial, industrial, technology, consumer staples, utility, healthcare, and consumer discretionary. The co-exceedance rule of the LM test is utilized to detect co-jumps among sectors. The shaded area in the figure denotes the financial crisis period from December 2007 to June 2009. See Section 4 and Section 7 for further details.
Jrfm 15 00334 g001
Figure 2. Number of sector-level idiosyncratic jumps: see notes in Figure 1. Figure 2 shows the number of idiosyncratic jumps in nine sectors in the S&P 500 market from January 2007 to December 2019. See Section 4 and Section 7 for further details.
Figure 2. Number of sector-level idiosyncratic jumps: see notes in Figure 1. Figure 2 shows the number of idiosyncratic jumps in nine sectors in the S&P 500 market from January 2007 to December 2019. See Section 4 and Section 7 for further details.
Jrfm 15 00334 g002
Figure 3. Number of sector-level co-jumps: see notes for Figure 1. Figure 3 shows the number of co-jumps in nine sectors in the S&P 500 market from January 2007 to December 2019. See Section 4 and Section 7 for further details.
Figure 3. Number of sector-level co-jumps: see notes for Figure 1. Figure 3 shows the number of co-jumps in nine sectors in the S&P 500 market from January 2007 to December 2019. See Section 4 and Section 7 for further details.
Jrfm 15 00334 g003
Table 1. Monte Carlo experiments—parameter settings *.
Table 1. Monte Carlo experiments—parameter settings *.
Jump Intensity λ = 0.1 , 0.2
Leverage Effect ρ = 0 , 0.5
Jump Distribution J t i . i . d N ( 0 , σ ) , where σ = 0.3 , 0.4
Sampling Frequency Δ n = 1 78
Other Parameters μ , κ v , θ v , ζ , N = 0.05 , 5 , 0.0144 , 0.5 , 10
* Notes: Entries in this table show the parameter settings in Equations (26) and (27) in the Monte Carlo experiments. N is the number of simulated log prices of assets in the data generating process(es) (DGP). Moreover, 1 78 indicates the 5-min sampling frequency. See Section 3 for complete details.
Table 2. Monte Carlo experiments—data generating process *.
Table 2. Monte Carlo experiments—data generating process *.
λ IJ λ CJ σ IJ σ CJ ρ Δ n μ κ v θ v ζ N
DGP1 0 1 78 0.0550.01440.510
DGP2 −0.5 1 78 0.0550.01440.510
DGP30.1 0.3 0 1 78 0.0550.01440.510
DGP40.2 0.3 0 1 78 0.0550.01440.510
DGP50.1 0.4 0 1 78 0.0550.01440.510
DGP60.2 0.4 0 1 78 0.0550.01440.510
DGP70.1 0.3 −0.5 1 78 0.0550.01440.510
DGP80.2 0.3 −0.5 1 78 0.0550.01440.510
DGP90.1 0.4 −0.5 1 78 0.0550.01440.510
DGP100.2 0.4 −0.5 1 78 0.0550.01440.510
DGP11 0.1 0.30 1 78 0.0550.01440.510
DGP12 0.2 0.30 1 78 0.0550.01440.510
DGP13 0.1 0.40 1 78 0.0550.01440.510
DGP14 0.2 0.40 1 78 0.0550.01440.510
DGP15 0.1 0.3−0.5 1 78 0.0550.01440.510
DGP16 0.2 0.3−0.5 1 78 0.0550.01440.510
DGP17 0.1 0.4−0.5 1 78 0.0550.01440.510
DGP18 0.2 0.4−0.5 1 78 0.0550.01440.510
DGP190.10.10.30.30 1 78 0.0550.01440.510
DGP200.20.20.30.30 1 78 0.0550.01440.510
DGP210.10.10.40.40 1 78 0.0550.01440.510
DGP220.20.20.40.40 1 78 0.0550.01440.510
DGP230.10.10.30.3−0.5 1 78 0.0550.01440.510
DGP240.20.20.30.3−0.5 1 78 0.0550.01440.510
DGP250.10.10.40.4−0.5 1 78 0.0550.01440.510
DGP260.20.20.40.4−0.5 1 78 0.0550.01440.510
* Notes: DGPs 1–2 are continuous processes. DGPs 3–10 are continuous processes plus idiosyncratic jumps characterized by various jump size densities. DGPs 11–18 are continuous processes plus co-jumps characterized by various jump size densities. DGPs 19–26 are continuous processes plus co-jumps and idiosyncratic jumps characterized by various jump size densities. λ C J is the co-jump intensity and λ I J is the idiosyncratic jump intensity. ρ is the leverage effect. See Section 3 for complete details.
Table 3. Empirical sizes of co-jump tests under continuous processes *.
Table 3. Empirical sizes of co-jump tests under continuous processes *.
CaseBLTJTLMLM–BNS
DGP10.06870.02200.04260.0014
DGP20.07180.01430.04350.0010
* Notes: Entries are rejection frequencies based on applications of the BLT co-jump test, JT co-jump tests, the co-exceedance rule of the LM test, and the co-exceedance rule of the LM–BNS test. Results for the 0.05 nominal size are reported. There are 1000 Monte Carlo replications. See Section 2.2 and Section 3 for complete details.
Table 4. Empirical sizes of co-jump tests under Monte Carlo simulations with idiosyncratic jumps *.
Table 4. Empirical sizes of co-jump tests under Monte Carlo simulations with idiosyncratic jumps *.
CaseBLTJTLMLM–BNS
DGP30.05650.18920.04780.0015
DGP40.05060.37510.05200.0022
DGP50.05410.24860.04800.0017
DGP60.04760.45100.05230.0028
DGP70.05540.17710.04660.0019
DGP80.04800.39820.05120.0027
DGP90.05260.21780.04680.0022
DGP100.04520.45430.05150.0034
* Notes: see notes in Table 3.
Table 5. Empirical power of co-jump tests under Monte Carlo simulations with co-jumps *.
Table 5. Empirical power of co-jump tests under Monte Carlo simulations with co-jumps *.
CaseBLTJTLMLM–BNS
DGP110.92550.94510.93620.1702
DGP120.91670.98360.94950.2121
DGP130.93620.96700.95210.2606
DGP140.92420.98910.96970.3409
DGP150.92020.94510.94150.2340
DGP160.92420.98360.96720.2020
DGP170.93090.95600.96280.3085
DGP180.92930.98910.97220.3409
* Notes: see notes in Table 3.
Table 6. Empirical power of co-jump tests under Monte Carlo simulations with co-jumps and idiosyncratic jumps *.
Table 6. Empirical power of co-jump tests under Monte Carlo simulations with co-jumps and idiosyncratic jumps *.
CaseBLTJTLMLM–BNS
DGP190.91260.93940.90290.2087
DGP200.84560.97980.93320.2166
DGP210.91750.96970.91260.3155
DGP220.86870.97980.95160.3433
DGP230.89320.95960.92230.2524
DGP240.86870.97470.94240.2097
DGP250.91260.96970.95150.3350
DGP260.88250.97980.95160.3295
* Notes: see notes in Table 3.
Table 7. Brief descriptions of forecasting models *.
Table 7. Brief descriptions of forecasting models *.
ModelDescription
HAR-RV-CFuture RV depends on the lags of the continuous component of the RV.
HAR-RV-C-CJFuture RV depends on the lags of the continuous component of RV and the co-jump variation.
HAR-RV-C-IJFuture RV depends on the lags of the continuous component of the RV and the idiosyncratic jump variation.
HAR-RV-C-CJ-IJFuture RV depends on the lags of the continuous component of RV, co-jump variation, and idiosyncratic jump variation.
* Notes: This table reports the models in the forecasting experiments. The benchmark model is the HAR-RV-C model. For complete details, refer to Section 4.
Table 8. Target variables in the forecasting experiment *.
Table 8. Target variables in the forecasting experiment *.
Target NameDescriptionTransformationFrequency
XLFFinancial Sector SPDR Fund l n ( X t ) l n ( X t 1 ) Daily
XLKTechnology Sector SPDR Fund l n ( X t ) l n ( X t 1 ) Daily
XLPConsumer Staples Select Sector SPDR Fund l n ( X t ) l n ( X t 1 ) Daily
XLUUtilities Select Sector SPDR Fund l n ( X t ) l n ( X t 1 ) Daily
XLYConsumer Discretionary Select Sector SPDR Fund l n ( X t ) l n ( X t 1 ) Daily
* Notes: This table reports the prediction targets. Data transformations used in the forecasting experiments are presented in the third column of the table. See Section 6 for further details.
Table 9. Descriptive statistics for co-jumps and idiosyncratic jumps (2007:1–2019:12) *.
Table 9. Descriptive statistics for co-jumps and idiosyncratic jumps (2007:1–2019:12) *.
SectorXLBXLEXLFXLIXLKXLPXLUXLVXLY
Co-jumpFrequency (%)30.10%24.45%25.31%30.78%24.48%32.73%29.55%29.49%29.03%
Prop. (%)79.36%70.08%72.79%87.46%83.52%79.76%64.37%78.23%86.09%
upside jump mean (%)0.63%0.78%0.72%0.56%0.57%0.43%0.61%0.49%0.55%
downside jump mean (%)−0.61%−0.76%−0.70%−0.54%−0.53%−0.43%−0.60%−0.47%−0.54%
St. dev.0.00740.00900.00900.00640.00640.00520.00720.00550.0067
Ratio of JV to RV (%)27.96%22.71%23.50%28.58%22.74%30.40%27.45%27.39%26.97%
Idiosyncratic jumpFrequency (%)15.01%17.42%14.82%10.39%10.36%15.19%24.05%16.56%10.36%
Prop. (%)20.64%29.92%27.21%12.54%16.48%20.24%35.63%21.77%13.91%
upside jump mean (%)0.50%0.59%0.67%0.42%0.38%0.31%0.43%0.39%0.40%
downside jump mean (%)−0.51%−0.59%−0.62%−0.43%−0.39%−0.31%−0.45%−0.40%−0.38%
St. dev.0.00560.00640.00770.00480.00420.00340.00470.00430.0043
Ratio of JV to RV (%)13.94%16.18%13.77%9.65%9.62%14.11%22.34%15.38%9.62%
* Notes: This table shows the descriptive statistics for co-jumps and idiosyncratic jumps identified from nine sector ETFs. The frequency denotes the percentage of days having co-jumps or idiosyncratic jumps. JV is the jump variation and RV is the total variation calculated by realized volatility. XLB is the material select sector ETF; XLE is the energy select sector ETF; XLF is the financial select sector EFT; XLI is the industrial select sector EFT; XLK is the technology select sector EFT; XLP is the consumer staples select sector EFT; XLU is the utility sector ETF; XLV is the healthcare sector ETF, and XLY is the consumer discretionary sector ETF. See Section 5 and Section 7 for further details.
Table 10. Model confidence set test *.
Table 10. Model confidence set test *.
SectorHAR-RV-CHAR-RV-C-CJHAR-RV-C-IJHAR-RV-C-CJ-IJ
Forecasting Period: 2008:11–2019:12
XLF0.6560 *1.0000 *0.6560 *0.6560 *
XLK0.2165 *0.4195 *0.1415 *1.0000 *
XLP0.1220 *0.1220 *0.08801.0000 *
XLU0.3860 *0.3860 *0.2975 *1.0000 *
XLY0.7780 *0.9385 *0.7780 *1.0000 *
* Notes: See notes in Table 9. This table reports the model confidence test (MCS) p-values associated with the MCS test discussed in Section 5. Moreover, the 90% model confidence set members are identified by models with asterisks.
Table 11. In-sample adjusted R 2 *.
Table 11. In-sample adjusted R 2 *.
SectorHAR-RV-CHAR-RV-C-CJHAR-RV-C-IJHAR-RV-C-CJ-IJ
Forecasting Period: 2008:11–2019:12
XLF0.77610.77650.77610.7768
XLK0.67480.67580.67580.6763
XLP0.64300.64410.64390.6443
XLU0.61570.61890.61750.6200
XLY0.74580.74700.74620.7470
* Notes: See notes in Table 10. Entries in this table include in-sample adjusted R 2 . All results are from forecasting experiments that utilized 5-min frequency data. See Section 5 and Section 7 for further details.
Table 12. Out-of-sample R 2 *.
Table 12. Out-of-sample R 2 *.
SectorHAR-RV-CHAR-RV-C-CJHAR-RV-C-IJHAR-RV-C-CJ-IJ
Forecasting Period: 2008:11–2019:12
XLF0.67650.68660.68230.6823
XLK0.56800.57720.56490.5817
XLP0.55670.56020.55560.5696
XLU0.49270.49700.49140.5047
XLY0.68040.68380.68020.6841
* Notes: See notes in Table 10. This table reports the out-of-sample R 2 for all forecasting experiments.
Table 13. Mean absolute forecasting error *.
Table 13. Mean absolute forecasting error *.
SectorHAR-RV-CHAR-RV-C-CJHAR-RV-C-IJHAR-RV-C-CJ-IJ
Forecasting Period: 2008:11–2019:12
XLF0.22930.22510.22760.2267
XLK0.23710.23480.23600.2320
XLP0.20620.20560.20630.2031
XLU0.20700.20610.20690.2042
XLY0.22250.22110.22140.2205
* Notes: See notes in Table 10. This table reports the mean absolute forecasting errors (MAFEs) for all forecasting experiments.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peng, W.; Yao, C. Co-Jumps, Co-Jump Tests, and Volatility Forecasting: Monte Carlo and Empirical Evidence. J. Risk Financial Manag. 2022, 15, 334. https://doi.org/10.3390/jrfm15080334

AMA Style

Peng W, Yao C. Co-Jumps, Co-Jump Tests, and Volatility Forecasting: Monte Carlo and Empirical Evidence. Journal of Risk and Financial Management. 2022; 15(8):334. https://doi.org/10.3390/jrfm15080334

Chicago/Turabian Style

Peng, Weijia, and Chun Yao. 2022. "Co-Jumps, Co-Jump Tests, and Volatility Forecasting: Monte Carlo and Empirical Evidence" Journal of Risk and Financial Management 15, no. 8: 334. https://doi.org/10.3390/jrfm15080334

APA Style

Peng, W., & Yao, C. (2022). Co-Jumps, Co-Jump Tests, and Volatility Forecasting: Monte Carlo and Empirical Evidence. Journal of Risk and Financial Management, 15(8), 334. https://doi.org/10.3390/jrfm15080334

Article Metrics

Back to TopTop