Next Article in Journal
Revisiting the Performance of MACD and RSI Oscillators
Previous Article in Journal
A Non-Parametric and Entropy Based Analysis of the Relationship between the VIX and S&P 500
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Testing for a Single-Factor Stochastic Volatility in Bivariate Series

1
Faculty of Engineering, Fukui University of Technology, 3-6-1 Gakuen, Fukui 910-8505, Japan
2
Faculty of Economics, Yokohama National University, 79-4 Tokiwadai, Yokohama 240-8501, Japan
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2013, 6(1), 31-61; https://doi.org/10.3390/jrfm6010031
Submission received: 4 October 2013 / Revised: 25 November 2013 / Accepted: 12 December 2013 / Published: 19 December 2013

Abstract

:
This paper proposes the Lagrange multiplier test for the null hypothesis that the bivariate time series has only a single common stochastic volatility factor and no idiosyncratic volatility factor. The test statistic is derived by representing the model in a linear state-space form under the assumption that the log of squared measurement error is normally distributed. The empirical size and power of the test are examined in Monte Carlo experiments. We apply the test to the Asian stock market indices.
Classification: PACS JEL classification:
C12; C32; C58

1. Introduction

The international co-movement of stock market volatility has attracted much attention after the Asian financial crisis in the late 1990s. It is often argued that national financial instability is quickly transmitted internationally through accelerated capital flows, especially in emerging markets and then international financial crisis is more difficult to cope with.
One type of co-movement analysis is to estimate volatility spillover effects, which is essentially transmission mechanism of volatility in stock markets. This line of research uses extensively various variants of ARCH model first proposed by Engle [1], such as the SWARCH model, which is the Markov switching model incorporated into the ARCH model. The existence of spillover effects has been recognized almost universally by empirical analyses.
Another type of co-movement analysis is to identify a volatility factor common to multiple series. This model is often used in explaining volatility contagion, namely increased correlation in high volatility period; correlation between two series increases when the volatility of their common factor is high in comparison with idiosyncratic volatilities. This phenomenon is also regarded as a stylized fact in financial markets.
In this paper we proposes the Lagrange multiplier (LM) test for a single common volatility factor model. The null hypothesis is an extreme case in that two markets have a single stochastic volatility process in common and have no idiosyncratic volatility factor. This situation is not unrealistic in that it expresses the situation where a single intermarket volatility factor overrides other minor volatility factors of two markets, for example, in the period that includes global financial crises. In such an “extreme" situation, other minor domestic volatility factors are negligible in size and markets looked as if they have a single common volatility process in common. The aim of our test is not to assert that there exists only one single volatility factor among markets, but that minor idiosyncratic volatility factors are overridden by a single intermarket volatility factors.
The use of the LM test is essential here. The Wald and Likelihood ratio test statistics cannot have χ 2 distribution asymptotically under the null hypothesis; the restricted parameter estimators, which is used by the Wald and Likelihood ratio test procedures, cannot have asymptotically standard distribution, since the null for our problem is on the boundary of the parameter space (General discussions of asymptotic tests in boundary situations have been provided by Chernoff [2], Moran [3], and Chant [4].), namely variance should be positive and correlation should be less than one in absolute value. Only the LM test statistic has χ 2 distribution asymptotically under the null hypothesis of our problem, since the LM test procedure uses only restricted parameter estimators.
We here consider the linearized form of the bivariate Stochastic volatility (SV) model and assume that this linearized form has a Gaussian disturbance term. This unconventional assumption is a cost to derive the test statistic for the number of volatility factors, which is an unprecedented challenge in the literature. By this assumption, the model is free from numerical integration to be used in estimating the original SV model (The original formulation of the SV model was introduced by Taylor [5], and has been applied successfully in financial analyses, for example, by Jacquier et al. [6] and Kim et al. [7]. For the recent development of the multivariate SV models, see Asai et al. [8].). A counterpart test in the ARCH framework was proposed by Engle and Susmel [9]. We believed that our test would be also useful to find a single overriding factor in bivariate state-space model in level.
The rest of the paper is organized as follows. Section 2 details the framework and notation of the model. Section 3 proposes the formula of the test statistic for the linear state-space model. Section 4 reports the results of our Monte Carlo experiments for the empirical size and power of the test. In Section 5, the test is applied to the empirical analysis of the Asian stock market indexes. Section 6 briefly concludes this study. The derivation of the score functions and the conditional moments of state variables are given in Appendix.

2. Model

We here consider the following observation equations of the bivariate SV model:
r 1 t = ψ 1 h 1 t 2 e 1 t , r 2 t = ψ 2 h 2 t 2 e 2 t
for series k = 1 , 2 and time t = 1 , , T . Note that r k t is the observed variable, which typically denotes the demeaned rate of return on assets, h k t is demeaned part of log volatility, where E ( h k t ) = 0 , e k t is an idiosyncratic noise with constant variance, standardized at 1, and log ψ k 2 is the mean of log volatility, and hence ψ k is, roughly speaking, overall standard deviation of r k t .
We linearize the observation equation to the linear state-space model (The log transformation breaks down when r k t = 0 , since log ( 0 ) = - . In practice, the use of log ( x + ϵ ) - ϵ / ( ϵ + x ) has been suggested in place of log x for small constant ϵ > 0 to avoid this problem. This transformation is useful in that ϵ has no effect for large x. See Breidt and Carriquiry [10] for details.); taking log of squares of Equation (1), we have
y 1 t = δ 1 + h 1 t + ξ 1 t , y 2 t = δ 2 + h 2 t + ξ 2 t
where y k t log r k t 2 , δ k log ψ k 2 + μ , ξ k t log e k t 2 - μ , and μ E log e k t 2 .
The main difference of our model from the original multivariate SV model is that Gaussianity of ( ξ 1 , ξ 2 ) is assumed in our model, namely
ξ 1 t ξ 2 t i.i.d. N 0 0 , η η γ η γ η
where γ is a correlation coefficient of ( ξ 1 , ξ 2 ) , and hence the distribution of e k t is no longer Gaussian. We are not asserting that the log squared returns follow a Gaussian distribution in reality (As suggested by an anonymous reviewer, the distribution of ( e 1 t , e 2 t ) under Gaussianity of ( ξ 1 t , ξ 2 t ) is, exactly speaking, inappropriate for market data. It can be easily shown that the marginal distribution of e k t is expressed as
f ( e k t ) = 1 2 π η e k t exp - ( log e k t 2 - μ ) 2 2 η , - < e k t <
noting that e k t 2 = exp ( μ + ξ k t ) follows a log-normal distribution. It is easy to see that this distribution has zero density function at the origin.).
Under the assumption of Gaussianity of ( e 1 t , e 2 t ) , the model (1) would be a nonlinear non-Gaussian multivariate state space model; even conventional hypothesis testing for this model has been rarely calculated because of its heavy computational burden and complexity, and hence constructing a test for the number of factors, which is our aim, seems almost unfeasible under the original assumption.
It seems also algebraically unfeasible to derive a test statistics from Equation (2) if Gaussianity of ( e 1 t , e 2 t ) is assumed. The properties of the linearized model under the Gaussianity assumption of ( e 1 t , e 2 t ) was extensively studied by Harvey et al. [11] and they reported that their analytical expression is extremely complicated; for example, the joint covariance of log e 1 t 2 and log e 2 t 2 is expressed as an infinite series; the mean and variance of log e k t 2 are expressed as μ = - 1.27 and η = π 2 / 2 4.93 , respectively.
On the other hand, the linearized model (2) is a linear state-space model under our assumption that the distribution of ( ξ 1 t , ξ 2 t ) is Gaussian and we can derive the score function and Fisher information, which are required to derive our test, though algebraically tedious, as shown in Appendix A.
The non-Gaussianity of ( e 1 t , e 2 t ) in our model is a serious drawback if we use this model for asset pricing, where Gaussianity of return is an essential assumption. However, in the context of spillover or contagion in international finance, Gaussianity is not always an essential element. We believe that our assumption of Gaussian ( ξ 1 t , ξ 2 t ) is a cost to derive the test for a single volatility factor. The usefulness of our assumption can be tested by checking the usefulness of our test.
The observation equation density, namely the conditional density of y 1 t and y 2 t given h 1 t and h 2 t , are expressed as
ν 1 t f ( y 1 t | h 1 t ) = 1 2 π η exp - ( y 1 t - δ 1 - h 1 t ) 2 2 η , ν 2 t f ( y 2 t | y 1 t , h 1 t , h 2 t ) = 1 2 π η ( 1 - γ 2 ) exp - [ ( y 2 t - δ 2 - h 2 t ) - γ ( y 1 t - δ 1 - h 1 t ) ] 2 2 η ( 1 - γ 2 ) .
We here assume that the log volatilities follow the bivariate first-order autoregressive (AR) process:
h 1 t = ϕ 1 h 1 , t - 1 + ω 1 1 / 2 u 1 , t - 1 , ϕ 1 < 1 h 2 t = ϕ 2 h 2 , t - 1 + λ ω 1 1 / 2 u 1 , t - 1 + ω 2 1 / 2 u 2 , t - 1 , ϕ 2 < 1
where the idiosyncratic noise u t = ( u 1 t , u 2 t ) follows Gaussian distribution, namely
u 1 t u 2 t i.i.d. N 0 0 , 1 0 0 1 .
The transition Equation (5), or the conditional densities of h 1 t and h 2 t given their past values, are expressed as
τ 1 t f ( h 1 t | h 1 , t - 1 ) = 1 2 π ω 1 exp - ( h 1 t - ϕ 1 h 1 , t - 1 ) 2 2 ω 1 , τ 2 t f ( h 2 t | h 1 t , h 1 , t - 1 , h 2 , t - 1 ) = 1 2 π ω 2 exp - [ ( h 2 t - ϕ 2 h 2 , t - 1 ) - λ ( h 1 t - ϕ 1 h 1 , t - 1 ) ] 2 2 ω 2 .
This formulation is not a general bivariate autoregressive process in that they are two univariate first-order AR processes with correlated innovations (Harvey et al. [11] set ϕ 1 , ϕ 2 to unity.). The general condition for the identity h 1 t = h 2 t is that the coefficient matrix of the multivariate first-order AR model has eigenvector ( 1 , 1 ) , in addition to the degeneracy of disturbances. However, we believe that our assumption is sufficient for practical use.
Then, the density function of y ˜ T ( y 11 , y 21 , , y 1 T , y 2 T ) , or the likelihood function, from which our LM test statistic is derived, is expressed as follows:
f ( y ˜ T ) = - - f ( y ˜ T , h ˜ T ) d h 11 d h 21 d h 1 T d h 2 T
by integrating out the state variables h ˜ T ( h 11 , h 21 , , h 1 T , h 2 T ) , where
f ( y ˜ T , h ˜ T ) = t = 1 T k = 1 2 ν k t τ k t .
We have two notes for the integration in Equation (8). First, our model defined by Equations (2), (3), (5), and (6) is linear and Gaussian, and hence the likelihood function can be evaluated and maximized, without numerical integration, whether by means of quadrature or Monte Carlo methods. Hence, we evaluate the likelihood by means of the Kalman filter method, which was suggested by Nelson [12] and Harvey et al. [11]. Second, the transition density is degenerate under the null hypothesis and the integration can be performed by integration by parts, as shown in Appendix A.
In the conventional multivariate SV model, which is different from ours, the observation equation disturbance term e t = ( e 1 t , e 2 t ) follows Gaussian distribution:
e 1 t e 2 t i.i.d. N 0 0 , 1 ρ ρ 1 .
Harvey et al. [11] showed that the maximum likelihood (ML) estimator of this transformed model, which is called the quasi-maximum likelihood estimator (QMLE), is consistent under the original assumption of the Gaussianity of e t . The estimation of SV model in the original form (9) can be carried out by various methods, for example, the generalized method of moments (GMM) used by Melino and Turnbull [13], the efficient method of moments (EMM) applied by Gallant et al. [14], and Markov Chain Monte Carlo (MCMC) procedures used by Jacquier et al. [6] and Kim et al. [7]. The construction of tests for the original model (9) based upon these methods is left for further research. However, little attention has been paid to hypothesis testing of this model.

3. Test Statistic

We propose the LM test for the hypothesis that the observation series y 1 t and y 2 t have the only stochastic volatility factor in common, and other idiosyncratic volatility factors are negligible, namely h 1 t = h 2 t for any t in Equation (5), against the alternative hypothesis that each series has different, possibly correlated, log volatility series h 1 t and h 2 t .
The null hypothesis that h 1 t = h 2 t is a first order approximation to the international synchronization of asset price fluctuation, found, or believed to be found, in the period that includes financial crises. In normal periods actual volatilities of markets are expressed as mixtures of many factors. In the period that includes financial crises, however, only a single intermarket volatility factor may override other minor volatility factors in markets. In such an extreme situation, other minor domestic volatility factors are negligible in size and some markets look as if they have a single common volatility factor. The aim of our test is not to assert there exists only a single volatility factor among markets, but is to check whether minor volatility factors are negligible in comparison with a big intermarket volatility factor.
In our Equation (7), the innovations are perfectly correlated if ω 2 = 0 , and mutually independent if ω 2 = . Our test checks that the two series have the same autoregressive coefficient and innovation variance, namely ϕ 1 = ϕ 2 , λ = 1 , as well as the degeneracy of the innovation, namely ω 2 = 0 . Therefore, in Equation (5), the null hypothesis is formalized as
H 0 : ϕ 1 = ϕ 2 , λ = 1 , ω 2 = 0 .
In our framework, in order to define the existence of a single overriding volatility factor in two markets, the volatility series of two markets should be exactly the same, which means ϕ 1 = ϕ 2 , λ = 1 , ω 2 = 0 , and hence h 1 t = h 2 t , when the volatility series is stationary; the logic is as follows: If the volatility processes of two market are not identical but correlated by assuming ω 2 0 , they have different volatility shocks and hence volatility series have different patterns especially when the absolute value of ω 2 is large, then we have no choice than to define the number of volatility factors in the market is two when ω 2 0 . Then the condition ω 2 = 0 , or very small absolute value of ω 2 , is a essential element in defining the single overriding volatility series. When ω 2 = 0 , we need the condition ϕ 1 = ϕ 2 , to define the single overriding volatility series, since ϕ 1 ϕ 2 implies different cyclicality or periodicity of volatility process.
It is arguable whether the condition λ = 1 is necessary. Under the assumption of ω 2 = 0 , the condition 0 < λ < 1 expresses smaller variance of h 2 t than that of h 1 t around 0 and hence, intuitively, smaller variance of exp ( h 2 t ) than that of exp ( h 1 t ) around 1, namely variance of volatility of r 1 t , disregarding proportionality constant exp ( δ 1 ) . Then λ can be regarded as “volatility of volatility” parameter and hence should be 1 if two markets follow the same volatility factor. If the same volatility process implies identical volatility of volatility, we need the condition λ = 1 .
We need the initial value condition h 10 = h 20 additionally for the log volatilities to be exactly identical stochastic processes under the conditions of (10). However, the log volatilities converge to zero quickly as long as they are stationary, even if their initial values are different, and the effect of ( h 10 , h 20 ) disappears in the vector moving average expression of h k t , as t increases. Then, we only consider the conditions (10) hereafter.
The Wald and Likelihood ratio principles are irrelevant to test for Equation (10), because the maximum likelihood estimator of ω 2 can have a singular distribution under the null hypothesis; the parameter value ω 2 = 0 under the null is on the boundary of the parameter space ω 2 0 , so that their estimators cannot have an asymptotically normal distribution. The LM test is free from the boundary value problem, since it forgoes estimation of the constrained parameter ω 2 ;
We now restate the null hypothesis (10), by the parameter transformation ρ 1 = ϕ 1 , ρ 2 = ϕ 2 - ϕ 1 , to
H 0 : ρ 2 = 0 , λ = 1 , ω 2 = 0 ;
the original parameter vector ϑ = ( δ 1 , δ 2 , η , γ , ϕ 1 , ϕ 2 , ω 1 , ω 2 , λ ) is transformed into the new parameter vector θ = ( δ 1 , δ 2 , η , γ , ρ 1 , ρ 2 , ω 1 , ω 2 , λ ) . The vectors of the restricted and unrestricted parameters are denoted by θ 1 = ( ρ 2 , ω 2 , λ ) , θ 2 = ( δ 1 , δ 2 , η , γ , ρ 1 , ω 1 ) respectively, and the log likelihood function for θ is denoted by L ( θ ) log f ( y ˜ T ) , using the likelihood function defined in Equation (8). Then, the LM test statistic for Equation (11) is defined as
L M = s 1 ( θ ^ ) I O P 11 ( θ ^ ) s 1 ( θ ^ )
where the derivative of L ( θ ) with respect to θ 1 is denoted by
s 1 ( θ ^ ) = L ( θ ^ ) θ 1
evaluated at θ ^ , namely the maximum likelihood estimator of θ under the null hypothesis, and I O P 11 ( θ ^ ) is the upper-left three-by-three submatrix of the inverse of the Fisher information matrix evaluated at θ ^ . Then, the LM test statistic follows asymptotically χ 2 distribution with three degrees of freedom, which is the number of constraints.
The Fisher information can be calculated by
I O P ( θ ) = E 0 ( t = 1 T t t )
where
t = log f ( y ˜ t ) θ - log f ( y ˜ t - 1 ) θ
as shown by Hamilton [15] from the identity
t = log f ( y ˜ t | y ˜ t - 1 ) θ .
In practice, I O P 11 ( θ ) is estimated by t = 1 T t t where the unknown parameters θ are estimated by the maximum likelihood method under the null hypothesis.
Hence, we have only to derive the score functions, or the derivatives of log of f ( y ˜ T ) in Equation (8) with respect to θ, to calculate the LM test statistics using Equations (13) and (14). Their derivation is straightforward, by means of integration by parts, but very lengthy so that we provide the detailed derivation in Appendix A.
Denoting E 0 as the expectation under the null hypothesis, their score functions evaluated under the null hypothesis are as follows:
L ( θ ^ ) ρ 1 = t = 2 T E 0 ( h t - 1 u 1 , t - 1 | y ˜ T ) ω ^ 1 , L ( θ ^ ) ρ 2 = t = 2 T j = t T ϕ ^ 2 j - t E 0 ( h t - 1 ξ 2 j | y ˜ T ) - γ ^ E 0 ( h t - 1 ξ 1 j | y ˜ T ) η ^ ( 1 - γ ^ 2 ) , L ( θ ^ ) λ = t = 2 T j = t T ϕ ^ 2 j - t E 0 ( u 1 , t - 1 ξ 2 j | y ˜ T ) - γ ^ E 0 ( u 1 , t - 1 ξ 1 j | y ˜ T ) η ^ ( 1 - γ ^ 2 ) , L ( θ ^ ) ω ^ 2 = t = 1 T k = 1 t ϕ ^ 2 2 ( k - 1 ) E 0 ( ξ 2 t 2 | y ˜ T ) - 2 γ ^ E 0 ( ξ 2 t ξ 1 t | y ˜ T ) + γ ^ 2 E 0 ( ξ 1 t 2 | y ˜ T ) η ^ 2 ( 1 - γ ^ 2 ) 2 - 1 η ^ ( 1 - γ ^ 2 ) + t = 1 T - 1 k = t T - 1 ϕ ^ 2 2 ( k - t ) l = k + 1 T ϕ ^ 2 l - k E 0 ( ξ 2 k ξ 2 l | y ˜ T ) - γ ^ E 0 ( ξ 2 k ξ 1 l | y ˜ T ) - γ ^ E 0 ( ξ 2 l ξ 1 k | y ˜ T ) + γ ^ 2 E 0 ( ξ 1 k ξ 1 l | y ˜ T ) η ^ 2 ( 1 - γ ^ 2 ) 2 ,
L ( θ ^ ) δ 1 = t = 1 T E 0 ( ξ 1 t | y ˜ T ) - γ ^ E 0 ( ξ 2 t | y ˜ T ) η ^ ( 1 - γ ^ 2 ) , L ( θ ^ ) δ 2 = t = 1 T E 0 ( ξ 2 t | y ˜ T ) - γ ^ E 0 ( ξ 1 t | y ˜ T ) η ^ ( 1 - γ ^ 2 ) , L ( θ ^ ) η = t = 1 T 1 2 E 0 ( ξ 2 t 2 | y ˜ ) - 2 γ ^ E 0 ( ξ 1 t ξ 2 t | y ˜ ) + E 0 ( ξ 1 t 2 ) η ^ 2 ( 1 - γ ^ 2 ) - 1 η ^ , L ( θ ^ ) γ = t = 1 T γ ^ 1 - γ ^ 2 - γ ^ E 0 ( ξ 1 t 2 | y ˜ T ) + γ ^ E 0 ( ξ 2 t 2 | y ˜ T ) - ( 1 + γ ^ 2 ) E 0 ( ξ 1 t ξ 2 t | y ˜ T ) η ^ ( 1 - γ ^ 2 ) 2 , L ( θ ^ ) ω ^ 1 = 1 2 t = 2 T E 0 ( u 1 , t - 1 2 | y ˜ T ) ω ^ 1 2 - 1 ω ^ 1 .
Denoting m t = E 0 ( h t | y ˜ T ) , c t s = E 0 ( h t h s | y ˜ T ) , the conditional moments of ξ 1 t , ξ 2 t , u 1 t , h t given y ˜ T under the null hypothesis can be expressed as
E 0 ( h t - 1 u 1 , t - 1 y ˜ T ) = c t - 1 , t - ϕ ^ 1 c t - 1 , t - 1 , E 0 ( u 1 , t - 1 2 y ˜ T ) = c t t - 2 ϕ ^ 1 c t - 1 , t + ϕ ^ 1 2 c t - 1 , t - 1 ,
E 0 ( h t - 1 ξ 1 j y ˜ T ) = m t - 1 ( y 1 j - δ ^ 1 ) - c t - 1 , j , E 0 ( h t - 1 ξ 2 j y ˜ T ) = m t - 1 ( y 2 j - δ ^ 2 ) - c t - 1 , j ,
E 0 ( ξ 1 t y ˜ T ) = y 1 j - δ ^ 1 - m t , E 0 ( ξ 2 t y ˜ T ) = y 2 j - δ ^ 2 - m t ,
E 0 ( ξ 1 k ξ 2 l y ˜ T ) = ( y 1 k - δ ^ 1 ) ( y 2 l - δ ^ 2 ) - m k ( y 2 l - δ ^ 2 ) - ( y 1 k - δ ^ 1 ) m l + c k l
from Equations (2) and (5). See Appendix B for their derivation.

4. Monte Carlo Experiments

To examine the performance of our proposed test, we conduct Monte Carlo experiments using the GAUSS programming language. We use the following data generating process:
  • Generate { e k t 2 } from multivariate log-normal random number generator from GAUSS library (The joint probability density function of e 1 t 2 and e 2 t 2 is expressed as
    f ( e 1 t 2 , e 2 t 2 ) = 1 ( 2 π ) η 2 ( γ 2 - 1 ) e 1 t 2 e 2 t 2 × exp η γ ( log e 1 t 2 - μ ) 2 - 2 η ( log e 1 t 2 - μ ) ( log e 2 t 2 - μ ) + η γ ( log e 2 t 2 - μ ) 2 η 2 ( γ 2 - 1 )
    for 0 < e 1 t 2 , e 2 t 2 < . Note that the library name of multivariate log-normal random number generator in GAUSS is r n d L o g N o r m ( · , · , · , · ) .). In addition, generate { u k t } from Gaussian random number generator.
  • Generate { r k t 2 } , { h k t } from { e k t 2 } , { u k t } . The processes of { r k t 2 } and { h k t } are obtained from the following model;
    r 1 t 2 r 2 t 2 = ψ 1 2 exp ( h 1 t ) 0 0 ψ 2 2 exp ( h 2 t ) e 1 t 2 e 2 t 2 ,
    h 1 t h 2 t = ϕ 1 0 0 ϕ 2 h 1 , t - 1 h 2 , t - 1 + ω 1 1 / 2 0 λ ω 1 1 / 2 ω 2 1 / 2 u 1 , t - 1 u 2 , t - 1
    for t = 1 , , T .
  • Calculate δ k from δ k = μ + log ψ k 2 and generate { y k t } , { ξ k t } from y k t = log r k t 2 , ξ k t = log e k t 2 - μ , respectively.
From the above data generating process, we can define the following linear Gaussian state-space model:
  • Observation equation
    y 1 t y 2 t = δ 1 δ 2 + h 1 t h 2 t + ξ 1 t ξ 2 t ,
  • Transition equation
    h 1 t h 2 t = ϕ 1 0 0 ϕ 2 h 1 , t - 1 h 2 , t - 1 + ω 1 1 / 2 0 λ ω 1 1 / 2 ω 2 1 / 2 u 1 , t - 1 u 2 , t - 1 ,
  • Structure of innovation distribution
    ξ 1 t ξ 2 t u 1 t u 2 t i.i.d. N 0 0 0 0 , η η γ 0 0 η γ η 0 0 0 0 1 0 0 0 0 1 , for t = 1 , , T .
We here report the simulation at μ = - 1.27 , η = π 2 / 2 , γ = 0.1 , ψ 1 = 1 , ψ 2 = 1 , ϕ 1 = 0.7 , ϕ 2 = 0.7 , ω 1 = 0.1 , ω 2 = 0 , λ = 1 and its neighborhood (Note that the parameters ( μ , η , ψ 1 , ψ 2 ) are fixed at ( - 1.27 , π 2 / 2 , 1 , 1 ) throughout this experiment. Then, δ k is expressed simply as δ k = μ + log 1 2 = - 1.27 .). In this experiments, the number of replications is 1000, the sample sizes are T = 500 , 750 , and the initial value of ( h 1 t , h 2 t , u 1 t , u 2 t ) are fixed at ( 1 , 1 , 0 , 0 ) respectively. The initial value of parameters for ML estimation are δ 1 ( 0 ) = 0 , δ 2 ( 0 ) = 0 , η ( 0 ) = π 2 / 2 , γ ( 0 ) = 0 , ϕ 1 ( 0 ) = 0.9 , ϕ 2 ( 0 ) = 0.9 , ω 1 ( 0 ) = 0.1 , ω 2 ( 0 ) = 0 , λ ( 0 ) = 1 . Finally, we generate T + 100 periods of data and use the last T periods by discarding the first 100 periods to reduce the effect of initial values. The procedure (See Durbin and Koopman [16] for a detailed discussion and review of the algorithm of the Kalman filter and smoother.) works as follows.
  • Apply the Kalman filtering to { y t } under the null hypothesis (10) and estimate the parameters by the maximum likelihood (ML) method.
  • Apply the smoothing algorithm to { y t } and obtain the conditional moments of { h t } .
  • Calculate the LM test statistic (12).
Table 1 gives the empirical size of the test, namely the rejection rate of the null hypothesis, under the null for some values of the correlation coefficient γ in Equation (4), the autocorrelation coefficient and variance of the transition equation, namely ϕ 1 and ω 1 in Equation (7). Table 2 reports the empirical power of the test, when the two series have different log volatility series, with different autocorrelation coefficients ϕ 1 and ϕ 2 , or nondegenerate innovations with nonzero value of ω 2 in Equation (7). Each table reports the percentage rate of rejection of the relevant null hypothesis at the asymptotic 5% and 1% significance level. Note that the rejection region of the tests are in the upper tail areas. Our main findings are summarized as follows.
  • Size of the test: For T = 500 , the actual size of the test deviates from the nominal sizes of 5 and 1 percent at most by 3.1 and 1.9 percent, respectively. In particular, when the estimated volatility series are strongly correlated and can be characterized as near unit root processes, the size distortion are increasing. As the sample size increases, namely when T = 750 , the actual rejection rate is closer to the nominal level, as is expected.
  • Power under the alternative: For T = 500 , the rejects rate of the null ranges from 12.9 percent to 58.9 percent for the nominal 5 percent significance level. As the sample size increases, the actual rejection rate is increasing. This results indicate that the absolute value of the difference between the null and the alternative is 0.2 or more, the test statistics can discriminate the hypothesis.
Table 1. Empirical size of the LM test by Monte Carlo experiments.
Table 1. Empirical size of the LM test by Monte Carlo experiments.
γ ϕ 1 ω 1 T = 500 T = 750
5%1%5%1%
0.10.70.16.11.25.50.9
00.70.15.81.95.31.2
−0.10.70.14.60.85.11.1
0.10.90.17.21.85.71.4
0.10.950.18.12.96.51.8
0.10.70.27.11.65.41.1
0.10.70.35.10.74.91.2
Note: The number of replications is 1000. We set ϕ 1 = ϕ 2 , λ = 1 , ω 2 = 0 throughout this experiment.
Table 2. Empirical power of the LM test by Monte Carlo experiments.
Table 2. Empirical power of the LM test by Monte Carlo experiments.
ϕ 2 λ ω 2 T = 500 T = 750
5%1%5%1%
0.51025.510.736.917.7
0.91039.018.452.324.6
0.70.8012.91.315.15.1
0.70.6031.514.144.420.9
0.710.217.95.524.99.0
0.710.458.932.379.152.3
The number of replications is 1000. We set γ = 0.1 , ϕ 1 = 0.7 , ω 1 = 0.1 throughout this experiment.
An important practical conclusion of our simulation is that a rather large sample size, such as more than T = 500 , is necessary to distinguish the null and the alternative. However, datasets with a sample size of 500 or more are easily available in financial analyses.
Summarizing, Monte Carlo results showed that the tests are reliable in terms of both size and power performance.

5. Empirical Analysis

5.1. Univariate Analysis

5.1.1. Data and Descriptive Statistics

We apply the proposed test to the Asian stock market indexes to study volatility spillover in this area. We use daily data (The data will be available from the first author upon request.) of the eleven stock indexes for an about 5-year period from February 23, 2007 to December 8, 2011, excluding the periods when the markets are closed, downloaded from Datastream. For each index, we have a series consisting of 1,250 days. Table 3 gives a complete list of their ticker symbols, index names, and the corresponding markets. These stock market indexes are transformed into daily rates of return. The daily return of index i at time t is calculated as
r i t = log p i t - log p i , t - 1
where p i t denotes the price of index i at time t. In Table 4, the univariate statistics for the returns in percentage terms are presented. The mean for each return is on average around 0.006 percent, while the standard deviation is around 1.704 percent. Mean return of Indonesia is higher than that of the other ten indexes. The return of Hong Kong has the highest volatility among the eleven, whereas the return of Malaysia has the lowest.
Table 3. Ticker symbols, Index names, and the corresponding markets.
Table 3. Ticker symbols, Index names, and the corresponding markets.
Index NumberTickerIndex nameMarket
1HKSPLCIS&P/HKEx Large Cap IndexHong Kong (HKG)
2IBOMBSEBSE 100 IndexIndia (IND)
3CHSASHRShanghai Stock Exchange A Share IndexShanghai (SHA)
4JAKCOMPJakarta Stock Exchange Composite IndexIndonesia (IDN)
5JAPDOWANikkei 225 Stock Average IndexJapan (JPN)
6KOR200IKorea Stock Exchange KOSPI 200 IndexKorea (KOR)
7FBMKLCIFTSE Bursa Malaysia KLCI IndexMalaysia (MYS)
8PSECOMPPhilippines Stock Exchange PSEi IndexPhilippines (PHL)
9TAIWGHTFTSE TWSE Taiwan 50 IndexTaiwan (TWN)
10BNGKS50Stock Exchange of Thailand SET 50 IndexThailand (THA)
11SNGPORIFTSE Straits Times IndexSingapore (SGP)
Table 5 also shows the Ljung-Box statistics to test the null of hypothesis absence of serial correlation in both the levels ( L B ) and in the squares ( L B 2 ) until the 15 lag. Note that the squared returns are used as an approximation to each index’s volatility. The returns in levels show a certain degree of serial correlation, since for eight out of eleven cases, we reject the null hypothesis. Furthermore, all the L B 2 test on the squared returns indicate the precense of serial correlation at any significance level and, therefore the existence of volatility clustering. In this case, the theoretical distribution of the Ljung-Box test is not correct and there is a tendency to over-reject the null.
Table 4. Univariate statistics for the daily index returns.
Table 4. Univariate statistics for the daily index returns.
MarketMaxMinMeanMedianStandard DeviationSkewnessKurtosis
Hong Kong13.353−13.459−0.0040.0572.0700.0888.696
India15.490−11.6890.0190.0961.9120.0229.343
Shanghai9.033−9.261−0.0220.1092.020−0.3665.448
Indonesia7.623−10.9540.0640.1811.740−0.6358.652
Japan13.235−12.111−0.0630.0151.876−0.49510.559
Korea11.540−10.9030.0230.1211.753−0.4228.303
Malaysia4.259−9.9790.0120.0540.958−1.22314.945
Philippines9.365−13.0890.0210.0561.550−0.77310.882
Taiwan6.525−6.735−0.0090.1101.549−0.3425.054
Thailand8.916−12.5630.0360.0531.782−0.4778.384
Singapore7.531−8.696−0.0140.0121.533−0.1186.576
Table 5. Specification tests for the daily returns.
Table 5. Specification tests for the daily returns.
Market LB ( 15 ) LB 2 ( 15 ) Market LB ( 15 ) LB 2 ( 15 )
Hong Kong 23.359 ( 0.077 ) 1101.201 ( 0.000 ) Malaysia 31.058 ( 0.009 ) 93.300 ( 0.000 )
India 30.264 ( 0.011 ) 310.849 ( 0.000 ) Philippines 40.040 ( 0.000 ) 302.482 ( 0.000 )
Shanghai 25.257 ( 0.047 ) 164.299 ( 0.000 ) Taiwan 36.679 ( 0.001 ) 378.213 ( 0.000 )
Indonesia 40.282 ( 0.000 ) 476.883 ( 0.000 ) Thailand 62.000 ( 0.000 ) 795.580 ( 0.000 )
Japan 14.605 ( 0.480 ) 1410.031 ( 0.000 ) Singapore 28.479 ( 0.019 ) 1180.016 ( 0.000 )
Korea 7.941 ( 0.926 ) 816.021 ( 0.000 )
Note: L B ( 15 ) and L B 2 ( 15 ) are the Ljung-Box statistics to test the null of absence of serial correlation in levels and in their squares, respectively, up to the 15th lag. The corresponding p-values are in parentheses. Values in italic represents significance at the 5% level.

5.1.2. Univariate SV Model

From the above discussion, we consider the univariate SV model. First, we transform the return r k t to obtain
y i t = log r i t - T - 1 t = 1 T r i t 2
for i = 1 , , 11 and t = 1 , , 1250 . Second, we consider the following linearlized univariate SV model in state-space form, namely
  • Observation equation
    y i t = δ + h i t + ξ i t ,
  • Transition equation
    h i t = ϕ h i , t - 1 + u i , t - 1 , ϕ < 1
  • Structure of innovation distribution
    ξ i t u i t i.i.d. N 0 0 , η 0 0 ω .
From the above specifications, we evaluate the likelihood and volatilities by means of the Kalman filter and associating smoothing method. Note that the parameter vector is ϑ = ( δ , η , ϕ , ω ) .
Table 6 reports estimation results (The initial value of parameters for ML estimation are δ ( 0 ) = 1 , η ( 0 ) = π 2 / 2 , ϕ ( 0 ) = 0.9 , ω ( 0 ) = 0.1 for all index i.). The estimates of η are larger than π 2 / 2 in almost all index; the only exception is Philippines. The estimates of ϕ indicate highly persistent own dynamics of h i t , with estimated own-lag coefficients; the mean of ϕ ^ is approximately 0.976. The estimates of ω appear relatively small; the mean of ω ^ is approximately 0.037. The estimates of ω in Shanghai, Malaysia, Taiwan, and Singapore are insignificantly different from 0.
Table 6. Estimation results for the univariate SV model.
Table 6. Estimation results for the univariate SV model.
Market δ ^ η ^ ϕ ^ ω ^ Market δ ^ η ^ ϕ ^ ω ^
Hong Kong 9.572 ( 0.321 ) 5.183 ( 0.211 ) 0.988 ( 0.006 ) 0.023 ( 0.009 ) Malaysia 10.095 ( 0.517 ) 5.481 ( 0.233 ) 0.998 ( 0.002 ) 0.008 ( 0.004 )
India 9.635 ( 0.275 ) 5.422 ( 0.232 ) 0.982 ( 0.010 ) 0.031 ( 0.017 ) Philippines 10.012 ( 0.138 ) 4.657 ( 0.225 ) 0.918 ( 0.031 ) 0.120 ( 0.053 )
Shanghai 9.034 ( 0.270 ) 5.354 ( 0.218 ) 0.999 ( 0.001 ) 0.002 ( 0.001 ) Taiwan 10.204 ( 0.230 ) 5.871 ( 0.240 ) 0.982 ( 0.011 ) 0.021 ( 0.013 )
Indonesia 9.943 ( 0.158 ) 5.211 ( 0.230 ) 0.950 ( 0.023 ) 0.062 ( 0.034 ) Thailand 9.842 ( 0.196 ) 4.993 ( 0.213 ) 0.969 ( 0.012 ) 0.041 ( 0.016 )
Japan 9.796 ( 0.238 ) 5.593 ( 0.243 ) 0.975 ( 0.011 ) 0.045 ( 0.019 ) Singapore 10.003 ( 0.265 ) 5.655 ( 0.238 ) 0.983 ( 0.011 ) 0.028 ( 0.016 )
Korea 10.059 ( 0.317 ) 5.205 ( 0.221 ) 0.986 ( 0.009 ) 0.022 ( 0.012 )
Note: The numbers in the parentheses are standard errors for the estimates. Values in italic represents statistical significance at the 5 percent level.
In Figure 1, we plot estimated volatilities of the eleven stock indexes for comparative assessment. Although the volatility series vary heavily over time for each markets, a similar behavior in the eleven indexes is apparent: the volatility processes display very persistence and are very high in the period of the global financial crisis of September 2008. In particular, the volatility processes in Indonesia, Philippines, Taiwan, Thailand, and Singapore display very similar patterns. However, we find that the volatility behavior in Shanghai is a quite distinct: the volatility of Shanghai shows a relatively smooth time process and maintains high value over time.
Figure 1. Stochastic volatilities for the eleven stock indexes. Notes: The figure depicts the stochastic volatilities for the log of squared returns: (from top left to the right and down) Hong Kong, India, Shanghai, Indonesia, Japan, Korea, Malaysia, Philippines, Taiwan, Thailand, and Singapore.
Figure 1. Stochastic volatilities for the eleven stock indexes. Notes: The figure depicts the stochastic volatilities for the log of squared returns: (from top left to the right and down) Hong Kong, India, Shanghai, Indonesia, Japan, Korea, Malaysia, Philippines, Taiwan, Thailand, and Singapore.
Jrfm 06 00031 g001

5.2. Bivariate Analysis

5.2.1. Correlation Analysis

We also examined the correlation matrix of the eleven returns, in levels and in squares, to find possible linkage among different indexes and their volatilities. Table 7 show the correlation matrix for the returns in levels in lower left triangle and for the squares in the upper right triangle.
The returns do exhibit positive and substantial correlations with each other not only in the levels, but also in the squares. Out of the fifty-five correlation coefficients in the levels, forty-five are higher than 0.4, twenty-one are bigger than 0.6, and, among them, one is bigger than 0.8. As might be expected, since their economies are more related to each other, there is a pattern of higher correlation between the neighboring markets; namely, Hong Kong and Korea, Hong Kong and Taiwan, Indonesia and Singapore, Japan and Korea, Korea and Taiwan. Though geographically distant, Hong Kong and Singapore have positive correlation, presumably because they are both important financial markets. Hence, this may not be surprising.
Table 7. Correlation Matrix for returns and squared returns.
Table 7. Correlation Matrix for returns and squared returns.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
Hong Kong10.380.280.550.590.530.270.510.480.610.76
India0.6510.120.340.400.330.190.180.220.340.47
Shanghai0.510.3210.170.140.150.150.220.250.200.20
Indonesia0.680.550.3310.530.450.330.380.540.560.60
Japan0.700.480.340.5610.690.240.420.430.560.50
Korea0.700.510.380.610.7310.260.320.540.360.62
Malaysia0.610.480.360.620.550.5510.310.270.290.31
Philippines0.510.350.270.530.540.480.5610.360.610.32
Taiwan0.650.450.360.580.640.740.560.5110.390.52
Thailand0.640.540.290.610.500.530.520.450.5010.63
Singapore0.840.660.380.690.670.710.630.440.640.641
Note: The table shows the correlation matrix for the daily stock index retuns in levels (lower left triangle) and for the squared returns (upper right triangle). Values in italic face represent correlation coefficients greater than or equal to 0.4.
The correlations between squared returns are naturally related to the correlations between the levels of returns, but can be helpful to discover possible co-movements in their volatilities. From the upper right triangle in Table 7, we can see that there are twenty-three correlation coefficients above 0.4 and almost all of them are between the neighboring markets; Indonesia and Singapore, Japan and Korea, Philippines and Thailand, and Thailand and Singapore. The other strong correlations outside the neighboring markets are those between Hong Kong and Singapore, Japan and Korea, and Korea and Singapore. All these results indicate that there are strong co-movements in the volatilities of the eleven Asian stock market indexes.

5.2.2. Bivariate SV Model

From the above discussion, we apply our test for the presence of the common volatility, using pairwise comparison similar to the methodology used by Engle and Susmel [9]’s; we test for a single-factor common SV, by checking the innovations are perfectly correlated.
First, we consider the following linearlized bivariate SV model in state-space form:
  • Observation equation
    y i t y j t = δ 1 δ 2 + h i t h j t + ξ i t ξ j t
  • Transition equation
    h i t h j t = ϕ 1 0 0 ϕ 2 h i , t - 1 h j , t - 1 + ω 1 1 / 2 0 λ ω 1 1 / 2 ω 2 1 / 2 u i , t - 1 u j , t - 1
  • Structure of innovation distribution
    ξ i t ξ j t u i t u j t i.i.d. N 0 0 0 0 , η η γ 0 0 η γ η 0 0 0 0 1 0 0 0 0 1
for i , j = 1 , , 11 , i j , and t = 1 , , 1250 . Note that the parameter vector is ϑ = ( δ 1 , δ 2 , η , γ , ϕ 1 , ϕ 2 , ω 1 , ω 2 , λ ) . Second, we apply the Kalman filtering to the bivariate SV model under the null hypothesis (10). Finally, we calculate the LM test statistics (12) from the smoothing recursions.
Note that the ordering of index i and j may affect the estimation results and the LM test statistics as per the set up in Equations (15)–(17). Hence, we conduct estimation and run the test for the index pair of i < j , as well as the pair of i > j . Further, we also conduct the estimation under the alternative hypothesis. However, due to space limitations, we do not mention the results in this paper (The results will be available from the first author upon request.).
Table 8 through 13 report estimation results under the null hypothesis (The initial value of parameters for ML estimation are δ 1 ( 0 ) = 0 , δ 2 ( 0 ) = 0 , η ( 0 ) = π 2 / 2 , γ ( 0 ) = 0 , ϕ 1 ( 0 ) = 0 . 9 , ϕ 1 ( 0 ) = 0 . 9 , ω 1 ( 0 ) = 0 . 1 , ω 2 ( 0 ) = 0 , λ ( 0 ) = 1 for all index pair i and j.). Table 11 shows the estimates of γ indicate the positive correlation between y i t and y j t ; the mean of γ ^ is approximately 0.11. Table 12 and Table 13 also show the estimates of ϕ 1 and ω 1 are close to that of ϕ and ω in the univariate SV models; the mean of ϕ ^ 1 is approximately 0.97, the mean of ω 1 is approximately 0.04.
Table 14 presents the result of the LM test proposed in this paper. We have some evidence of the single common volatility for the pairs under investigation; Out of the fifty-five pairs in the upper right triangle in Table 14, we can see that seventeen pairs are not rejected the null hypothesis. In addition, from the lower left triangle in Table 14, twenty-one pairs are not rejected the null.
Table 15 shows the list which cannot reject the null at the 5% significance level for the pair of i < j , as well as the pair of i > j in Table 14. These pairs have a possibility to have a single common volatility. Table 15 suggests that the markets in Southeast Asia are frequently overridden by a single intermarket volatility; the null hypothesis is not rejected for the pair of Indonesia and Philippines, Indonesia and Thailand, Indonesia and Singapore, Malaysia and Singapore, Philippines and Thailand, and Philippines and Singapore. One possible explanation for the results is that these pairs share the same information, since these markets lie within the same time zone. Moreover, we cannot reject the null for the pair of Hong Kong and India, Hong Kong and Singapore, India and Singapore, Korea and Singapore. This suggests that these important financial markets are also overridden by a single intermarket volatility factor in the period of the financial crisis.
Table 8. Estimation results of δ 1 in the bivariate SV model.
Table 8. Estimation results of δ 1 in the bivariate SV model.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
HKG 9.43 ( 0.30 ) 9.48 ( 0.36 ) 9.70 ( 0.48 ) 9.52 ( 0.26 ) 9.66 ( 0.31 ) 9.13 ( 0.42 ) 9.43 ( 0.26 ) 9.64 ( 0.27 ) 9.60 ( 0.25 ) 9.46 ( 0.28 )
IND 9.53 ( 0.30 ) 9.54 ( 0.30 ) 9.75 ( 0.18 ) 9.70 ( 0.27 ) 9.75 ( 0.22 ) 9.00 ( 0.65 ) 9.63 ( 0.23 ) 9.73 ( 0.21 ) 9.76 ( 0.22 ) 9.68 ( 0.25 )
SHA 9.40 ( 0.30 ) 9.31 ( 0.29 ) 9.51 ( 0.29 ) 9.39 ( 0.25 ) 9.68 ( 0.27 ) 8.66 ( 0.35 ) 8.84 ( 0.39 ) 9.65 ( 0.25 ) 9.63 ( 0.25 ) 9.33 ( 0.30 )
IDN 9.89 ( 0.23 ) 9.91 ( 0.17 ) 9.87 ( 0.22 ) 9.93 ( 0.20 ) 9.94 ( 0.19 ) 9.79 ( 0.25 ) 9.92 ( 0.14 ) 9.94 ( 0.17 ) 9.97 ( 0.18 ) 9.92 ( 0.19 )
JPN 9.72 ( 0.29 ) 9.76 ( 0.25 ) 9.72 ( 0.25 ) 9.82 ( 0.20 ) 9.91 ( 0.22 ) 9.56 ( 0.31 ) 9.77 ( 0.18 ) 9.86 ( 0.21 ) 9.87 ( 0.20 ) 9.75 ( 0.24 )
KOR 9.93 ( 0.29 ) 9.91 ( 0.23 ) 9.96 ( 0.25 ) 9.88 ( 0.21 ) 9.90 ( 0.23 ) 9.71 ( 0.27 ) 9.83 ( 0.18 ) 10.47 ( 0.40 ) 10.01 ( 0.23 ) 9.87 ( 0.24 )
MYS 10.60 ( 0.39 ) 10.33 ( 0.50 ) 10.32 ( 0.37 ) 10.96 ( 0.22 ) 10.78 ( 0.29 ) 10.99 ( 0.28 ) 9.68 ( 0.79 ) 10.89 ( 0.27 ) 10.98 ( 0.20 ) 10.43 ( 0.55 )
PHL 9.88 ( 0.23 ) 9.91 ( 0.22 ) 9.36 ( 0.33 ) 10.02 ( 0.15 ) 9.98 ( 0.16 ) 10.02 ( 0.18 ) 8.57 ( 0.63 ) 10.03 ( 0.18 ) 9.97 ( 0.16 ) 9.93 ( 0.22 )
TWN 10.20 ( 0.30 ) 10.09 ( 0.25 ) 10.23 ( 0.25 ) 10.13 ( 0.20 ) 10.15 ( 0.21 ) 10.81 ( 0.34 ) 9.90 ( 0.27 ) 10.09 ( 0.17 ) 10.24 ( 0.22 ) 10.21 ( 0.26 )
THA 9.84 ( 0.26 ) 9.78 ( 0.24 ) 9.91 ( 0.24 ) 9.85 ( 0.20 ) 9.83 ( 0.21 ) 10.80 ( 0.44 ) 9.66 ( 0.21 ) 9.74 ( 0.15 ) 10.54 ( 0.38 ) 9.83 ( 0.22 )
SGP 9.93 ( 0.30 ) 10.02 ( 0.25 ) 9.92 ( 0.28 ) 10.08 ( 0.19 ) 10.05 ( 0.24 ) 10.17 ( 0.24 ) 9.44 ( 0.65 ) 9.98 ( 0.22 ) 10.09 ( 0.22 ) 10.13 ( 0.21 )
Table 9. Estimation results of δ 2 in the bivariate SV model.
Table 9. Estimation results of δ 2 in the bivariate SV model.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
HKG 9.43 ( 0.30 ) 9.48 ( 0.36 ) 9.70 ( 0.48 ) 9.52 ( 0.26 ) 9.66 ( 0.31 ) 9.13 ( 0.42 ) 9.43 ( 0.26 ) 9.64 ( 0.27 ) 9.60 ( 0.25 ) 9.46 ( 0.28 )
IND 9.41 ( 0.30 ) 9.22 ( 0.30 ) 9.89 ( 0.18 ) 9.69 ( 0.27 ) 9.85 ( 0.22 ) 10.31 ( 0.65 ) 9.87 ( 0.23 ) 10.06 ( 0.22 ) 9.72 ( 0.22 ) 9.96 ( 0.25 )
SHA 9.46 ( 0.30 ) 9.54 ( 0.29 ) 9.88 ( 0.29 ) 9.68 ( 0.25 ) 9.96 ( 0.27 ) 10.22 ( 0.35 ) 9.27 ( 0.38 ) 10.23 ( 0.25 ) 9.86 ( 0.25 ) 9.87 ( 0.30 )
IDN 9.55 ( 0.23 ) 9.71 ( 0.17 ) 9.43 ( 0.22 ) 9.78 ( 0.19 ) 9.84 ( 0.19 ) 10.93 ( 0.25 ) 10.01 ( 0.15 ) 10.10 ( 0.17 ) 9.80 ( 0.18 ) 10.04 ( 0.19 )
JPN 9.50 ( 0.29 ) 9.63 ( 0.25 ) 9.35 ( 0.25 ) 9.89 ( 0.20 ) 9.90 ( 0.22 ) 10.82 ( 0.31 ) 9.92 ( 0.18 ) 10.11 ( 0.21 ) 9.79 ( 0.20 ) 10.02 ( 0.24 )
KOR 9.62 ( 0.29 ) 9.69 ( 0.23 ) 9.59 ( 0.25 ) 9.92 ( 0.21 ) 9.78 ( 0.23 ) 10.95 ( 0.27 ) 9.96 ( 0.18 ) 10.70 ( 0.40 ) 9.87 ( 0.23 ) 10.05 ( 0.24 )
MYS 9.15 ( 0.39 ) 8.97 ( 0.50 ) 8.71 ( 0.37 ) 9.81 ( 0.22 ) 9.52 ( 0.29 ) 9.73 ( 0.27 ) 8.69 ( 0.79 ) 9.92 ( 0.27 ) 9.65 ( 0.20 ) 9.41 ( 0.55 )
PHL 9.41 ( 0.23 ) 9.59 ( 0.22 ) 8.81 ( 0.33 ) 9.88 ( 0.16 ) 9.73 ( 0.16 ) 9.78 ( 0.18 ) 9.63 ( 0.63 ) 10.06 ( 0.18 ) 9.70 ( 0.16 ) 9.97 ( 0.22 )
TWN 9.72 ( 0.30 ) 9.72 ( 0.25 ) 9.63 ( 0.25 ) 9.97 ( 0.20 ) 9.83 ( 0.21 ) 10.57 ( 0.34 ) 10.93 ( 0.27 ) 10.00 ( 0.17 ) 9.86 ( 0.22 ) 10.16 ( 0.25 )
THA 9.60 ( 0.26 ) 9.73 ( 0.24 ) 9.55 ( 0.24 ) 9.92 ( 0.20 ) 9.80 ( 0.21 ) 10.86 ( 0.44 ) 10.99 ( 0.21 ) 9.94 ( 0.15 ) 10.87 ( 0.38 ) 10.07 ( 0.22 )
SGP 9.43 ( 0.30 ) 9.62 ( 0.25 ) 9.26 ( 0.28 ) 9.92 ( 0.19 ) 9.71 ( 0.23 ) 9.86 ( 0.24 ) 10.44 ( 0.65 ) 9.90 ( 0.22 ) 10.09 ( 0.23 ) 9.77 ( 0.21 )
Note: Paraemter estimates are calculated under the null hypohtesis H 0 : ϕ 1 = ϕ 2 , λ = 1 , ω 2 = 0 . The numbers in the parentheses are standard errors for the estimates. Values in italic represents statistical significance at the 5 percent level.
Table 10. Estimation results of η in the bivariate SV model.
Table 10. Estimation results of η in the bivariate SV model.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
HKG 5.38 ( 0.17 ) 5.45 ( 0.16 ) 5.08 ( 0.16 ) 5.49 ( 0.16 ) 5.26 ( 0.15 ) 5.46 ( 0.16 ) 5.07 ( 0.15 ) 5.79 ( 0.17 ) 5.21 ( 0.16 ) 5.08 ( 0.15 )
IND 5.26 ( 0.16 ) 5.86 ( 0.16 ) 5.67 ( 0.17 ) 6.10 ( 0.19 ) 5.57 ( 0.17 ) 5.84 ( 0.18 ) 5.36 ( 0.16 ) 5.74 ( 0.17 ) 5.42 ( 0.17 ) 5.53 ( 0.17 )
SHA 5.46 ( 0.17 ) 5.50 ( 0.16 ) 5.43 ( 0.16 ) 5.70 ( 0.17 ) 5.74 ( 0.17 ) 5.60 ( 0.16 ) 5.29 ( 0.17 ) 5.73 ( 0.16 ) 5.60 ( 0.16 ) 5.48 ( 0.16 )
IDN 5.50 ( 0.16 ) 5.50 ( 0.17 ) 5.40 ( 0.16 ) 5.60 ( 0.18 ) 5.39 ( 0.15 ) 5.40 ( 0.16 ) 5.33 ( 0.16 ) 5.59 ( 0.17 ) 5.24 ( 0.16 ) 5.23 ( 0.15 )
JPN 5.40 ( 0.16 ) 5.75 ( 0.17 ) 5.56 ( 0.16 ) 5.38 ( 0.16 ) 5.56 ( 0.16 ) 5.62 ( 0.17 ) 5.28 ( 0.16 ) 6.01 ( 0.18 ) 5.55 ( 0.17 ) 5.54 ( 0.16 )
KOR 5.17 ( 0.15 ) 5.63 ( 0.17 ) 5.59 ( 0.16 ) 5.33 ( 0.16 ) 5.47 ( 0.16 ) 5.56 ( 0.17 ) 5.46 ( 0.17 ) 5.36 ( 0.15 ) 5.26 ( 0.16 ) 5.38 ( 0.15 )
MYS 5.54 ( 0.16 ) 5.66 ( 0.17 ) 5.47 ( 0.17 ) 5.54 ( 0.17 ) 5.56 ( 0.17 ) 5.54 ( 0.16 ) 5.15 ( 0.15 ) 5.71 ( 0.16 ) 5.50 ( 0.17 ) 5.37 ( 0.16 )
PHL 5.22 ( 0.16 ) 5.29 ( 0.16 ) 5.34 ( 0.15 ) 5.08 ( 0.16 ) 5.39 ( 0.17 ) 5.29 ( 0.16 ) 5.28 ( 0.16 ) 5.67 ( 0.18 ) 4.94 ( 0.16 ) 5.24 ( 0.16 )
TWN 5.41 ( 0.16 ) 5.48 ( 0.16 ) 5.68 ( 0.16 ) 5.49 ( 0.16 ) 6.07 ( 0.17 ) 5.38 ( 0.16 ) 5.58 ( 0.16 ) 5.45 ( 0.17 ) 5.75 ( 0.17 ) 5.65 ( 0.16 )
THA 5.17 ( 0.15 ) 5.37 ( 0.16 ) 5.66 ( 0.17 ) 5.26 ( 0.16 ) 5.51 ( 0.16 ) 5.25 ( 0.15 ) 5.48 ( 0.16 ) 5.03 ( 0.16 ) 5.64 ( 0.17 ) 5.28 ( 0.16 )
SGP 5.08 ( 0.15 ) 5.50 ( 0.17 ) 5.64 ( 0.16 ) 5.32 ( 0.15 ) 5.64 ( 0.18 ) 5.52 ( 0.16 ) 5.52 ( 0.16 ) 5.15 ( 0.15 ) 5.88 ( 0.17 ) 5.50 ( 0.16 )
Table 11. Estimation results of γ in the bivariate SV model.
Table 11. Estimation results of γ in the bivariate SV model.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
HKG 0.15 ( 0.03 ) 0.11 ( 0.03 ) 0.20 ( 0.00 ) 0.20 ( 0.00 ) 0.20 ( 0.00 ) 0.13 ( 0.03 ) 0.09 ( 0.03 ) 0.18 ( 0.03 ) 0.17 ( 0.03 ) 0.20 ( 0.00 )
IND 0.17 ( 0.03 ) 0.01 ( 0.03 ) 0.16 ( 0.03 ) 0.06 ( 0.02 ) 0.06 ( 0.03 ) 0.08 ( 0.02 ) 0.01 ( 0.03 ) 0.07 ( 0.02 ) 0.08 ( 0.03 ) 0.15 ( 0.03 )
SHA 0.10 ( 0.02 ) 0.02 ( 0.03 ) 0.02 ( 0.02 ) 0.05 ( 0.03 ) 0.05 ( 0.02 ) 0.02 ( 0.02 ) 0.02 ( 0.03 ) 0.04 ( 0.03 ) 0.01 ( 0.02 ) 0.05 ( 0.03 )
IDN 0.20 ( 0.00 ) 0.16 ( 0.03 ) 0.02 ( 0.03 ) 0.11 ( 0.03 ) 0.17 ( 0.03 ) 0.16 ( 0.03 ) 0.08 ( 0.03 ) 0.18 ( 0.02 ) 0.10 ( 0.03 ) 0.20 ( 0.00 )
JPN 0.20 ( 0.00 ) 0.07 ( 0.02 ) 0.04 ( 0.03 ) 0.10 ( 0.03 ) 0.20 ( 0.00 ) 0.11 ( 0.03 ) 0.11 ( 0.03 ) 0.16 ( 0.02 ) 0.06 ( 0.03 ) 0.16 ( 0.03 )
KOR 0.20 ( 0.00 ) 0.07 ( 0.03 ) 0.07 ( 0.03 ) 0.18 ( 0.02 ) 0.20 ( 0.00 ) 0.16 ( 0.03 ) 0.07 ( 0.03 ) 0.20 ( 0.00 ) 0.09 ( 0.03 ) 0.20 ( 0.00 )
MYS 0.15 ( 0.03 ) 0.10 ( 0.03 ) 0.06 ( 0.03 ) 0.16 ( 0.03 ) 0.13 ( 0.02 ) 0.17 ( 0.03 ) 0.05 ( 0.03 ) 0.11 ( 0.03 ) 0.08 ( 0.03 ) 0.15 ( 0.02 )
PHL 0.09 ( 0.03 ) 0.04 ( 0.03 ) 0.01 ( 0.03 ) 0.10 ( 0.03 ) 0.10 ( 0.03 ) 0.05 ( 0.03 ) 0.06 ( 0.03 ) 0.11 ( 0.03 ) 0.01 ( 0.03 ) 0.06 ( 0.03 )
TWN 0.20 ( 0.00 ) 0.06 ( 0.02 ) 0.05 ( 0.02 ) 0.17 ( 0.03 ) 0.18 ( 0.03 ) 0.20 ( 0.00 ) 0.10 ( 0.03 ) 0.09 ( 0.03 ) 0.07 ( 0.02 ) 0.17 ( 0.03 )
THA 0.16 ( 0.03 ) 0.11 ( 0.03 ) 0.00 ( 0.03 ) 0.10 ( 0.03 ) 0.07 ( 0.03 ) 0.09 ( 0.03 ) 0.07 ( 0.03 ) 0.00 ( 0.03 ) 0.07 ( 0.03 ) 0.16 ( 0.03 )
SGP 0.20 ( 0.00 ) 0.18 ( 0.03 ) 0.06 ( 0.03 ) 0.20 ( 0.00 ) 0.14 ( 0.03 ) 0.19 ( 0.03 ) 0.13 ( 0.03 ) 0.06 ( 0.03 ) 0.17 ( 0.02 ) 0.13 ( 0.02 )
Note: Paraemter estimates are calculated under the null hypohtesis H 0 : ϕ 1 = ϕ 2 , λ = 1 , ω 2 = 0 . The numbers in the parentheses are standard errors for the estimates. Values in italic represents statistical significance at the 5 percent level.
Table 12. Estimation results of ϕ 1 in the bivariate SV model.
Table 12. Estimation results of ϕ 1 in the bivariate SV model.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
HKG 0.98 ( 0.01 ) 0.99 ( 0.00 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.01 ) 0.99 ( 0.01 ) 0.97 ( 0.01 ) 0.99 ( 0.01 ) 0.98 ( 0.01 ) 0.98 ( 0.01 )
IND 0.99 ( 0.01 ) 0.99 ( 0.01 ) 0.96 ( 0.02 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.01 ) 0.97 ( 0.01 ) 0.98 ( 0.01 ) 0.97 ( 0.01 ) 0.98 ( 0.01 )
SHA 0.99 ( 0.00 ) 1.00 ( 0.00 ) 0.99 ( 0.01 ) 0.99 ( 0.01 ) 0.99 ( 0.00 ) 1.00 ( 0.00 ) 1.00 ( 0.00 ) 1.00 ( 0.00 ) 0.99 ( 0.01 ) 1.00 ( 0.00 )
IDN 0.97 ( 0.01 ) 0.96 ( 0.01 ) 0.99 ( 0.01 ) 0.97 ( 0.01 ) 0.97 ( 0.01 ) 0.98 ( 0.01 ) 0.92 ( 0.03 ) 0.97 ( 0.01 ) 0.96 ( 0.01 ) 0.97 ( 0.01 )
JPN 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.01 ) 0.97 ( 0.01 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.96 ( 0.02 ) 0.98 ( 0.01 ) 0.97 ( 0.01 ) 0.98 ( 0.01 )
KOR 0.99 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.00 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.01 ) 0.96 ( 0.02 ) 0.98 ( 0.01 ) 0.97 ( 0.01 ) 0.98 ( 0.01 )
MYS 0.99 ( 0.01 ) 0.99 ( 0.01 ) 1.00 ( 0.00 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.01 ) 0.96 ( 0.02 ) 0.99 ( 0.01 ) 0.97 ( 0.01 ) 0.99 ( 0.01 )
PHL 0.98 ( 0.01 ) 0.97 ( 0.01 ) 1.00 ( 0.00 ) 0.94 ( 0.02 ) 0.96 ( 0.02 ) 0.96 ( 0.02 ) 0.96 ( 0.02 ) 0.96 ( 0.02 ) 0.95 ( 0.02 ) 0.97 ( 0.01 )
TWN 0.99 ( 0.01 ) 0.98 ( 0.01 ) 1.00 ( 0.00 ) 0.97 ( 0.01 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.01 ) 0.96 ( 0.02 ) 0.97 ( 0.01 ) 0.99 ( 0.00 )
THA 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.01 ) 0.97 ( 0.01 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.95 ( 0.02 ) 0.98 ( 0.01 ) 0.98 ( 0.01 )
SGP 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.00 ) 0.97 ( 0.01 ) 0.98 ( 0.01 ) 0.98 ( 0.01 ) 0.99 ( 0.00 ) 0.97 ( 0.01 ) 0.98 ( 0.01 ) 0.97 ( 0.01 )
Table 13. Estimation results of ω 1 in the bivariate SV model.
Table 13. Estimation results of ω 1 in the bivariate SV model.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
HKG 0.03 ( 0.01 ) 0.01 ( 0.00 ) 0.03 ( 0.01 ) 0.03 ( 0.01 ) 0.02 ( 0.01 ) 0.02 ( 0.01 ) 0.04 ( 0.01 ) 0.02 ( 0.01 ) 0.03 ( 0.01 ) 0.03 ( 0.01 )
IND 0.02 ( 0.01 ) 0.01 ( 0.01 ) 0.06 ( 0.02 ) 0.04 ( 0.01 ) 0.04 ( 0.01 ) 0.02 ( 0.01 ) 0.04 ( 0.01 ) 0.03 ( 0.01 ) 0.03 ( 0.01 ) 0.04 ( 0.02 )
SHA 0.01 ( 0.00 ) 0.00 ( 0.00 ) 0.01 ( 0.00 ) 0.01 ( 0.01 ) 0.01 ( 0.00 ) 0.01 ( 0.00 ) 0.01 ( 0.00 ) 0.00 ( 0.00 ) 0.01 ( 0.00 ) 0.01 ( 0.00 )
IDN 0.04 ( 0.01 ) 0.05 ( 0.02 ) 0.01 ( 0.01 ) 0.04 ( 0.01 ) 0.03 ( 0.01 ) 0.02 ( 0.01 ) 0.11 ( 0.04 ) 0.03 ( 0.01 ) 0.05 ( 0.02 ) 0.04 ( 0.01 )
JPN 0.03 ( 0.01 ) 0.04 ( 0.01 ) 0.01 ( 0.01 ) 0.04 ( 0.01 ) 0.03 ( 0.01 ) 0.02 ( 0.01 ) 0.05 ( 0.02 ) 0.03 ( 0.01 ) 0.04 ( 0.01 ) 0.04 ( 0.01 )
KOR 0.02 ( 0.01 ) 0.04 ( 0.01 ) 0.01 ( 0.00 ) 0.03 ( 0.01 ) 0.04 ( 0.01 ) 0.01 ( 0.01 ) 0.05 ( 0.03 ) 0.02 ( 0.01 ) 0.03 ( 0.01 ) 0.03 ( 0.01 )
MYS 0.02 ( 0.01 ) 0.01 ( 0.01 ) 0.00 ( 0.00 ) 0.02 ( 0.01 ) 0.02 ( 0.01 ) 0.02 ( 0.01 ) 0.06 ( 0.03 ) 0.01 ( 0.01 ) 0.03 ( 0.02 ) 0.01 ( 0.01 )
PHL 0.04 ( 0.01 ) 0.04 ( 0.02 ) 0.00 ( 0.00 ) 0.08 ( 0.03 ) 0.06 ( 0.02 ) 0.05 ( 0.02 ) 0.05 ( 0.02 ) 0.04 ( 0.02 ) 0.07 ( 0.02 ) 0.04 ( 0.02 )
TWN 0.01 ( 0.01 ) 0.02 ( 0.01 ) 0.00 ( 0.00 ) 0.03 ( 0.01 ) 0.03 ( 0.01 ) 0.02 ( 0.01 ) 0.01 ( 0.01 ) 0.04 ( 0.02 ) 0.03 ( 0.01 ) 0.02 ( 0.01 )
THA 0.02 ( 0.01 ) 0.03 ( 0.01 ) 0.01 ( 0.00 ) 0.04 ( 0.01 ) 0.03 ( 0.01 ) 0.02 ( 0.01 ) 0.02 ( 0.01 ) 0.06 ( 0.02 ) 0.02 ( 0.01 ) 0.02 ( 0.01 )
SGP 0.03 ( 0.01 ) 0.04 ( 0.02 ) 0.01 ( 0.01 ) 0.05 ( 0.01 ) 0.04 ( 0.01 ) 0.04 ( 0.02 ) 0.02 ( 0.01 ) 0.04 ( 0.01 ) 0.02 ( 0.01 ) 0.03 ( 0.01 )
Note: Paraemter estimates are calculated under the null hypohtesis H 0 : ϕ 1 = ϕ 2 , λ = 1 , ω 2 = 0 . The numbers in the parentheses are standard errors for the estimates. Values in italic represents statistical significance at the 5 percent level.
Table 14. LM test statistics for the bivariate SV model with a single common volatility.
Table 14. LM test statistics for the bivariate SV model with a single common volatility.
MarketHKGINDSHAIDNJPNKORMYSPHLTWNTHASGP
Hong Kong 5.87 ( 0.12 ) 16.15 ( 0.00 ) 15.85 ( 0.00 ) 7.87 ( 0.05 ) 10.50 ( 0.01 ) 12.25 ( 0.01 ) 17.62 ( 0.00 ) 11.79 ( 0.01 ) 20.50 ( 0.00 ) 2.17 ( 0.54 )
India 4.76 ( 0.19 ) 13.47 ( 0.00 ) 9.12 ( 0.03 ) 12.40 ( 0.01 ) 9.69 ( 0.02 ) 12.17 ( 0.01 ) 10.27 ( 0.02 ) 8.54 ( 0.04 ) 8.14 ( 0.04 ) 3.14 ( 0.37 )
Shanghai 15.71 ( 0.00 ) 11.57 ( 0.01 ) 14.82 ( 0.00 ) 28.21 ( 0.00 ) 25.44 ( 0.00 ) 13.56 ( 0.00 ) 12.06 ( 0.01 ) 11.62 ( 0.01 ) 30.07 ( 0.00 ) 10.22 ( 0.02 )
Indonesia 16.05 ( 0.00 ) 16.45 ( 0.00 ) 19.05 ( 0.00 ) 21.95 ( 0.00 ) 11.65 ( 0.01 ) 11.46 ( 0.01 ) 0.82 ( 0.85 ) 2.99 ( 0.39 ) 1.36 ( 0.71 ) 7.23 ( 0.06 )
Japan 5.25 ( 0.15 ) 11.56 ( 0.01 ) 23.54 ( 0.00 ) 12.09 ( 0.01 ) 9.81 ( 0.02 ) 20.83 ( 0.00 ) 16.69 ( 0.00 ) 1.50 ( 0.68 ) 18.92 ( 0.00 ) 11.72 ( 0.01 )
Korea 14.87 ( 0.00 ) 8.98 ( 0.03 ) 26.10 ( 0.00 ) 11.86 ( 0.01 ) 11.65 ( 0.01 ) 12.38 ( 0.01 ) 8.60 ( 0.04 ) 4.11 ( 0.25 ) 3.73 ( 0.29 ) 6.07 ( 0.11 )
Malaysia 13.21 ( 0.00 ) 7.57 ( 0.06 ) 12.86 ( 0.00 ) 12.20 ( 0.01 ) 20.56 ( 0.00 ) 8.26 ( 0.04 ) 10.59 ( 0.01 ) 4.77 ( 0.19 ) 19.65 ( 0.00 ) 3.09 ( 0.38 )
Philippines 13.26 ( 0.00 ) 8.04 ( 0.05 ) 4.72 ( 0.19 ) 1.42 ( 0.70 ) 21.23 ( 0.00 ) 14.92 ( 0.00 ) 9.26 ( 0.03 ) 9.22 ( 0.03 ) 3.84 ( 0.28 ) 6.35 ( 0.10 )
Taiwan 13.56 ( 0.00 ) 4.97 ( 0.17 ) 8.19 ( 0.04 ) 0.28 ( 0.96 ) 5.75 ( 0.12 ) 1.04 ( 0.79 ) 1.59 ( 0.66 ) 1.54 ( 0.67 ) 5.47 ( 0.14 ) 2.84 ( 0.42 )
Thailand 12.78 ( 0.01 ) 8.48 ( 0.04 ) 30.63 ( 0.00 ) 4.17 ( 0.24 ) 20.87 ( 0.00 ) 5.11 ( 0.16 ) 24.85 ( 0.00 ) 2.83 ( 0.42 ) 7.43 ( 0.06 ) 10.34 ( 0.02 )
Singapore 3.67 ( 0.30 ) 0.79 ( 0.85 ) 14.34 ( 0.00 ) 4.63 ( 0.20 ) 10.31 ( 0.02 ) 1.59 ( 0.66 ) 2.95 ( 0.40 ) 9.66 ( 0.02 ) 4.67 ( 0.20 ) 8.49 ( 0.04 )
Note: The variables in each cell are the LM test statistics under the null hypohtesis H 0 : ϕ 1 = ϕ 2 , λ = 1 , ω 2 = 0 . The corresponding p-values are paretheses. Values in italic represents statistical significance at the 5 percent level.
Table 15. The pairs which cannot the null hypothesis of the single common volatility.
Table 15. The pairs which cannot the null hypothesis of the single common volatility.
No.Pair of marketsNo.Pair of markets
1Hong Kong and India9Korea and Taiwan
2Hong Kong and Singapore10Korea and Thailand
3India and Singapore11Korea and Singapore
4Indonesia and Philippines12Malaysia and Taiwan
5Indonesia and Taiwan13Malaysia and Singapore
6Indonesia and Thailand14Philippines and Thailand
7Indonesia and Singapore15Taiwan and Thailand
8Japan and Taiwan16Taiwan and Singapore
It is worth noting that the pairs between Taiwan and the other six indexes cannot reject the null; namely, Indonesia and Taiwan, Japan and Taiwan, Korea and Taiwan, Malaysia and Taiwan, Taiwan and Thailand, Taiwan and Singapore. This suggests that the idiosyncratic volatility factor of Taiwan are negligible in size in the period of the financial crisis. Then, the pairs between Taiwan and the other six indexes are looked as if they have a single common volatility process in common. On the other hand, we strongly reject the null for all the pairs between Shanghai and the other ten indexes. This confirms that the stock market in Shanghai has a powerful idiosyncratic volatility.
Summarizing, the findings in empirical analysis are as follows; the volatility process displays a high degree of similarities in the Asian stock market indexes. The only exception is Shanghai. Moreover, we cannot the null hypothesis of the single common volatility in the pairs of Southeast Asia. This suggests that the markets in Southeast Asia are frequently overridden by a single intermarket volatility.

6. Conclusions

This paper provided the Lagrange multiplier test statistic for the null hypothesis that the bivariate time series has only a single stochastic volatility factor. The test statistic is constructed using smoothing algorithm for linear state-space models under the assumption that the log of squared measurement error is normally distributed. The finite sample properties of the test has reliable size and power, when the sample size is moderately large. In an empirical analysis we have shown that the test can capture key features of the Asian Stock markets.

Acknowledgements

The authors are grateful to Professor Daisuke Nagakura and seminar participants for helpful comments and discussions at Hokkaido University. Two anonymous referees also made valuable comments that helped the authors to improve the paper. The work was partly supported by Grant-in-Aid for Scientific Research (No. 21530195, 23730219) from Japan Society for the Promotion of Science. Needless to say, the authors are solely responsible for any errors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. R.F. Engle. “Autoregressive conditional heteroskedasticity with estimates of the variance of U.K. inflation.” Econometrica 50 (1982): 987–1008. [Google Scholar] [CrossRef]
  2. H. Chernoff. “On the distribution of the likelihood ratio.” Ann. Stat. 25 (1954): 573–578. [Google Scholar] [CrossRef]
  3. P.A.P. Moran. “Maximum likelihood estimation in non-standard conditions.” Proc. Camb. Philos. Soc. 70 (1971): 441–450. [Google Scholar] [CrossRef]
  4. D. Chant. “On asymptotic tests of composite hypotheses in non-standard conditions.” Biometrica 61 (1974): 291–299. [Google Scholar] [CrossRef]
  5. S.J. Taylor. “Financial Returns Modelled by the Product of Two Stochastic Processes—a Study of Daily Sugar Prices.” In Time Series Analysis: Theory and Practice 1. Amsterdam, The Netherlands: North-Holland, 1982. [Google Scholar]
  6. E. Jacquier, N.G. Polson, and P.E. Rossi. “Bayesian analysis of stochastic volatility models.” J. Bus. Econ. Stat. 12 (1994): 371–389. [Google Scholar]
  7. S. Kim, N. Shephard, and S. Chib. “Stochastic volatility: Likelihood inference and comparison with ARCH models.” Rev. Econ. Stud. 65 (1998): 361–393. [Google Scholar] [CrossRef]
  8. M. Asai, M. McAleer, and J. Yu. “Multivariate stochastic volatileity: A review.” Econ. Rev. 25 (2006): 145–175. [Google Scholar] [CrossRef]
  9. R.F. Engle, and R. Susmel. “Common volatility in international equity markets.” J. Bus. Econ. Stat. 11 (1993): 167–176. [Google Scholar]
  10. F.J. Breidt, and A.L. Carriquiry. “Improved Quasi-Maximum Likelihood Estimation for Stochastic Volatility Models.” In Modelling and Prediction: Honoring Seymour Geisser. Berlin/Heidelberg, Germany: Springer, 1996, pp. 228–247. [Google Scholar]
  11. A. Harvey, E. Ruiz, and N. Shephard. “Multivariate stochastic variance models.” Rev. Econ. Stud. 61 (1994): 247–264. [Google Scholar] [CrossRef]
  12. D.B. Nelson. “The Time Series Behaviour of Stock Market Volatility and Returns.” Ph.D. Dissertation, Massachusetts Institute of Technology, Cambridge, MA, USA, 1988. [Google Scholar]
  13. A. Melino, and S.M. Turnbull. “Pricing foreign currency options with stochastic volatility.” J. Econ. 45 (1990): 239–265. [Google Scholar] [CrossRef]
  14. A.R. Gallant, D.A. Hsieh, and G.E. Tauchen. “Estimation of stochastic volatility models with diagnostics.” J. Econ. 81 (1997): 159–192. [Google Scholar] [CrossRef]
  15. J.D. Hamilton. “Specification testing in Markov-switching time-series models.” J. Econ. 70 (1996): 127–157. [Google Scholar] [CrossRef]
  16. J. Durbin, and S.J. Koopman. Time Series Analysis by State Space Methods. Oxford, UK: Oxford University Press, 2001. [Google Scholar]

Appendix

A. Derivation of Scores

It is shown here that we can obtain the score function under the null hypothesis using integration by parts, even though the integrated function is degenerate under the null. Hereafter, we denote y ˜ T , h ˜ T by y ˜ , h ˜ , respectively, the multiple integration simply by
d h - - d h 11 d h 21 d h 1 T d h 2 T
and the integration region ( - , ) is suppressed for the sake of notational simplicity, where there is no fear of ambiguity.

A.1. Integration by Parts Formula

We now derive a useful formula for calculating score functions in the following:
- ϖ ( x ) x f ( x ) d x = - - ϖ ( x ) f ( x ) x d x
where ϖ ( x ) is a normal density function and f ( · ) denotes an arbitrary function whose n-th derivative f ( n ) ( x ) for n 0 converges to zero as | x | . We have this formula easily noting that the left-hand side is equal to
- - ϖ ( x ) f ( x ) x d x + ϖ ( ) f ( ) - ϖ ( - ) f ( - ) .
We also have that
- 2 ϖ ( x ) x 2 f ( x ) d x = - ϖ ( x ) 2 f ( x ) x 2 d x
easily by using the formula (18) iteratively.

A.2. Scores of ρ 1 and ρ 2

We obtain the scores of ρ 1 = ϕ 1 and ρ 2 = ϕ 2 - ϕ 1 , where ϕ 1 and ϕ 2 are autocorrelation coefficients of the transition Equation (5). The score of ϕ 1 is expressed as the derivative of the log likelihood L ( ϑ ) evaluated the constrained maximum likelihood estimate of ϑ = ( δ 1 , δ 2 , η , γ , ϕ 1 , ϕ 2 , ω 1 , ω 2 , λ ) , say ϑ ^ , namely
L ( ϑ ^ ) ϕ 1 = 1 f ( y ˜ ) t = 1 T ( A 1 | t + A 2 | t )
where
A 1 | t τ 1 t ϕ 1 τ 1 t - 1 f ( y ˜ , h ˜ ) d h , A 2 | t τ 2 t ϕ 1 τ 2 t - 1 f ( y ˜ , h ˜ ) d h
and τ 1 t = f ( h 1 t | h 1 , t - 1 ) , τ 2 t = f ( h 2 t | h 2 , t - 1 , h 1 t , h 1 , t - 1 ) are transition densities defined in (7), since ϕ 1 is included only in τ 1 t and τ 2 t . Note that, under the null hypothesis, τ 2 t is a degenerate density function and
τ 2 t - 1 τ 2 t ϕ 1 = [ ( h 2 t - h 2 , t - 1 ) - ( h 1 t - h 1 , t - 1 ) ] λ h 1 , t - 1 ω 2
in the integrand of Equation (20) is intractable; the denominator and numerator converges to zero as ω 2 converges to zero under the null, namely ω 2 + 0 , λ = 1 , ϕ 1 = ϕ 2 , and hence h 1 t = h 2 t . We show in the following that this integral can be evaluated by means of integration by parts.
To evaluate Equation (20), we first change the derivative with respect to ϕ 1 to the derivative with respect to h 2 t in the integrand and obtain
A 2 | t = λ h 1 , t - 1 τ 2 t h 2 t τ 2 t - 1 f ( y ˜ , h ˜ ) d h
using the identity
τ 2 t ϕ 1 = λ h 1 , t - 1 τ 2 t h 2 t
which can be shown easily by direct differentiation of Equation (7). Then, using Equation (18), we obtain
A 2 | t = λ h 1 , t - 1 τ 2 t h 2 t τ 2 t - 1 f ( y ˜ , h ˜ ) d h = - λ h 1 , t - 1 ( ν 2 t τ 2 , t + 1 ) h 2 t ν 2 t τ 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h = - λ h 1 , t - 1 ν 2 t h 2 t ν 2 t - 1 f ( y ˜ , h ˜ ) d h - λ h 1 , t - 1 τ 2 , t + 1 h 2 t τ 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h
noting that, from Equation (7), h 2 t is included in ν 2 t and τ 2 , t + 1 as well as τ 2 t . Note that we cannot evaluate the second term of Equation (21) numerically, because τ 2 , t + 1 - 1 ( τ 2 , t + 1 / h 2 t ) is intractable again under the null hypothesis. Then, using
τ 2 , t + 1 h 2 t = - ϕ 2 τ 2 , t + 1 h 2 , t + 1
and Equation (18), we have
A 2 | t = - λ h 1 , t - 1 ν 2 t h 2 t ν 2 t - 1 f ( y ˜ , h ˜ ) d h + λ ϕ 2 h 1 , t - 1 τ 2 , t + 1 h 2 , t + 1 τ 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h .
Using Equation (18) and Equation (22) in Equation (23) iteratively, we finally have
A 2 | t = - λ i = 0 T - t ϕ 2 i h 1 , t - 1 ν 2 , t + i h 2 , t + i ν 2 , t + i - 1 f ( y ˜ , h ˜ ) d h = - λ j = t T ϕ 2 j - t h 1 , t - 1 ν 2 j h 2 j ν 2 j - 1 f ( y ˜ , h ˜ ) d h .
We now express the scores by means of the conditional moments given the observation vector y ˜ under the null hypothesis, noting that f ( y ˜ | h ˜ ) = f ( y ˜ , h ˜ ) / f ( h ˜ ) . Let us denote h t = h 1 t = h 2 t for the sake of notational simplicity, noting that h 1 t = h 2 t under the null hypothesis. Then, under the null hypothesis, we have (Note that λ in Equation (24) equals to one under the null hypothesis.)
L ( ϑ ^ ) ϕ 1 = t = 1 T A ˜ 1 | t - j = t T ϕ ^ 2 j - t A ˜ 2 | j t
where
A ˜ 1 | t τ 1 t ϕ 1 τ 1 t - 1 f ( h ˜ | y ˜ ) d h = E 0 ( h t - 1 u 1 , t - 1 | y ˜ ) ω ^ 1 , A ˜ 2 | j t h 1 , t - 1 ν 2 j h 2 j ν 2 j - 1 f ( h ˜ | y ˜ ) d h = E 0 h t - 1 ξ 2 j | y ˜ - γ ^ E 0 ( h t - 1 ξ 1 j | y ˜ ) η ^ ( 1 - γ ^ 2 )
since
τ 1 t ϕ 1 τ 1 t - 1 = h 1 , t - 1 u 1 , t - 1 ω 1 , ν 2 j h 2 j ν 2 j - 1 = ξ 2 j - γ ξ 1 j η ( 1 - γ 2 ) .
The score of ϕ 2 is expressed as
L ( ϑ ) ϕ 2 = 1 f ( y ˜ ) t = 1 T B t
where
B t τ 2 t ϕ 2 τ 2 t - 1 f ( y ˜ , h ˜ ) d h
since ϕ 2 appears only in τ 2 t as shown in Equation (7). Similarly to the derivative with respect to ϕ 1 , this expression is rewritten as the derivative with respect to h 2 t in the integrand; we have
B t = - h 2 , t - 1 τ 2 t h 2 t τ 2 t - 1 f ( y ˜ , h ˜ ) d h
using
τ 2 t ϕ 2 = - h 2 , t - 1 τ 2 t h 2 t .
For t = T , we have
B T = h 2 , T - 1 ν 2 T h 2 t ν 2 T - 1 f ( y ˜ , h ˜ ) d h
using Equation (22), since h 2 T is included in ν 2 T , as well as τ 2 T , as shown in Equations (4) and (7).
For t < T , we have
B t = h 2 , t - 1 ν 2 t h 2 t ν 2 t - 1 f ( y ˜ , h ˜ ) d h - ϕ 2 h 2 , t - 1 τ 2 , t + 1 h 2 , t + 1 τ 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h = i = 0 T - t ϕ 2 i h 2 , t - 1 ν 2 , t + i h 2 , t + i ν 2 , t + i - 1 f ( y ˜ , h ˜ ) d h = j = t T ϕ 2 j - t h 2 , t - 1 ν 2 j h 2 j ν 2 j - 1 f ( y ˜ , h ˜ ) d h .
Then, under the null hypothesis, we have
L ( ϑ ^ ) ϕ 2 = t = 1 T j = t T ϕ ^ 2 j - t B ˜ j t
where
B ˜ j t = E 0 ( h t - 1 ξ 2 j | y ˜ ) - γ ^ E 0 ( h t - 1 ξ 1 j | y ˜ ) η ^ ( 1 - γ ^ 2 )
and that
A ˜ 2 | j t = B ˜ j t
since ϕ 1 = ϕ 2 and h 1 t = h 2 t = h t . Then, we see that the scores of ρ 1 and ρ 2 are
L ( θ ^ ) ρ 1 = L ( ϑ ^ ) ϕ 1 · ϕ 1 ρ 1 + L ( ϑ ^ ) ϕ 2 · ϕ 2 ρ 1 = t = 1 T A ˜ 1 | t , L ( θ ^ ) ρ 2 = L ( ϑ ^ ) ϕ 1 · ϕ 1 ρ 2 + L ( ϑ ^ ) ϕ 2 · ϕ 2 ρ 2 = L ( ϑ ) ϕ 2 = t = 1 T j = t T ϕ ^ 2 j - t A ˜ 2 | j t
from ϕ 1 = ρ 1 , ϕ 2 = ρ 1 + ρ 2 , and Equation (25).

A.3. Score of λ

The score of λ is expressed as
L ( ϑ ) λ = 1 f ( y ˜ ) t = 1 T C t
where
C t τ 2 t λ τ 2 t - 1 f ( y ˜ , h ˜ ) d h
since λ is only included in the transition density τ 2 t . We express Equation (26) using a derivative with respect to h 2 t in the integrand as
C t = - u 1 , t - 1 τ 2 t h 2 t τ 2 t - 1 f ( y ˜ , h ˜ ) d h
since
τ 2 t λ = - u 1 , t - 1 τ 2 t h 2 t .
Then, for t = T , we have
C T = - u 1 , T - 1 τ 2 T h 2 T τ 2 T - 1 f ( y ˜ , h ˜ ) d h
using Equation (18). For t < T , using Equations (18) and (22) iteratively, we have
C t = u 1 , t - 1 ν 2 t h 2 t ν 2 t - 1 f ( y ˜ , h ˜ ) d h - ϕ 2 u 1 , t - 1 τ 2 , t + 1 h 2 , t + 1 τ 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h = i = 0 T - t ϕ 2 i u 1 , t - 1 τ 2 , t + i h 2 , t + i τ 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h = j = t T ϕ 2 j - t u 1 , t - 1 τ 2 j h 2 j τ 2 j - 1 f ( y ˜ , h ˜ ) d h .
Then, under the null hypothesis, we have
L ( θ ^ ) λ = t = 1 T j = t T ϕ ^ 2 j - t C ˜ j t
where
C ˜ j t u 1 , t - 1 ξ 2 j - γ ^ ξ 1 j η ^ ( 1 - γ ^ 2 ) f ( h ˜ | y ˜ ) d h = E 0 ( u 1 , t - 1 ξ 2 j y ˜ ) - γ ^ E 0 ( u 1 , t - 1 ξ 1 j y ˜ ) η ^ ( 1 - γ ^ 2 ) .

A.4. Score of ω 2

The score of ω 2 is expressed as
L ( θ ) ω 2 = 1 f ( y ˜ ) t = 1 T D t
where
D t τ 2 t ω 2 τ 2 t - 1 f ( y ˜ , h ˜ ) d h .
We reexpress Equation (27) as
D t = 1 2 2 τ 2 t h 2 t 2 τ 2 t - 1 f ( y ˜ , h ˜ ) d h
using
τ 2 t ω 2 = 1 2 2 τ 2 t h 2 t 2
which can be easily shown from Equation (7). For t = T , using the formula (19), we have
D T = E T
where
E T 1 2 2 ν 2 T h 2 T 2 ν 2 T - 1 f ( y ˜ , h ˜ ) d h .
For t < T , we have
D t = E t + F t + H t
where
F t τ 2 , t + 1 h 2 t τ 2 , t + 1 - 1 ν 2 t h 2 t ν 2 t - 1 f ( y ˜ , h ˜ ) d h , H t 1 2 2 τ 2 , t + 1 h 2 t 2 τ 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h
because h 2 t is included in ν 2 t and ν 2 , t + 1 as well as τ 2 t . The second term F t is expressed as
F t = - ϕ 2 τ 2 , t + 1 h 2 , t + 1 τ 2 , t + 1 - 1 ν 2 t h 2 t ν 2 t - 1 f ( y ˜ , h ˜ ) d h
using Equation (22). We then have that
F t = ϕ 2 ν 2 t h 2 t ν 2 t - 1 ν 2 , t + 1 h 2 , t + 1 ν 2 , t + 1 - 1 f ( y ˜ , h ˜ ) d h - ϕ 2 2 τ 2 , t + 2 h 2 , t + 2 τ 2 , t + 2 - 1 ν 2 t h 2 t ν 2 t - 1 f ( y ˜ , h ˜ ) d h
using Equation (18), since h 2 t is included in ν 2 t , τ 2 , t + 1 , as well as τ 2 t . Then, using Equations (19) and (22) iteratively, F t can be expressed as
F t = i = 1 T - t ϕ 2 i G t , t + i
where
G t , t + i ν 2 t h 2 t ν 2 t - 1 ν 2 , t + i h 2 , t + i ν 2 , t + i - 1 f ( y ˜ , h ˜ ) d h
for i = 1 , , T - t . The third term H t can be expressed as
H t = ϕ 2 2 D t + 1
from the identity
2 τ 2 , t + 1 h 2 t 2 = ϕ 2 2 2 τ 2 , t + 1 h 2 , t + 1 2 .
Then, since we have obtained E t , F t , H t , we can express D t as
D t = t = 1 T ϕ 2 2 ( k - t ) E k + k = t T - 1 ϕ 2 2 ( k - t ) l = k + 1 T ϕ 2 l - k G k l .
Thus, under the null hypothesis, we have the score of ω 2 as
L ( θ ^ ) ω 2 = t = 1 T k = 1 t ϕ ^ 2 2 ( k - 1 ) E ˜ t + t = 1 T - 1 k = t T - 1 ϕ ^ 2 2 ( k - t ) l = k + 1 T ϕ ^ 2 l - k G ˜ k l
where
E ˜ t 1 2 2 ν 2 t h 2 t 2 ν 2 t - 1 f ( h ˜ | y ˜ ) d h , G ˜ k l ν 2 t h 2 t ν 2 t - 1 ν 2 j h 2 j ν 2 j - 1 f ( h ˜ | y ˜ ) d h .
Then, we have that
E ˜ t = E 0 ( ξ 2 t 2 | y ˜ ) - 2 γ ^ E 0 ( ξ 2 t ξ 1 t | y ˜ ) + γ ^ 2 E 0 ( ξ 1 t 2 | y ˜ ) η ^ 2 ( 1 - γ ^ 2 ) 2 - 1 η ^ ( 1 - γ ^ 2 ) , G ˜ k l = E 0 ( ξ 2 k ξ 2 l | y ˜ ) - γ ^ E 0 ( ξ 2 k ξ 1 l | y ˜ ) - γ ^ E 0 ( ξ 2 l ξ 1 k | y ˜ ) + γ ^ 2 E 0 ( ξ 1 k ξ 1 l | y ˜ ) η ^ 2 ( 1 - γ ^ 2 ) 2
since
2 ν 2 t h 2 t 2 ν 2 t - 1 = ξ 2 t - γ ξ 1 t η ( 1 - γ 2 ) 2 - 1 η ( 1 - γ 2 ) , ν 2 t h 2 t ν 2 t - 1 = ξ 2 t - γ ξ 1 t η ( 1 - γ 2 ) .

A.5. Scores of δ 1

The score of δ 1 is expressed as
L ( θ ) δ 1 = t = 1 T ν 1 t δ 1 ν 1 t - 1 f ( h ˜ | y ˜ ) d h + ν 2 t δ 1 ν 2 t - 1 f ( h ˜ | y ˜ ) d h
Then, under the null hypothesis, we have
L ( θ ^ ) δ 1 = t = 1 T E 0 ( ξ 1 t | y ˜ ) η ^ - γ ^ E 0 ( ξ 2 t | y ˜ ) - γ ^ E 0 ( ξ 1 t | y ˜ ) η ^ ( 1 - γ ^ 2 )
since
ν 1 t δ 1 ν 1 t - 1 = ξ 1 t η , ν 2 t δ 1 ν 2 t - 1 = - γ ξ 2 t - γ ξ 1 t η ( 1 - γ 2 ) .

A.6. Score of δ 2

Under the null hypothesis, the score of δ 2 is expressed as
L ( θ ^ ) δ 2 = t = 1 T ν 2 t δ 2 ν 2 t - 1 f ( h ˜ | y ˜ ) d h = t = 1 T E 0 ( ξ 2 t | y ˜ ) - γ ^ E 0 ( ξ 1 t | y ˜ ) η ( 1 - γ ^ 2 )
since
ν 2 t δ 2 ν 2 t - 1 = ξ 2 t - γ ξ 1 t η ( 1 - γ 2 ) .

A.7. Score of η

The score of η is expressed as
L ( θ ) η = t = 1 T ν 1 t η ν 1 t - 1 f ( h ˜ | y ˜ ) d h + ν 2 t η ν 2 t - 1 f ( h ˜ | y ˜ ) d h .
Then, under the null hypothesis, we have
L ( θ ^ ) η = t = 1 T 1 2 E 0 ( ξ 2 t 2 | y ˜ ) - 2 γ ^ E 0 ( ξ 1 t ξ 2 t | y ˜ ) + E 0 ( ξ 1 t 2 ) η ^ 2 ( 1 - γ ^ 2 ) - 1 η ^
since
ν 1 t η ν 1 t - 1 = 1 2 ξ 1 t 2 η 2 - 1 η , ν 2 t η ν 2 t - 1 = 1 2 ( ξ 2 t - γ ξ 1 t ) 2 η 2 ( 1 - γ 2 ) - 1 η .

A.8. Score of γ

Under the null hypothesis, the score of γ is expressed as
L ( θ ^ ) γ = t = 1 T ν 2 t γ ^ ν 2 t - 1 f ( h ˜ | y ˜ ) d h = t = 1 T γ ^ 1 - γ ^ 2 - γ ^ E 0 ( ξ 1 t 2 | y ˜ ) + γ ^ E 0 ( ξ 2 t 2 | y ˜ ) - ( 1 + γ ^ 2 ) E 0 ( ξ 1 t ξ 2 t | y ˜ ) η ^ ( 1 - γ ^ 2 ) 2
since
ν 2 t γ ν 2 t - 1 = γ 1 - γ 2 - ( ξ 2 t - γ ξ 1 t ) ( γ ξ 2 t - ξ 1 t ) η ( 1 - γ 2 ) 2 .

A.9. Score of ω 1

Under the null hypothesis, the score of ω 1 is expressed as
L ( θ ^ ) ω 1 = t = 2 T τ 1 t ω 1 τ 1 t - 1 f ( h ˜ | y ˜ ) d h = 1 2 t = 2 T E 0 ( u 1 , t - 1 2 | y ) ω ^ 1 2 - 1 ω ^ 1
since
τ 1 t ω 1 τ 1 t - 1 = 1 2 u 1 , t - 1 2 ω 1 2 - 1 ω 1 .

B. Conditional Moments of State Variables

As we derived in Appendix A, the scores are expressed as the conditional moments of the state variables h 1 t , h 2 t , given the observation y ˜ . We here show that these moments can be easily obtained by applying the completing square formula to the likelihood. Under the null hypothesis, from the transition and measurement density functions (4) and (7),
f ( y ˜ h ˜ ) exp - 1 2 ( y ˜ - δ - Z h ˜ ) Σ - 1 ( y ˜ - δ - Z h ˜ ) , f ( h ˜ ) exp - 1 2 ( Γ h ˜ ) ( Γ h ˜ )
where
δ = 1 T δ 1 δ 2 , Z = I T 1 1 , Σ = I T η η γ η γ η
Γ = 1 - ϕ 1 2 0 0 0 0 - ϕ 1 1 0 0 0 0 - ϕ 1 1 0 0 0 0 0 1 0 0 0 0 - ϕ 1 1
and 1 T is a T × 1 vector of ones. Then, we have
f ( y ˜ , h ˜ ) exp - 1 2 ( y ˜ - δ - Z h ˜ ) Σ - 1 ( y ˜ - δ - Z h ˜ ) + ( Γ h ˜ ) ( Γ h ˜ ) exp - 1 2 ( h ˜ - m ) Ω 1 ( h ˜ - m ) + ( y ˜ - δ ) Ω 2 ( y ˜ - δ )
where m Γ Γ + Z Σ - 1 Z - 1 Z Σ - 1 ( y ˜ - δ ) , Ω 1 Γ Γ + Z Σ - 1 Z , Ω 2 Σ + Z Γ Γ - 1 Z - 1 , applying the completing square formula. Then, we have that
f ( h ˜ | y ˜ ) exp - 1 2 ( h ˜ - m ) Ω 1 ( h ˜ - m ) , f ( y ˜ ) exp - 1 2 ( y ˜ - δ ) Ω 2 ( y ˜ - δ ) .
The conditional expected value and variance-covariance matarix of h ˜ given the value of y ˜ are expressed as E 0 ( h ˜ | y ˜ ) = m , Var 0 ( h ˜ | y ˜ ) = Ω 1 - 1 , respectively. Hereafter, we denote m t = E 0 ( h t | y ˜ ) , c t s = E 0 ( h t h s | y ˜ ) . We also have that
E 0 ( ξ 1 t | y ˜ ) = y 1 t - δ 1 - m t , E 0 ( ξ 2 t | y ˜ ) = y 2 t - δ 2 - m t
E 0 ( u 1 , t - 1 | y ˜ ) = E 0 ( u 2 , t - 1 | y ˜ ) = m t - ϕ 1 m t - 1
from (2) and (5) under the null hypothesis (10).

Share and Cite

MDPI and ACS Style

Chiba, M.; Kobayashi, M. Testing for a Single-Factor Stochastic Volatility in Bivariate Series. J. Risk Financial Manag. 2013, 6, 31-61. https://doi.org/10.3390/jrfm6010031

AMA Style

Chiba M, Kobayashi M. Testing for a Single-Factor Stochastic Volatility in Bivariate Series. Journal of Risk and Financial Management. 2013; 6(1):31-61. https://doi.org/10.3390/jrfm6010031

Chicago/Turabian Style

Chiba, Masaru, and Masahito Kobayashi. 2013. "Testing for a Single-Factor Stochastic Volatility in Bivariate Series" Journal of Risk and Financial Management 6, no. 1: 31-61. https://doi.org/10.3390/jrfm6010031

APA Style

Chiba, M., & Kobayashi, M. (2013). Testing for a Single-Factor Stochastic Volatility in Bivariate Series. Journal of Risk and Financial Management, 6(1), 31-61. https://doi.org/10.3390/jrfm6010031

Article Metrics

Back to TopTop