Next Article in Journal
Exploring Industry-Distress Effects on Loan Recovery: A Double Machine Learning Approach for Quantiles
Previous Article in Journal
Comparing the Conditional Logit Estimates and True Parameters under Preference Heterogeneity: A Simulated Discrete Choice Experiment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Building Multivariate Time-Varying Smooth Transition Correlation GARCH Models, with an Application to the Four Largest Australian Banks

by
Anthony D. Hall
1,
Annastiina Silvennoinen
1 and
Timo Teräsvirta
2,3,*
1
National Centre for Econometric Research (NCER), Queensland University of Technology, Brisbane, QLD 4000, Australia
2
Aarhus BSS, Aarhus University, DK-8210 Aarhus V, Denmark
3
Center for Applied Statistics and Economics (C.A.S.E.), Humboldt-Universität zu Berlin, DE-10178 Berlin, Germany
*
Author to whom correspondence should be addressed.
Econometrics 2023, 11(1), 5; https://doi.org/10.3390/econometrics11010005
Submission received: 15 November 2022 / Revised: 26 January 2023 / Accepted: 27 January 2023 / Published: 6 February 2023

Abstract

:
This paper proposes a methodology for building Multivariate Time-Varying STCC–GARCH models. The novel contributions in this area are the specification tests related to the correlation component, the extension of the general model to allow for additional correlation regimes, and a detailed exposition of the systematic, improved modelling cycle required for such nonlinear models. There is an R-package that includes the steps in the modelling cycle. Simulations demonstrate the robustness of the recommended model building approach. The modelling cycle is illustrated using daily return series for Australia’s four largest banks.

1. Introduction

Recently, Silvennoinen and Teräsvirta (2021) introduced a new multivariate GARCH model called the Multivariate Time-Varying Smooth Transition GARCH model (MTV model). This is a model that explicitly accounts for nonstationarities that are common in daily return series. The authors considered maximum likelihood (ML) estimation of the parameters of the model and, under suitable conditions, proved the consistency and asymptotic normality of the resulting ML estimators.
Before actually estimating an MTV model, however, the model builder has to make a number of data-driven decisions needed for specifying the parametric structure of the model. Further, the estimated structure has to be evaluated by statistical tests to reveal its potential weaknesses. Silvennoinen and Teräsvirta (2021) did not, however, discuss any model building issues, leaving them for further research. The present work is intended to fill this void.
As with many other multivariate GARCH models, the MTV model is based on the decomposition of the conditional covariance matrix Bollerslev (1990), in which the conditional covariance is decomposed to conditional variances and a conditional correlation matrix. In the MTV model, however, it is assumed that the conditional variances can be nonstationary, while a nested special case, (weak) stationarity, is a testable hypothesis.
Likewise, the correlations in this model are time-varying such that its time-varying correlation matrix nests a constant correlation matrix. Due to the parametric structure of this nonlinear correlation matrix, the constancy of correlations has to be tested (and rejected) before fitting a model with time-varying correlations.
Since, as will be explained later, these testing situations, both in the variances and the correlation matrix, are nonstandard, specification of the MTV model is an important issue in building MTV models. Furthermore, the estimated MTV model has to be evaluated before using it, from which it follows that techniques for this part of the model building process have to be examined as well.
In order to illustrate the MTV model building, we consider the Australian banking sector. This is an oligopoly dominated by four banks, commonly called the ‘Big Four’. In early 2020, they represented approximately 19% of the market value of the ASX200 share index and held about 80% of the home loan market in Australia; see Figure 1. Consequently, the banking sector is a major component for many Australian superannuation and other investment funds. As to the Big Four daily returns, their volatility cannot automatically be assumed to be stationary.
Furthermore, the correlations, even when time-varying, cannot a priori be assumed to fluctuate around a constant level, which is one of the assumptions in many popular multivariate GARCH models. Applying the flexible MTV model to these return series is, therefore, an interesting exercise. An in-depth analysis of the Australian banking sector is beyond the scope of this paper; however, modelling the daily returns of the Big Four over a period of almost 30 years serves as a useful example of how our MTV model building techniques work and are applied in practice.
The modelling process is data driven, requiring user input and consists of several steps. For this reason, we developed an R-package to help users build MTV models. The version of R used is 4.1.0. The package, called mtvgarch, includes, among other things, the estimation routines as well as the necessary specification and evaluation tests. The code is maintained in a private GitHub repository and can be obtained upon request.
The plan of the paper is as follows. The MTV model is introduced in Section 2, followed by details of the stages and procedures related to the model building in Section 3. Model specification is considered in Section 4, estimation in Section 5 and evaluation in Section 6. Section 7 is devoted to the illustration of the complete modelling cycle on the Big Four volatilities and correlations. Our conclusions can be found in Section 8. There are also appendices containing material, such as relevant test statistics, simulation results, details of the estimation algorithm, and estimated equations.

2. The MTV Model

The MTV model used in this paper belongs to the family of multivariate GARCH models introduced by Bollerslev (1990). In the original model, the conditional correlations were constant, hence, the name Constant Conditional Correlation (CCC-GARCH) model. This assumption that made the resulting model rather parsimonious was later found to be too restrictive in applications, and time-varying correlations were simultaneously introduced by Engle (2002) (dynamic conditional correlations, DCC) and Tse and Tsui (2002) (varying correlations, VC). In these models, conditional variance components are typically assumed to be stationary, and correlations are assumed, at least implicitly, to fluctuate around a constant level.
In order to consider the MTV model as defined in Silvennoinen and Teräsvirta (2021), we introduce certain notation. The observable stochastic N × 1 vector ε t is decomposed in a customary fashion as
ε t = H t 1 / 2 z t = S t D t P t 1 / 2 ζ t ,
where H t = S t D t P t D t S t is an N × N covariance matrix, and ζ t iid ( 0 , I N ) . We also define z t = P t 1 / 2 ζ t , a vector of independent random variables with E z t = 0 and a positive definite deterministically varying covariance matrix cov( z t ) = P t . The structure of P t will be defined later. The deterministic matrix S t = diag ( g 1 t 1 / 2 , , g N t 1 / 2 ) has positive diagonal elements for all t, and D t = diag ( h 1 t 1 / 2 , , h N t 1 / 2 ) contains the conditional standard deviations of the elements of S t 1 ε t = ( ε 1 t / g 1 t 1 / 2 , , ε N t / g N t 1 / 2 ) . As in Silvennoinen and Teräsvirta (2021) and earlier univariate papers, beginning with Amado and Teräsvirta (2008), and in the multivariate time-varying GARCH article by Amado and Teräsvirta (2014), the diagonal elements of S t 2 are defined as follows:
g i t = g i ( t / T ) = δ i 0 + j = 1 r i δ i j G i j ( t / T , γ i j , c i j ) ,
i = 1 , , N , where δ i 0 > 0 is a known constant, δ i j 0 , j = 1 , , r i , and the (generalised) logistic function
G i j ( t / T , γ i j , c i j ) = ( 1 + exp { γ i j k = 1 K i j ( t / T c i j k ) } ) 1 ,
where γ i j > 0 and c i j = ( c i j 1 , , c i j K i j ) such that c i j 1 c i j K i j . Both γ i j > 0 , c i j 1 c i j K i j , and δ i j 0 , j = 1 , , r i are identification restrictions. Assuming δ i 0 in (2) is known is another one. Furthermore, to prevent exchangeability of the components in (2), restrictions are needed on c i j . As an example, if K i j = 1 for j = 1 , , r i , one can assume (for instance) that c i 11 < < c i r 1 1 .
As discussed in earlier papers, the idea of g i t is to normalise or rescale the observations. Left-multiplying (1) by S t 1 yields
ϕ t = S t 1 ε t = D t z t ,
where each element of ϕ t is assumed to have a standard weakly stationary GARCH representation. In our work, the conditional variances have a GARCH or GJR-GARCH(1,1) structure; see Glosten et al. (1993) for the latter:
h i t = α i 0 + α i 1 ϕ i , t 1 2 + κ i 1 I ( ϕ i , t 1 < 0 ) ϕ i , t 1 2 + β i 1 h i , t 1 ,
where I ( A ) is an indicator function: I ( A ) = 1 when A occurs, zero otherwise. A higher-order structure is possible, although there do not seem to exist applications of the GJR-GARCH model of order greater than one.
The conditional covariance matrix E { ϕ t ϕ t | F t 1 } = D t P t D t . In order to describe the correlation structure, we employ the Double Smooth Transition Conditional Correlation (DSTCC) model by Silvennoinen and Teräsvirta (2009). In that model, assuming that the transition variable is t / T throughout, the time-varying correlation matrix P t is defined as
P t = ( 1 G 2 ( t / T , γ 2 , c 2 ) ) { ( 1 G 1 ( t / T , γ 1 , c 1 ) ) P ( 11 ) + G 1 ( t / T , γ 1 , c 1 ) P ( 21 ) } + G 2 ( t / T , γ 2 , c 2 ) { ( 1 G 1 ( t / T , γ 1 , c 1 ) ) P ( 12 ) + G 1 ( t / T , γ 1 , c 1 ) P ( 22 ) } ,
where P ( i j ) , i , j = 1 , 2 , are four positive definite correlation matrices not equal to each other, and
G i ( t / T , γ i , c i ) = ( 1 + exp { γ i k = 1 K i ( t / T c i k ) } ) 1 , γ i > 0
where c i = ( c i 1 , , c i K i ) ,   c i 1 < < c i K i ,   i = 1 , 2 . This variant of the DSTCC model is called the Time-Varying Correlation (TVC) model to emphasise its deterministic rather than stochastic nature—hence, removing the term ‘Conditional’ from its name. For the Big Four application, we simplify the definition (5) slightly by assuming P ( 12 ) = P ( 22 ) , and therefore (5) becomes
P t = ( 1 G 2 ( t / T , γ 2 , c 2 ) ) { ( 1 G 1 ( t / T , γ 1 , c 1 ) ) P ( 1 ) + G 1 ( t / T , γ 1 , c 1 ) P ( 2 ) } + G 2 ( t / T , γ 2 , c 2 ) P ( 3 ) ,
where re-indexing the matrices highlights the interpretation that there are two transitions over time. One is from P ( 1 ) to P ( 2 ) , and the other one is from a convex combination of these two to P ( 3 ) . Since P ( 1 ) , P ( 2 ) , and P ( 3 ) are positive definite, P t is positive definite as a convex combination of the three matrices. This simplified version of the TVC model is especially useful when modelling correlations that shift from one state to the next as a function of time. To that end, the obvious extension to n such transitions is best expressed as a recursion
P t ( 0 ) = P ( 1 ) P t ( n ) = ( 1 G n ( t / T , γ n , c n ) ) P t ( n 1 ) + G n ( t / T , γ n , c n ) P ( n + 1 ) .
When G 2 ( t / T , γ 2 , c 2 ) 1 and N = 2 , (5) and (7) collapse into the smooth transition correlation GARCH model by Berben and Jansen (2005) or, if the transition variable in G 1 is stochastic and N 2 , into the smooth transition conditional correlation GARCH model of Silvennoinen and Teräsvirta (2005, 2015). An MTV-Conditional Correlation GARCH model with GARCH equations similar to the ones here but differently defined stochastic P t was discussed in Amado and Teräsvirta (2014). It may be noted that Feng (2006) introduced another multivariate Conditional Correlation type GARCH model with deterministically varying correlations. In this model, the variation is described nonparametrically, and the model can be viewed as a generalisation of the univariate model in Feng (2004).

3. The Three Stages of Model Building

The MTV model is rather general and nests many models. To take one example, fitting an MTV model when a nested CCC-GARCH model generates the data leads to inconsistent parameter estimates. For this reason, building adequate MTV models requires care, and a systematic approach is necessary. Selecting a candidate from this family of models is a data-driven process, and statistical inference has to be used to obtain an acceptable model such that it passes the available misspecification tests.
In this work, we follow the classical approach to model building advocated by Box and Jenkins (1970) and later applied to nonlinear models of the conditional mean; see, for example, Teräsvirta et al. (2010, Ch. 16). It has also been applied to building single-equation MTV-GARCH models; see Amado and Teräsvirta (2017) and Amado et al. (2017). The idea is to first specify the model (select a member from the family of MTV models) and, once this has been done, to estimate its parameters. At the evaluation stage, the estimated model is subjected to a battery of misspecification tests. These three stages, specification, estimation, and evaluation, will be considered in the next three sections. The emphasis will be on specification and evaluation as maximum likelihood estimation of the parameters of the MTV model was already considered in Silvennoinen and Teräsvirta (2021).

4. Specification of the MTV Model

4.1. Specification of the Univariate Variance Equations

Specification of the MTV model is begun by specifying the univariate volatility equations. This was first discussed in Amado and Teräsvirta (2017). The idea is to begin with a GARCH(1,1) model by Bollerslev (1986) or the GJR-GARCH model by Glosten et al. (1993) and to test the hypothesis that the multiplicative deterministic component is constant. The single-equation MTV-GARCH model has the following form:
ε i t = z i t h i t 1 / 2 g i t 1 / 2 ,
where z i t iid ( 0 , 1 ) , the conditional variance h i t is defined as in (4) with ϕ i t = ε i t / g i t 1 / 2 , and the deterministic positive-valued function g i t = g i ( t / T ) is defined as in (2) and (3).
Positivity of (2) imposes the following restrictions on δ i j , j = 1 , , r i :
δ i 0 + j = 1 r i δ i j G i j ( r , γ i j , c i j ) > 0
for all r [0,1].
Typically in applications, K i j = 1 , 2 in (3). There are two specification issues, determining r i and choosing K i j , j = 1 , , r i . It is possible that g i ( t / T ) = δ 0 > 0 —that is, g i ( t / T ) is a positive constant. In this case, the MTV-GARCH model collapses into a standard GARCH or GJR-GARCH equation.
Amado and Teräsvirta (2017) solved the problem of choosing r i by first estimating the GARCH model and testing the hypothesis of a constant g i ( t / T ) against the alternative r i = 1 in (2) thereafter using a Lagrange multiplier type test. The test can be viewed as a misspecification test of the estimated GARCH model. If the null hypothesis is rejected, an MTV-GARCH model with a single transition is estimated, and the hypothesis r i = 1 is tested against r i = 2 . Sequential testing continues until the first non-rejection of the null hypothesis.
The number of transitions is determined in this order because of an identification problem: the model with r i + 1 transitions is not identified if the true number of transitions is r i . The shape of the logistic function, controlled by the parameter K i j , can be determined using the sequence of tests familiar from the specification of smooth transition autoregressive (STAR) models; see Teräsvirta (1994) or Teräsvirta et al. (2010, Ch. 16). Details can be found in Amado and Teräsvirta (2017).
More recently, Silvennoinen and Teräsvirta (2016) considered testing the constancy of g i ( t / T ) before estimating the GARCH model—that is, assuming h i t = 1 in (9). The details are laid out in Appendix A.1. This implies that the size of the test is distorted because conditional heteroskedasticity is ignored, so it has to be adjusted by simulation. It turns out that, by doing this, the power of the size-adjusted test considerably improves compared to the case where the test is a standard misspecification test. Reasons for this improvement are discussed in Silvennoinen and Teräsvirta (2016).
A major difficulty with this approach is that, while in simulations, the parameters of the conditional variance component h i t under the null hypothesis are known—in practice, this is not the case. The underlying ‘null’ GARCH process has to be generated artificially. In so doing, special attention is to be placed on the persistence of the (GJR-)GARCH process, measured by α i 1 + κ i 1 / 2 + β i 1 in (4) when g i t 1 . In fact, the asymmetry parameter has no practical importance for the purpose of calibrating the test statistic distribution, and it is therefore sufficient to restrict attention to the standard GARCH process. Other features, such as implied kurtosis or relative sizes of α and β corresponding to a particular level of persistence only have a negligible effect on the performance of the test.
A practical problem is that it is not possible to estimate this measure of persistence when the null hypothesis does not hold—that is, when g i t is not constant over time. How this difficulty is handled has an effect on the power of the test. We study two approaches that are discussed more in detail in Appendix B.1. The first one consists of visually identifying a period of time where there appears to be no change in the overall level of baseline volatility. A standard GARCH(1,1) is estimated over this subperiod. The second approach is to use rolling window variance targeting. This means that the intercept in the GARCH equation is time-varying, and its value at each point in time is calculated such that it matches the unconditional variance obtained from a window around that point in time.
Simulations discussed in Appendix B.1 experiment with the choice of window size. Both of these methods provide GARCH parameter and persistence estimates that are used for calibrating the null distribution of the test statistic and calculating p-values.

4.2. Specification of Time-Varying Correlations

After the MTV-GARCH equations have been specified and estimated assuming the errors are uncorrelated, the next step is to specify the time-varying correlation structure. This is done by sequential testing. First, the constancy of correlation tested against the model with a single transition, i.e., G 2 ( t / T , γ 2 , c 2 ) 1 in (5). The null hypothesis is that the model is a MTV-Constant Correlation GARCH model as in Bollerslev (1990), except that the GARCH equations are MTV-GARCH equations. If this model is rejected, the one-transition model is estimated and tested against (5) or (7). If this specification is also rejected, the alternative with two transitions is estimated. This is repeated until no further evidence for time-variation in the correlations is detected.
As discussed in Silvennoinen and Teräsvirta (2005, 2015), the MTV model with one transition is only identified under the alternative, which invalidates the standard asymptotic inference. The identification problem can be circumvented by approximating the transition function (6) by its Taylor expansion around the null hypothesis, H0: γ 1 = 0 . The form of the expansion depends on the order of the exponent in (6).
The test can be constructed along the lines presented in the appendix of Silvennoinen and Teräsvirta (2005).1 See also Silvennoinen and Teräsvirta (2021). To derive the test statistic, consider the first-order Taylor expansion of (6) around γ 1 = 0 assuming K 1 = 2 . It has the following form:
G 1 ( t / T , γ 1 , c 1 ) = ( 1 + exp { γ 1 k = 1 K 1 ( t / T c 1 k ) } ) 1 = 1 2 + 1 4 ( t / T c 11 ) ( t / T c 12 ) γ 1 + R 2 ( t / T ; γ 1 ) ,
where R 2 ( t / T ; γ 1 ) is the remainder. Using (10), (5) becomes
P t = ( P ( 1 ) P ( 2 ) ) ( 1 2 + γ 1 c 11 c 12 4 ) + P ( 2 ) ( t / T ) ( P ( 1 ) P ( 2 ) ) γ 1 ( c 11 + c 12 ) 4 + ( t / T ) 2 ( P ( 1 ) P ( 2 ) ) γ 1 4 + ( P ( 1 ) P ( 2 ) ) R 2 ( t / T ; γ 1 ) = P ( A 0 ) + ( t / T ) P ( A 1 ) + ( t / T ) 2 P ( A 2 ) + ( P ( 1 ) P ( 2 ) ) R 2 ( t / T ; γ 1 ) ,
where P ( A 0 ) = ( P ( 1 ) P ( 2 ) ) ( 1 / 2 + γ 1 c 11 c 22 / 4 ) + P ( 2 ) , P ( A 1 ) = ( P ( 1 ) P ( 2 ) ) γ 1 ( c 11 + c 12 ) / 4 , P ( A 2 ) = ( P ( 1 ) P ( 2 ) ) γ 1 / 4 , and P ( 1 ) P ( 2 ) . The main diagonals of P ( A 1 ) and P ( A 2 ) consist of zeroes. Setting ρ A = ( ρ A 0 , ρ A 1 , ρ A 2 ) , where ρ A i = vecl ( P ( A i ) ) , i = 0 , 1 , 2 , the new null hypothesis is H0: ρ A 1 = ρ A 2 = 0 N ( N 1 ) / 2 .2
Note that a simpler version of the test assumes K i = 1 and yields a similar approximation although without the term ( t / T ) 2 P ( A 2 ) . The new null in this case is H0: ρ A 1 = 0 N ( N 1 ) / 2 . This version of the test is more powerful than the former in the case that time-variation in the correlations is monotonic. However, and especially with longer time horizons, this may not always be the case, and the square term of the expansion is able to capture at least some nonmonotonic changes.
The details of the ensuing LM-type test statistic for the test of constant correlations is presented in Appendix A.3, and the test for an additional transition in correlations is laid out in Appendix A.4.

5. Estimation of the MTV Model

After specifying the deterministic components of the model, both in GARCH equations and correlations, one can estimate the complete model with conditional heteroskedasticity included. The log-likelihood of the MTV-STCC-GARCH model has the form
ln f ( ζ t | θ ) ( 1 / 2 ) i = 1 N ln g i t ( θ g i ) ( 1 / 2 ) i = 1 N ln h i t ( θ h i ) ( 1 / 2 ) ln | P t ( θ P ) | ( 1 / 2 ) ε t { S t ( θ g ) D t ( θ g , θ h ) P t ( θ P ) D t ( θ g , θ h ) S t ( θ g ) } 1 ε t ,
where the full parameter vector θ = ( θ h , θ g , θ P ) is partitioned according to the relevant functions: the conditional variance in (4) θ h = ( θ h 1 , , θ h N ) with θ h i = ( α i 0 , α i 1 , κ i 1 , β i 1 ) , i = 1 , , N ; the deterministic variance component in (2) and (3) θ g = ( θ g 1 , , θ g N ) with θ g i = ( δ i 1 , γ i 1 , c i 1 , , δ i r i , γ i r i , c i r i ) , i = 1 , , N ; and correlations θ P =   ( vecl P ( 1 ) , , γ 1 , , c 1 , ) , where the number of matrices as well as transition functions (6) are determined by the choice of the model, (5), (7), or (8).
We make the following assumptions; see Silvennoinen and Teräsvirta (2021):
AN1.
In (4), α i 0 > 0 , either α i 1 > 0 and α i 1 + κ i 1 0 or α i 1 0 and α i 1 + κ i 1 > 0 , β i 1 0 , and α i 1 + κ i 1 / 2 + β i 1 < 1 for i = 1 , , N .
AN2.
The parameter subspaces { α i 0 × α i 1 × κ i 1 × β i 1 } , i = 1 , , N , are compact, the whole space Θ h is compact, and the true parameter value θ h 0 is an interior point of Θ h .
AN3.
ζ t iid N ( 0 , I N ) .
AN1 is the necessary and sufficient weak stationarity condition for the ith first-order GJR-GARCH equation. Assumption AN2 is a standard regularity condition required for proving the asymptotic normality of maximum likelihood estimators of θ h i , i = 1 , , N . AN3 (normality) is a strong condition; however, it is needed here for the proofs to go through; see Silvennoinen and Teräsvirta (2021). These assumptions are sufficient for the maximum likelihood estimators of the GARCH parameters in single-equation GARCH models to be consistent and asymptotically normal.
The parameters are estimated in turn: first estimate θ g i to obtain starting-values for the joint estimation of θ g and θ P . This is done assuming h i t ( θ h i ) 1 , i = 1 , , N . Amado and Teräsvirta (2013) showed in the single-equation GJR-GARCH case that, under regularity conditions, the maximum likelihood estimator of θ g i is consistent and asymptotically normal. Silvennoinen and Teräsvirta (2021) generalised this result to MTV models. That means that joint estimation of θ g and θ P by maximum likelihood produces consistent estimates of these parameter vectors.
If θ ^ g and θ ^ P are consistent and Assumptions AN1, AN2, and AN3 hold, then, by Theorem 3.3 of Song et al. (2005), the maximum likelihood estimator of θ h is consistent and asymptotically normal. After estimating θ h , the parameter vectors θ g and θ P are re-estimated. Iteration continues until convergence. Song et al. (2005) showed that the final maximum likelihood estimator of θ is consistent and asymptotically normal. A more detailed description of the maximisation by parts applied to the present situation can be found in Appendix C; see also Silvennoinen and Teräsvirta (2021).

6. Evaluation of the MTV Model

Once the model has been specified and estimated, it has to be evaluated in order to find potential misspecifications. The tests in Section 4.1 were used to guide the choice of the functional form of the deterministic component, and a rejection of the null was seen as evidence of the current model still lacking in its specification. In that sense, the tests in Section 4.1 are seen as both specification and evaluation tests. It is worth reiterating that these specification tests were constructed at the stage when the GARCH part was not yet specified—that is, h t = 1 in (4). However, when the deterministic part passes these tests, and an MTV-GARCH equation is subsequently estimated, there is room for additional checks in terms of model misspecification, beyond the presently final model specification.
The tests in Amado and Teräsvirta (2017) are available for this purpose. They fall into three categories. In the first one, the deterministic component is additively misspecified. In the context of the current MTV-GARCH model, the relevant case is a test for yet another transition in (2). The second test assesses the GARCH equation for additive misspecification. The concern here is the validity of the maximum lags p or q. The final test is the ‘test of no remaining ARCH’, which is based on the idea of a sufficiently well-specified model managing to clear any autocorrelation from the squared standardised residuals. The test that suits each of these situations (or its robustified version to avoid the assumption of normality) is conveniently performed following a set of steps outlined in Appendix A.2.
It is worth stating that the tests here are applied to P ^ t 1 / 2 ε t , one series at a time. Efforts towards completing the tests in the complete N-variate system simultaneously would open up a vast number of permutations of various misspecification options. To manage the task, the recommendation is to focus on the univariate specifications one at a time, even with the acknowledgment of some potential for deviating from the asymptotically exact results. Simulations in Section Appendix B.2 indicate that applying the tests on the pre-filtered data has very little impact on the distributions of the test statistics. While the standard form of the misspecification tests suffers from minor oversizing, this is mostly corrected when using the robust version of the test.
The test for an additional transition in the correlations in Section 4.2 may also be used as an evaluation test. It is based on the completely specified univariate and correlation components, and therefore its role as a misspecification test of a complete model is justified. The number of degrees of freedom in this test quickly becomes large with increasing N. One way of restricting this growth would be to assume that, under the alternative, only the eigenvalues of the correlation matrix are changing over time. The alternative would be a correlation matrix only if all correlations were identical (see Engle and Kelly (2012)); however, an LM test can nevertheless be built on this assumption.
Write the correlation matrix as P t = Q t Λ t Q t , where P t is defined as in (5), Λ t is the diagonal matrix of eigenvalues and Q t contains the corresponding eigenvectors. Simplify this by assuming Q t = Q and approximate Λ t by Ψ t = k = 0 K Ψ k ( t / T ) k . Under the null hypothesis, K = 0 . The resulting test statistic is derived, and its small-sample properties are studied in Kang et al. (2022).

7. Big Four Results

7.1. Main Features of the Australian Banking Sector 1990–2020

In order to provide some background for our empirical results, we shall now draw attention to a number of interesting features of the Australian banking sector between the years 1990 and 2020. In 1990, the Australian government adopted an intervention policy called ‘six pillars’. It covered the four biggest Australian banks commonly referred to as the ‘Big Four’, the Commonwealth Bank of Australia (CBA), Westpac Banking Corporation (WBC), National Australia Bank (NAB) and Australia and New Zealand Banking Group (ANZ), listed in descending order of market share, as well as two insurers (AMP Limited and National Mutual).
This policy stated that further mergers of these institutions would not be accepted. The basic idea was to ensure a competitive banking market. In 1997, the policy became ‘four pillars’ as the insurers were left outside the arrangement. Since its establishment, it has mostly enjoyed the support of the two main political parties, and the proponents of the policy have argued that it has contributed to the stability and strength of the Australian financial sector. The government also had sympathetic policy settings, which allowed the banks to recapitalize in the 1990s and 2000s.
During this period, there were financial losses at Westpac (a $1.6 billion loss in 1992 and close to insolvency), ANZ (poorly executed international expansions) and subsequently with NAB (purchase of the US mortgage originator and servicer Homeside led to $2.2 billion in losses 2002). Although the Big Four were not allowed to merge with each other, larger financial concentration due to mergers with other financial institutions was seen as acceptable: in 2008, Westpac and CBA acquired St. George and BankWest and then the fifth and the sixth largest Australian banks, respectively. The impact can be seen in the right panel of Figure 1 as the Big Four’s share of the ASX200 Financials Index drastically increased.
The recent history of the pillars policy coincides with a few major incidents and changes. These include not only the dot-com boom in the late 1990s and early 2000s and the global financial crisis (GFC) nearly ten years later but also events that have had more localised impacts, such as a number of of regulatory changes (Basel guidelines), the most recent mining boom that started around 2005 and was interrupted by the GFC, and technology-driven market disruptions (non-bank lenders and payment providers). Since the GFC, the banks have enjoyed substantial government support, including a deposit guarantee and, as already noted, have come to dominate the home mortgage market with an 80% market share.
During the first decade of the millennium, a few of the above-mentioned events have positioned the Big Four in an increasingly competitive environment. The stagnation of the housing market and the removal of barriers to changing mortgage providers may also have contributed towards this trend; however, the most notable event was the announcement in 2003 that Basel II was to be implemented in Australia by end of 2007. The idea behind the updated accord was to level inequalities amongst the internationally active banks and to set expectations regarding capital adequacy requirements. The Australian Prudential Regulation Authority, in charge of overseeing the uptake of the accord, worked extensively with numerous Authorised Deposit-Taking Institutions, the industry and other relevant bodies during 2005 to 2007, aiming at ensuring that the adoption of Basel II included all relevant aspects of the implementation process, its goals, and impacts.
Fear of being subjected to a competitive disadvantage relative to their international counterparts, both within international and domestic operations, coupled with an opportunity for a reduced regulatory capital, incentivised the banks to signal early their preference to conform with the accord.3 As a result, Australia was amongst the first nations to have fully implemented the framework on 1 January 2008.

7.2. Modelling the Error Variances

The daily return series for the Big Four used in this paper extend from 2 January 1992 to 31 January 2020. As suggested in the Introduction, these series may not be adequately described by a weakly stationary GARCH or GJR-GARCH model. From the plots in Figure 2, it is seen that the amplitude of clusters varies for all four banks, in particular during and after the financial crisis beginning in 2008. The crisis was preceded by a rather tranquil period between 2003 and 2008. This variation also shows in the autocorrelation functions of the squared returns in Figure 3. In all four cases, the autocorrelations decay very slowly as a function of the lag length, which suggests nonstationarity.
For this reason, modelling the returns has to be initiated by testing the stationarity hypothesis. As discussed in Section 4.1, the slow moving ‘baseline’ volatility is specified first, followed by the inclusion of the GJR-GARCH component. The test statistic from Appendix A.1 is calibrated using methods described in Appendix B.1. In the first one, based on choosing a period during which the amplitude of clusters seems constant, the selected period extends from November 2003 to October 2007. In the other, where a rolling window is moved over the observation period, the window size chosen by simulations is 400, as outlined in Appendix B.1.
To enhance the performance of the test, the entire sample of over 7000 observations is broken into subsections. Once one transition is found and estimated, the test is applied to both before and after this transition to determine if there is another transition on either side. The process is continued until the null of no transition is not rejected. After the deterministic component has been specified and its parameters tentatively estimated, the GJR-GARCH equations are estimated together with the time-varying component to form a complete TV-GJR-GARCH model. The estimated equations are then checked for signs of misspecification using the tests from Amado and Teräsvirta (2017).
Estimated TV-GJR-GARCH equations results appear in Table 1, and the deterministic components are in Appendix D. For comparison, the GJR-GARCH equations without the deterministic component are also reported in this table. In all four cases, the persistence strongly decreases after rescaling the returns with the TV component.
This is also indirectly obvious from the autocorrelations in Figure 4 that are considerably smaller and decay faster than the ones in Figure 3. The main cause for this decrease lies in the coefficient of the lagged conditional variance whose estimate shrinks in the process. Estimates of the asymmetry parameter κ i 1 slightly increase, and thus asymmetry becomes more pronounced when nonstationarity is properly modelled. Table 1 also contains the kurtosis estimates for the two GJR-GARCH processes, obtained using definitions in He and Teräsvirta (1999). It is seen that, in all four cases, rescaling lowers the kurtosis to values close to three.
Figure 5 contains the estimated transitions. There are two conspicuous features in these graphs. One is the downward shift around 2004, which, for WBC, is a long and rather smooth decline. This coincides with the local events discussed in the previous Section. The other is that, for all four banks, the deterministic component remains higher after 2010 than it was before 2008.
For WBC, the deterministic component slowly but steadily declines after the crisis, whereas, for the three others, it remains constant. This can also be seen from the estimated equations behind the figures reported in Appendix D.
Effects of the deterministic component g i t on the GARCH equations also become obvious by comparing the conditional variances from the GJR-GARCH equations in Figure 6 with the TV-GJR-GARCH ones in Figure 7. Clearly, for all four banks, the nonstationarity around 2008–2010 in the former figure is no longer visible in the latter.

7.3. Modelling the Error Correlations

The stability of correlations over time, as discussed in Section 4.2, is tested using the test statistic (A5) in Appendix A.3. The p-value of the test is very close to zero, and thus the null hypothesis is rejected. An TVC model is then estimated with a single monotonic time transition. The adequacy of this specification is tested with the test for an additional transition.
The resulting p-value is 0.467; therefore, the single time-transition is deemed sufficient. The estimation results of the TVC component of the model are presented in Table 2 (see also Figure 8) indicate that the correlations between the standardised residuals are stable, around 0.49–0.60 depending on the bank, from 1992 until the mid-to-late 2006. At that point, the correlations begin their steady increase to the range of 0.78–0.83, which they reach by early 2008. The final correlations are not only large but also remarkably similar.
In addition to considering the complete Big Four system, we repeated the analysis using bivariate models. The outcome was similar: a single transition was sufficient for all pairs. Furthermore, inspection of the pairwise correlation estimates in Table 3 reveals strong similarity between the four-variate and bivariate estimates. This is illustrated in Figure 9. It should be noted that the shift in the correlation of the ANZ–NAB pair is estimated as a step function.
This happens when the speed of transition increases without a bound due to the likelihood function effectively becoming ‘flat’ with respect to that particular parameter. The solution is to fix the speed to a large value and proceed with the estimation of the remaining parameters.
The steady increase of the correlations among the Big Four over time as well as the maintained high correlation state since the GFC could be linked to a few changes in the Australian financial sector. For instance, the time-varying correlation structure may reflect the four pillars policy, the financial concentration, which was further strengthened by the acquisition of the next largest banks by CBA and WBC in 2008. Further contributing factors that may have made the banks look more similar from the investors’ point of view include the regulatory requirements brought by the implementation of Basel II, the easing of restrictions that directly impacted the home mortgage market that the Big Four now dominated and the stagnation of the housing credit market. Furthermore, since the GFC, all four banks have enjoyed substantial government support.
The fact that serious effort has been made to model the volatilities and correlations separately allows for observations on the timing and magnitude of those features without the cross-contamination occurring if covariances were examined instead. It is often noted that correlations increase during turbulent times. In the case of the Big Four, this is not exactly the case. The calm volatility period (from 2003 until late 2007) overlaps with the period of smoothly increasing correlations (16 months leading to early 2008). Furthermore, it is notable that the GFC has a tremendous impact on the volatilities, whereas the correlations have by then settled to their high levels and exhibit no further change.

8. Conclusions

In this paper, a data-driven modelling cycle for building MTV-GARCH models for asset returns was constructed and illustrated with an empirical example. The paper complements Silvennoinen and Teräsvirta (2021), which presented the asymptotic theory for this model, however, as already mentioned in the Introduction, does not contain any discussion on practical model-building issues. All three phases of the cycle: model specification, estimation, and evaluation, were covered. Specification includes testing the constancy of the GARCH equations against multiplicatively time-varying GARCH. This is a nonstandard testing problem as the MTV-GARCH model is not identified when the null hypothesis holds.
Furthermore, constructing the null model requires new techniques due to the fact that the conditional variance is not observed, and two such alternatives for solving this problem are presented. Simulations reported in an online appendix show that the proposed testing procedure had reasonable small-sample properties. In specifying the correlation structure, the constancy of correlations has to be tested, and a relevant test for this nonstandard testing situation was developed.
The GARCH equations and the correlation structure, both nonlinear, were estimated jointly, and the technical details of this process were presented. Misspecification tests for the estimated model were derived to be used for the evaluation of the estimated model. The application to the four main Australian banks, the Big Four, demonstrated the use of the modelling cycle. The GARCH equations were found to be multiplicatively time-varying, and the correlations also changed over time. The estimation results indicated that the amplitude of the volatility clusters declines to a lower level around 2004 and temporarily rises during the GFC.
The (positive) correlations were found to be nonconstant, and they increased to a fairly high level already before the GFC. There is no single reason for this shift that is also established by estimating pairwise models for the six pairs of banks separately. Generally speaking, it may be argued that, whatever the aforementioned events before 2008 have meant to the ability of the banks to compete with each other, from the investors’ viewpoint, they have become increasingly similar. Given their size, these four banks may represent a systemic risk to the Australian financial sector or the Australian economy in general.
Finally, the paper is accompanied by an R-package entitled mtvgarch, which is maintained in a private GitHub repository and contains all the econometric tools necessary for building MTV-GARCH models.
There are appendices contain additional material to the paper. Appendix A provides details of the TVV-model specification, the MTV-GARCH model evaluation, the test of constant correlations, and finally the test for an additional transition in the correlations. The simulations studies in Appendix B explore aspects of the specification and evaluation of the GARCH equations, and the size and sensitivity of the test of constant correlations. Appendix C presents the details of maximisation by parts. The estimated deterministic components of the Four Banks’ transition equations are presented in Appendix D.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/econometrics11010005/s1.

Author Contributions

Conceptualization, A.D.H., A.S. and T.T.; methodology, T.T. and A.S.; software, A.S.; validation, A.S.; formal analysis, A.S. and T.T.; data curation, A.D.H.; writing A.S. and T.T.; visualization, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Supplementary Material that includes data for the application, source code for the estimation and simulations, as well as the MTVGARCH package version used in the production of this paper can be found at https://econ.au.dk/research/researchcentres/creates/research/creates-research-papers/supplementary-downloads/rp-2021-13 (accessed on 26 January 2023).

Acknowledgments

This research was supported by the Center for Research in Econometric Analysis of Time Series (CREATES). An earlier version of this paper with the title ‘Four Australian Banks and the Multivariate Time-Varying Smooth Transition Correlation GARCH model’ appeared as CREATES Research Papers No. 2021–2013. Part of the work was conducted when the third author was visiting School of Economics and Finance of Queensland University of Technology, Brisbane, whose kind hospitality is gratefully acknowledged. We would also like to thank Glen Wade for their work with the R-package. Material from this paper was presented at the 26th Annual Symposium of the Society for Nonlinear Dynamics and Econometrics, Tokyo, March 2018; the workshop ‘Frontiers in Econometrics’, Queensland University of Technology, Brisbane, July 2018; the Quantitative Methods in Finance 2018 Conference (UTS), Sydney, December 2018; the 13th International Conference of the ERCIM WG on Computational and Methodological Statistics, London, December 2020; the International Association of Applied Econometrics Conference, Rotterdam, July 2021; and the Workshop on Financial Econometrics, Örebro University, November 2021. Comments from the participants are gratefully acknowledged. The responsibility for any errors and shortcomings in this work remains ours.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Test Statistics

Appendix A.1. Test Statistic for TVV-Model Specification

In order to specify g t , we not only test constancy but even specify the number of transitions before estimating the GARCH component of the model. Amado and Teräsvirta (2013) showed that maximum likelihood estimators of the corresponding time-varying variance (TVV) model, assuming that there is no conditional heteroskedasticity, are consistent and asymptotically normal. This forms the base for constructing Lagrange multiplier type tests for testing r against r + 1 transitions. For notational simplicity, consider testing one transition against two. Omitting the subscript i for simplicity, the TVV model is (9) with h t = 1 , and
g t = δ 0 + δ 1 G 1 ( t / T , γ 1 , c 1 ) + δ 2 G 2 ( t / T , γ 2 , c 2 ) , γ i > 0 , i = 1 , 2 .
The null hypothesis is γ 2 = 0 —in which case, G 2 ( t / T , γ 2 , c 2 ) 1 / 2 . To circumvent the identification problem (the model with one transition is only identified when the alternative γ 2 > 0 is true), we follow Luukkonen et al. (1988) and approximate the second transition by a third-order Taylor expansion around the null hypothesis. After reparametrisation, this yields
g t = δ 0 * 0 + δ 1 G 1 ( t / T , γ 1 , c 1 ) + ψ 1 t / T + ψ 2 ( t / T ) 2 + ψ 3 ( t / T ) 3 , γ 1 > 0 .
We may call (9) with (A1) the auxiliary TVV model. The parameters ψ i = γ 2 ψ ˜ i , where ψ ˜ i 0 , i = 1 , 2 , 3 . The new null hypothesis in (A1) equals H 0 : ψ 1 = ψ 2 = ψ 3 = 0 . The remainder term of the expansion can be ignored because when we construct a Lagrange multiplier test, the model is only estimated under H 0 (or H 0 ) , and, under this hypothesis, the order of the Taylor expansion equals zero. The remainder is present only under the alternative, and thus ignoring it when H 0 is valid does not affect the asymptotic size of the test. It does make a positive contribution to the power of the test when H 0 does not hold.
Assume (again, for notational simplicity) that K 1 = 1 in (A1), so c 1 = c 1 (a scalar). The log-likelihood for observation t of the auxiliary TVV model equals
t = k ( 1 / 2 ) ln g t ( 1 / 2 ) ε t 2 g t
and the corresponding element of the score is
t θ 1 = 1 2 ( ε t 2 g t 1 ) 1 g t g t θ 1 ,
where θ 1 = ( δ 0 * 0 , δ 1 , γ 1 , c 1 , ψ 1 , ψ 2 , ψ 3 ) . Denoting G 1 ( t / T ) = G 1 ( t / T , γ 1 , c 1 ) , the partial derivative in (A2) is g t / θ 1 = ( g 1 ( t / T ) , τ t ) where
g 1 ( t / T ) = ( 1 , G 1 ( t / T ) , G 1 γ ( t / T ) , G 1 c ( t / T ) G 1 γ ( t / T ) )
with G 1 γ ( t / T ) = G 1 ( t / T ) ( 1 G 1 ( t / T ) ) ( t / T c 1 ) , G 1 c ( t / T ) = γ 1 G 1 ( t / T ) ( 1 G 1 ( t / T ) ) , and τ t = ( t / T , ( t / T ) 2 , ( t / T ) 3 ) . Define the true parameter vector under H 0 as θ 1 0 = ( δ 0 * 0 , δ 1 0 , γ 1 0 , c 1 0 , 0 , 0 , 0 ) . If z t is normally distributed, the corresponding element of the information matrix under H 0 has the form
B t = 1 4 E ( ε t 2 g t 1 ) 2 B 11 t B 12 t B 21 t B 22 t = 1 2 B 11 t B 12 t B 21 t B 22 t ,
where, letting g t 0 = g 0 ( t / T ) = δ 0 * 0 + δ 1 0 G 1 ( t / T , γ 1 0 , c 1 0 ) and denoting G 1 0 ( t / T ) = G 1 ( t / T , γ 1 0 , c 1 0 ) ,
B 11 t = 1 2 ( g 0 ( t / T ) ) 2 g 1 0 ( t / T ) g 1 0 ( t / T ) , B 12 t = 1 2 ( g 0 ( t / T ) ) 2 g 1 0 ( t / T ) τ t
and
B 22 t = 1 2 ( g 0 ( t / T ) ) 2 τ t τ t .
Let
g 1 0 ( r ) = ( 1 , G 1 0 ( r ) , G 1 γ 0 ( r ) , G 1 c 0 ( r ) )
and r = ( r , r 2 , r 3 ) . We state the following lemma:
Lemma A1.
Under the null hypothesis and assuming z t iid N ( 0 , 1 ) , the information matrix
B = lim T 1 2 T t = 1 T B 11 t B 12 t B 21 t B 22 t = 1 2 B 11 B 12 B 21 B 22 ,
where
B 11 = 1 2 0 1 ( g 0 ( r ) ) 2 g 1 0 ( r ) g 1 0 ( r ) d r , B 12 = 1 2 0 1 ( g 0 ( r ) ) 2 g 1 0 ( r ) r d r
and
B 22 = 1 2 0 1 ( g 0 ( r ) ) 2 r r d r .
Proof. 
The ‘sample’ information matrix with T observations equals
1 4 T t = 1 T B t = 1 4 T t = 1 T B 11 t B 12 t B 21 t B 22 t .
Consider the ( 1 , 2 ) element of B 11 ( T ) = ( 1 / T ) t = 1 T B 11 t :
[ B 11 ( T ) ] 12 = ( 1 / T ) t = 1 T ( g 0 ( t / T ) ) 2 G 1 0 ( t / T ) ,
which is an average of T values of the logistic cumulative distribution function. Let [ T r ] = t be the integer closest to t. Then,
( 1 / T ) t = 1 T ( g 0 ( t / T ) ) 2 G 1 0 ( t / T ) = t = 1 T t / T ( t + 1 ) / T ( g * 0 ( [ T r ] / T ) ) 2 G 1 0 ( [ T r ] / T ) d r = 1 / T ( T + 1 ) / T ( g * 0 ( [ T r ] / T ) ) 2 G 1 0 ( [ T r ] / T ) d r 0 1 ( g * 0 ( r ) ) 2 G 1 0 ( r ) d r
as T . The other elements of B 11 = lim T B 11 ( T ) , are derived in a similar fashion. In matrix form,
B 11 = 1 2 0 1 ( g * 0 ( r ) ) 2 g 1 0 ( r ) g 1 0 ( r ) d r .
The blocks B 12 and B 22 are obtained similarly.   □
Since the maximum likelihood estimators of the parameters of the auxiliary TVV model under H 0 are consistent, we may construct the LM test for the hypothesis H 0 : ψ = ( ψ 1 , ψ 2 , ψ 3 ) = 0 . Denoting the relevant block of the score by
s 2 ( θ ^ 1 ) = 1 2 T t = 1 T ( ε t 2 g ^ t 1 ) 1 g ^ t g t ψ ,
where g t / ψ = τ t and
g ^ t = δ ^ 0 + δ ^ 1 ( 1 + exp { γ ^ 1 ( t / T c ^ 1 ) } ) 1 .
Then, assuming z t = ε t / g t 1 / 2 is standard normal under H 0 , the test statistic has the following form:
L M T 1 = ( T / 2 ) s 2 ( θ ^ 1 ) ( B 22 B 21 B 11 1 B 12 ) 1 s 2 ( θ ^ 1 ) ,
where θ ^ 1 = ( δ ^ 0 * , δ ^ 1 , γ ^ 1 , c ^ 1 , 0 , 0 , 0 ) ; see, for example, Godfrey (1988, p. 14). In order to make (A3) operational, the blocks of B are replaced by their consistent counterparts.
When constancy of the error variance is tested against a single transition, g t δ 0 , g 1 ( t / T ) = 1 (scalar), and g t / ψ = τ t as before. Then, B 11 = ( 2 ( δ 0 0 ) 2 ) 1 ,
B 12 = 1 2 ( δ 0 0 ) 2 0 1 r d r and B 22 = 1 2 ( δ 0 0 ) 2 0 1 r r d r .
The test statistic (A3) becomes
L M T 0 = T ( δ 0 0 ) 2 2 δ ^ 0 2 s 2 ( θ ^ 1 ) ( B 22 B 21 B 11 1 B 12 ) 1 s 2 ( θ ^ 1 ) ,
where
s 2 ( θ ^ 1 ) = 1 2 T t = 1 T ( ε t 2 δ ^ 0 1 ) τ t .
When the elements of the covariance matrix are replaced by their consistent estimators in (A4), the ratio ( δ 0 0 ) 2 / δ ^ 0 2 equals unity.
As already mentioned, conditional heteroskedasticity is ignored in setting up the test. For this reason, the test statistic (A3) is likely to be size distorted when applied to financial time series of sufficiently high frequency—that is, when GARCH-type volatility clustering is present. In applications, its size has to be adjusted by calibrating its distribution to reflect the persistence of the GARCH effect present in the data. This is the topic of discussion in Appendix B.1.

Appendix A.2. Test Statistic for MTV-GARCH Model Evaluation

In this section, the evaluation tests of the univariate MTV-GARCH equations are presented in an easy to implement fashion. The full details can be found in Amado and Teräsvirta (2017).
The test statistic is computed based on the following components: ζ ^ t , r 1 t , and r 2 t . ζ ^ t = ε t / h ^ t g ^ t are the residuals, r 1 t contains the derivatives of the functions h t and g t with respect to the parameters that govern the MTV-GARCH model under the null, θ g and θ h :
r 1 t = ( g ^ t 1 g t θ ^ g + h ^ t 1 h t θ ^ g , h ^ t 1 h t θ ^ h ) ,
evaluated at the estimated parameters θ ^ g and θ ^ h . These are recursively calculated, and depend on the prevailing MTV-GARCH model under the null. For example, one could have g t = δ 0 + δ 1 G 1 ( t / T , γ 1 , c 1 ) + δ 2 G 2 ( t / T , γ 2 , ( c 21 , c 22 ) ) and h t = α 0 + α 1 ε t 1 2 / g t 1 + κ 1 I ( ε t 1 < 0 ) ε t 1 2 / g t 1 + β 1 h t 1 . Here, θ g = ( δ 1 , γ 1 , c 1 , δ 2 , γ 2 , c 21 , c 22 ) and θ h = ( α 0 , α 1 , κ 1 , β 1 ) . Then,
g t θ g = ( G 1 , δ 1 G 1 γ 1 , δ 1 G 1 c 1 , G 2 , δ 2 G 2 γ 2 , δ 2 G 2 c 21 , δ 2 G 2 c 22 )
where
G 1 γ 1 = G 1 ( 1 G 1 ) ( t / T c 1 ) G 1 c 1 = G 1 ( 1 G 1 ) γ 1 G 2 γ 2 = G 2 ( 1 G 2 ) ( t / T c 21 ) ( t / T c 22 ) G 2 c 21 = G 2 ( 1 G 2 ) γ 2 ( t / T c 22 ) G 2 c 22 = G 2 ( 1 G 2 ) γ 2 ( t / T c 21 ) .
The GARCH equation derivatives are formed recursively as
h t θ g = g t 1 ( α 1 ε t 1 2 / g t 1 + κ 1 I ( ε t 1 < 0 ) ε t 1 2 / g t 1 ) g t 1 θ g + β 1 h t 1 θ g
and
h t θ h = ( ε t 1 2 / g t 1 , I ( ε t 1 < 0 ) ε t 1 2 / g t 1 , h t 1 ) + β 1 h t 1 θ h
From this example, it should be easy to extend the null model to include more additive deterministic terms and/or have a higher order GARCH equation with or without asymmetric terms.
One extension regarding the deterministic part should be mentioned. It is often convenient to replace the slope parameter γ with e η . In this case, θ g = ( δ 1 , η 1 , c 1 , δ 2 , η 2 , c 21 , c 22 ) , and
G 1 η 1 = G 1 ( 1 G 1 ) e η 1 ( t / T c 1 ) G 1 c 1 = G 1 ( 1 G 1 ) e η 1 G 2 η 2 = G 2 ( 1 G 2 ) e η 2 ( t / T c 21 ) ( t / T c 22 ) G 2 c 21 = G 2 ( 1 G 2 ) e η 2 ( t / T c 22 ) G 2 c 22 = G 2 ( 1 G 2 ) e η 2 ( t / T c 21 ) .
Vector r 2 t contains the derivatives of the misspecified part. Details in the most commonly encountered situations will be given shortly. The number of variables (columns) in r 2 t defines the degrees of freedom in the χ 2 -distribution for the test statistic under the null.
Given the three components, the LM-test is performed as follows:
  • Compute the S S R 0 = t = 1 T ( ζ ^ t 2 1 ) 2 .
  • Regress ζ ^ t 2 1 on ( r 1 t , r 2 t ) , and form the sum of squared residuals S S R 1 .
  • Compute the test statistic L M = T S S R 0 S S R 1 S S R 0 .
The robust version that does not rely on the normality of the error term is formed as follows:
  • Regress r 2 t on r 1 t and obtain residuals w t . When r 2 t has more than one variable, run the regression for each of them separately and, thereby, obtain a set of residuals w t .
  • Regress 1 on ( ζ ^ t 2 1 ) w t and form the sum of squared residuals S S R .
  • Compute the test statistic L M R = T S S R .
The first case seeks to find evidence of misspecification of the determininstic part of the MTV-GARCH model. The conditional variance is of the form
σ t 2 = h t ( g t + f t ) ,
where the additive term f t is zero under the null of the model being correctly specified. The case that we consider here is the one of testing r against r + 1 transitions in the deterministic part. The additive term is linearised and reparameterised, after which, it becomes
f t = δ 0 * + δ 1 * t / T + δ 2 * ( t / T ) 2 + δ 3 * ( t / T ) 3 .
The derivative component for the alternative is then
r 2 t = g ^ t 1 ( 1 , t / T , ( t / T ) 2 , ( t / T ) 3 ) .
The second case addresses misspecification in the GARCH part:
σ t 2 = ( h t + f t ) g t ,
where the additive term f t is again zero under the null. A common scenario is when f t may increase either the ARCH or the GARCH order (but not both). An example of the former is GARCH(1,1) vs. GARCH(2,1), in which case, f t = α 2 ε t 2 2 / g t 2 , and therefore
r 2 t = h ^ t 1 ε t 2 2 / g ^ t 2 .
If the model is a GJR one, and the potential increase in the order of the ARCH term extends to the asymmetric terms as well, then f t = α 2 ε t 2 2 / g t 2 + κ 2 I ( ε t 2 < 0 ) ε t 2 2 / g t 2 , and
r 2 t = h ^ t 1 ( ε t 2 2 / g ^ t 2 , I ( ε t 2 < 0 ) ε t 2 2 / g ^ t 2 ) .
An example of the latter is GARCH(1,1) vs. GARCH(2,1), which leads to f t = β 2 h t 2 , and thus
r 2 t = h ^ t 1 h ^ t 2 .
The third case is the test of no remaining ARCH. This is a test against multiplicative misspecification,
σ t 2 = h t g t f t ,
where f t = 1 under the null. If the alternative is that there is ARCH of order m left unaccounted for, then
r 2 t = ( ζ ^ t 1 2 , , ζ ^ t m 2 ) .

Appendix A.3. Test of Constant Correlations

The log-likelihood of the auxiliary MTV model for observation t assuming K = 2 equals
ln f A ( ζ t | θ ) = ( 1 / 2 ) i = 1 N ln g i t ( 1 / 2 ) i = 1 N ln h i t ( 1 / 2 ) ln | P A t | ( 1 / 2 ) ε t { S t D t P A t D t S t } 1 ε t ,
where
P A t = P ( A 0 ) + ( t / T ) P ( A 1 ) + ( t / T ) 2 P ( A 2 )
and g i t = δ i 0 + δ i 1 G i 1 ( t / T , γ i 1 , c i 1 ) ; only one transition for notational simplicity, and h i t is as in (4). The first sub-block of the score corresponding to the deterministic variance component under H 0 becomes
s t ( θ g i ) = 1 2 ( g i t 1 g i t θ g i + h i t 1 h i t θ g i ) ( 1 e i P ( A 0 ) 1 z t z t e i ) ,
where e i = ( 0 i 1 , 1 , 0 N i ) , i = 1 , , N , and 0 0 is an empty set. The sub-block corresponding to the GARCH parameters under H 0 is
s t ( θ h i ) = 1 2 ( h i t 1 h i t θ h i ) ( 1 e i P ( A 0 ) 1 z t z t e i ) ,
i = 1 , , N . The remaining sub-blocks under H 0 equal
s t ( ρ A j ) = 1 2 vec ( P A t ) ρ A j { vec ( P ( A 0 ) 1 ) ( P ( A 0 ) 1 P ( A 0 ) 1 ) vec ( z t z t ) } = 1 2 ( t / T ) j U { vec ( P ( A 0 ) 1 ) ( P ( A 0 ) 1 P ( A 0 ) 1 ) vec ( z t z t ) } ,
j = 0 , 1 , 2 , where U = vec( P ( A j ) ) / ρ A j consists of zeroes and ones and is identical for all j. The N 2 × N ( N 1 ) / 2 matrix U is a column-wise collection of vectorised indicator matrices that identify the locations of the particular correlation parameters within the matrix P A t . For example, the first correlation parameter in ρ A j is located in positions (2,1) and (1,2) in P A t . An indicator matrix corresponding to this parameter has ones in those positions and zeros elsewhere. This vectorised indicator matrix is then the first column of matrix U and so on. Consequently, the 3 N ( N 1 ) / 2 × N 2 matrix ∂vec ( P A t ) / ρ A equals
vec ( P A t ) ρ A = 1 ( t / T ) ( t / T ) 2 U .
The information matrix for observation t under H 0 is quite similar to, but simpler than, the corresponding one in Silvennoinen and Teräsvirta (2021). In order to give the matrix a proper expression, we need the commutation matrix K , an N 2 × N 2 matrix whose ( i , j ) block equals e j e i —that is, [ K ] i j = e j e i ; see, for example, Lütkepohl (1996, pp. 115–18). Let the superscript 0 indicate that the corresponding entity is evaluated under H 0 (for example, g i t 0 equals g i t | H 0 , and g i t 0 / θ g i equals g i t / θ g i | H 0 ). The matrix is defined in the following lemma.
Lemma A2.
The expectations of the nine blocks of the information matrix at (rescaled) time t / T under H 0 : ρ A 1 = ρ A 2 = 0 N ( N 1 ) / 2 are
B t 0 = E s t ( θ 0 ) s t ( θ 0 ) = E s t ( θ g 0 ) s t ( θ g 0 ) s t ( θ g 0 ) s t ( θ h 0 ) s t ( θ g 0 ) s t ( ρ A ) s t ( θ h 0 ) s t ( θ g 0 ) s t ( θ h 0 ) s t ( θ h 0 ) s t ( θ h 0 ) s t ( ρ A ) s t ( ρ A ) s t ( θ g 0 ) s t ( ρ A ) s t ( θ h 0 ) s t ( ρ A ) s t ( ρ A ) .
The ( i , j ) sub-block of B 11 t = E s t ( θ g 0 ) s t ( θ g 0 ) , i j , equals
[ B 11 t ] i j = E s t ( θ g i 0 ) s t ( θ g j 0 ) = 1 4 ( 1 g i t 0 g i t 0 θ g i + 1 h i t 0 h i t 0 θ g i ) ( 1 g j t 0 g j t 0 θ g j + 1 h j t 0 h j t 0 θ g j ) e i P ( A 0 ) 1 e j e i P ( A 0 ) e j .
When i = j ,
[ B 11 t ] i i = E s t ( θ g i 0 ) s t ( θ g i 0 ) = 1 4 ( 1 g i t 0 g i t 0 θ g i + 1 h i t 0 h i t 0 θ g i ) ( 1 g i t 0 g i t 0 θ g i + 1 h i t 0 h i t 0 θ g i ) ( 1 + e i P ( A 0 ) 1 e i ) .
The ( i , j ) sub-block of B 22 t = E s t ( θ h 0 ) s t ( θ h 0 ) , i j , equals
[ B 22 t ] i j = E s t ( θ h i 0 ) s t ( θ h j 0 ) = 1 4 ( 1 h i t 0 h i t 0 θ h i ) ( 1 h j t 0 h j t 0 θ h j ) e i P ( A 0 ) 1 e j e i P ( A 0 ) e j .
When i = j ,
[ B 22 t ] i i = E s t ( θ h i 0 ) s t ( θ h i 0 ) = 1 4 ( 1 h i t 0 h i t 0 θ h i ) ( 1 h i t 0 h i t 0 θ h i ) ( 1 + e i P ( A 0 ) 1 e i ) .
The ( i , j ) sub-block of B 12 t = E s t ( θ g 0 ) s t ( θ h 0 ) , i j , equals
[ B 12 t ] i j = E s t ( θ g i 0 ) s t ( θ h j 0 ) = 1 4 ( 1 g i t 0 g i t 0 θ g i + 1 h i t 0 h i t 0 θ g i ) ( 1 h j t 0 h j t 0 θ h j ) e i P ( A 0 ) 1 e j e i P ( A 0 ) e j .
When i = j ,
[ B 12 t ] i i = E s t ( θ g i 0 ) s t ( θ h i 0 ) = 1 4 ( 1 g i t 0 g i t 0 θ g i + 1 h i t 0 h i t 0 θ g i ) ( 1 h i t 0 h i t 0 θ h i ) ( 1 + e i P ( A 0 ) 1 e i ) .
Furthermore, the ( i , j ) sub-block of E s t ( θ g 0 ) s t ( ρ A ) equals
[ B 13 t ] i j = E s t ( θ g i 0 ) s t ( ρ A j ) = 1 4 ( t / T ) j ( 1 g i t 0 g i t 0 θ g i + 1 h i t 0 h i t 0 θ g i ) { ( e i e i ) ( P ( A 0 ) 1 I N ) + ( e i e i ) ( I N P ( A 0 ) 1 ) } U ,
i = 1 , , N ; j = 0 , 1 , 2 . The ( i , j ) sub-block of E s t ( θ h 0 ) s t ( ρ A ) equals
[ B 23 t ] i j = E s t ( θ h i 0 ) s t ( ρ A j ) = 1 4 ( t / T ) j ( 1 h i t 0 h i t 0 θ h i ) { ( e i e i ) ( P ( A 0 ) 1 I N ) + ( e i e i ) ( I N P ( A 0 ) 1 ) } U ,
i = 1 , , N ; j = 0 , 1 , 2 . Finally, the ( i , j ) sub-block of the last block is equal to
[ B 33 t ] i j = E s t ( ρ A i ) s t ( ρ A j ) = 1 4 ( t / T ) i + j U M A U ,
i , j = 0 , 1 , 2 , where
M A = P ( A 0 ) 1 P ( A 0 ) 1 + ( P ( A 0 ) 1 I N ) K ( P ( A 0 ) 1 I N ) .
Proof. 
In order to define the test statistic, let B 13 · j be the ( i , j ) blocks of B 13 where i { 1 , , N } —that is,
B 13 · j = [ B 13 ] j 1 , , [ B 13 ] j N , j = 0 , 1 , 2
and define B 23 · j similarly. Partition the matrix B as follows:
B ˜ 11 = B 11 B 12 B 13 · 0 B 12 B 22 B 23 · 0 B 13 · 0 B 23 · 0 [ B 33 ] 00 ,
B ˜ 12 = B 13 · 1 B 13 · 2 [ B 33 ] 01 [ B 33 ] 02
and
B ˜ 33 = [ B 33 ] 11 [ B 33 ] 12 [ B 33 ] 12 [ B 33 ] 22 .
Next, define
x ^ j t = 1 2 ( t T ) j U { vec ( P ( A 0 ) 1 ) ( P ^ ( A 0 ) 1 P ^ ( A 0 ) 1 ) vec ( z ^ t z ^ t ) } ,
j = 1 , 2 , where z ^ t and P ^ ( A 0 ) equal z t and P ( A 0 ) estimated under H 0 , respectively. The test statistic
L M T = T ( 1 T t = 1 T x ^ 1 t , 1 T t = 1 T x ^ 2 t ) { B ˜ 22 B ˜ 12 ( B ˜ 11 0 ) 1 B ˜ 12 } 1 ( 1 T t = 1 T x ^ 1 t , 1 T t = 1 T x ^ 2 t )
has an asymptotic χ 2 -distribution with N ( N 1 ) degrees of freedom when H 0 holds. To make the test statistic operational, the sub-blocks of the information matrix in (A5) have to be replaced by consistent plug-in estimators.

Appendix A.4. Test for an Additional Transition in the Correlations

The test statistic for an additional transition is constructed in the same way as in Appendix A.3, and the blocks related to the volatility components are identical. However, all blocks related to the correlation need modifications to include the parameters governing the time-varying correlation that exists under the null. This includes both the parameters in the correlation matrices under the null, and their corresponding transition parameters.
Let us define x h i t = h i t 1 h i t θ h i , x g i t = g i t 1 g i t θ g i + h i t 1 h i t θ g i . Let us also partition the linearised correlation model as P A t = P A t 0 + t / T P ( A 1 ) + ( t / T ) 2 P ( A 2 ) , where P A t 0 contains the time-varying correlation model under the null. When testing L transitions against L + 1 transitions, P A t 0 contains L + 1 correlation matrices P ( 1 ) , , P ( L + 1 ) and L transition functions G l ( t / T , γ l , c l ) (here, we assume K l = 1 for simplicity), l = 1 , , L . The information matrix is approximated by its consistent estimator
B ^ = T 1 t = 1 T E t 1 [ s t ( θ 0 ) s t ( θ 0 ) ] ,
where θ 0 = ( θ g , θ h , θ G , θ ρ A 0 , θ ρ A 1 ) , where θ G contains the transition parameters from the L transitions that are present under the null, θ ρ A 0 = ( ρ ( 1 ) , , ρ ( L + 1 ) ) and θ ρ A 1 = ( ρ ( A 1 ) , ρ ( A 2 ) ) . From here on, the expressions are evaluated at the true parameter values under the null (we omit the additional superscripts of 0 to keep the notation simple).
With this notation, the ( i , j ) sub-block of B ^ 11 , i j , equals
[ B ^ 11 ] i j = 1 4 T t = 1 T x g i t x g j t e i P ( A t 0 ) 1 e j e i P ( A t 0 ) e j .
When i = j ,
[ B ^ 11 ] i i = 1 4 T t = 1 T x g i t x g i t ( 1 + e i P ( A t 0 ) 1 e i ) .
Similarly, the ( i , j ) sub-block of B ^ 22 , i j , is equal to
[ B ^ 22 ] i j = 1 4 T t = 1 T x h i t x h j t e i P ( A t 0 ) 1 e j e i P ( A t 0 ) e j .
When i = j ,
[ B ^ 22 ] i i = 1 4 T t = 1 T x h i t x h i t ( 1 + e i P ( A t 0 ) 1 e i ) .
The ( i , j ) sub-block of B ^ 12 , i j , equals
[ B ^ 12 ] i j = 1 4 T t = 1 T x g i t x h j t e i P ( A t 0 ) 1 e j e i P ( A t 0 ) e j .
When i = j ,
[ B ^ 12 ] i i = 1 4 T t = 1 T x g i t x h i t ( 1 + e i P ( A t 0 ) 1 e i ) .
The next blocks deal with the transition parameters. Define x G t = vec P A t θ G . The lth block of vec P A t θ G is
i = l + 1 L ( 1 G i ) G i θ G i vec ( P ( l + 1 ) P t ( l 1 ) )
using the recursion in (8). The block B ^ 33 is equal to
B ^ 33 = 1 4 T t = 1 T x G t M A x G t .
The ith sub-block of B ^ 13 equals
[ B ^ 13 ] i = 1 4 T t = 1 T x g i t ( e i P ( A t 0 ) 1 e i + e i e i P ( A t 0 ) 1 ) x G t
and the ith sub-block of B ^ 23 equals
[ B ^ 23 ] i = 1 4 T t = 1 T x h i t ( e i P ( A t 0 ) 1 e i + e i e i P ( A t 0 ) 1 ) x G t .
Next, we will consider the blocks related to the correlations. The matrix ∂vec ( P A t ) / θ ρ A 0 equals
vec ( P A t ) θ ρ A 0 = v s . U = l = 1 L ( 1 G l t ) G 1 t l = 2 L ( 1 G l t ) G 2 t l = 3 L ( 1 G l t ) G l 1 , t ( 1 G L t ) G L t U
and the matrix ∂vec ( P A t ) / θ ρ A 1 equals
vec ( P A t ) ρ A 1 = t / T ( t / T ) 2 U .
The block B ^ 44 is equal to
B ^ 44 = 1 4 T t = 1 T v t v t U M A U
and the ( i , j ) sub-block of B ^ 55 is equal to
[ B ^ 55 ] i j = 1 4 T t = 1 T ( t / T ) i + j U M A U
for i , j = 1 , 2 . The ith sub-block of B ^ 14 is equal to
[ B ^ 14 ] i = 1 4 T t = 1 T x g i t ( e i P ( A t 0 ) 1 e i + e i e i P ( A t 0 ) 1 ) ( v t U )
and the ( i , j ) sub-block of B ^ 15 is equal to
[ B ^ 15 ] i j = 1 4 T t = 1 T ( t / T ) j x g i t ( e i P ( A t 0 ) 1 e i + e i e i P ( A t 0 ) 1 ) U ,
j = 1 , 2 . The corresponding sub-blocks of B ^ 24 and B ^ 25 are
[ B ^ 24 ] i = 1 4 T t = 1 T x h i t ( e i P ( A t 0 ) 1 e i + e i e i P ( A t 0 ) 1 ) ( v t U )
and the ( i , j ) sub-block of B ^ 15 is equal to
[ B ^ 25 ] i j = 1 4 T t = 1 T ( t / T ) j x h i t ( e i P ( A t 0 ) 1 e i + e i e i P ( A t 0 ) 1 ) U ,
j = 1 , 2 . The block B ^ 34 equals
B ^ 34 = 1 4 T t = 1 T x G t M A ( v U ) .
The ith sub-block of B ^ 35 equals
[ B ^ 35 ] i = 1 4 T t = 1 T ( t / T ) i x G t M A U ,
i = 1 , 2 Finally, the ith sub-block of B ^ 45 is equal to
[ B ^ 45 ] i = 1 4 T t = 1 T ( t / T ) i ( v s . U ) M A U .
Next, define
x ^ j t = 1 2 ( t T ) j U { vec ( P ( A 0 t ) 1 ) ( P ^ ( A 0 t ) 1 P ^ ( A 0 t ) 1 ) vec ( z ^ t z ^ t ) } ,
j = 1 , 2 , where z ^ t and P ^ ( A 0 t ) equal z t and P ( A 0 t ) estimated under H 0 , respectively. The test statistic
L M T = T ( 1 T t = 1 T x ^ 1 t , 1 T t = 1 T x ^ 2 t ) [ B ^ 1 ] S W ( 1 T t = 1 T x ^ 1 t , 1 T t = 1 T x ^ 2 t ) ,
where [ B ^ 1 ] S W is the N ( N 1 ) × N ( N 1 ) block in the south-west corner of the inverse of B ^ . As the matrix B ^ can have a large dimension, its inverse could be obtained by using block inversion methods, perhaps applying them recursively. The test statistic has an asymptotic χ 2 -distribution with N ( N 1 ) degrees of freedom when H 0 holds.

Appendix B. Simulations of Test Statistics

Appendix B.1. Tests of GARCH Equations

The test for slow moving baseline volatility has a statistic whose distribution is sensitive to the high frequency, GARCH, volatility. For this reason, one cannot use the asymptotic distribution, rather the distribution must be generated via simulation. Further, Silvennoinen and Teräsvirta (2016) showed that the size of the test was distorted if the GARCH parameterisation deviates from the true one. For this reason, a few alternative approaches to estimate the GARCH parameters, and especially the persistence, have been investigated. It should be noted that estimating GARCH without taking the nonstationarity into account will yield overestimated persistence, thereby, impacting the null distribution of the test statistic and thus rendering the test outcomes unreliable. These estimates are given in Table A1.
The baseline volatility may be very different in different series. Therefore, one should not ignore visual inspection of the returns nor rely on general rules of thumb. If there are sufficiently long sections of data where the general level of volatility remains constant, it is advisable to estimate the GARCH parameters over such subsample. In the present case, there are a couple of relatively constant volatility sections—for example, one from November 2003 until October 2007.
The parameter estimates for that calm subperiod are in Table A1. Comparison with the estimates from the entire period GARCH model makes it clear that the neglected nonstationarity has biased the estimates, resulting in high persistence and kurtosis. As the data set has a sufficiently long span of GARCH-type clustering without (visually) significant movement in the general baseline level, relevant estimates are obtained by using that subsample only.
Another approach consists of estimating the GARCH equation over a rolling window such that the intercept is time-varying, targeting the unconditional volatility over each window, while the other parameters are assumed to be constant over the entire sample period and estimated in the usual way. The choice of the window length should consider the general recommendations regarding the sample size when attempting GARCH estimation. Too long a window will be impacted by the slowly changing baseline volatility level, whereas too short a window will yield very uncertain GARCH estimates.
To investigate the properties of this approach, we ran a simulation experiment with a few different baseline volatilities. The window widths varied from 250 to 1000 observations. Figure A1, Figure A2 and Figure A3 depict the distributions of the GARCH estimates and the derived persistence and kurtosis measures, as explained in He and Teräsvirta (1999), for a selection of baseline volatilities and window widths. Based on these experiments, we concluded that a window width of 400 observations yields sufficiently robust results for our application.
The resulting GARCH estimates are reported in Table A1, and they are quite similar to the ones obtained for the aforementioned calm period. This can be interpreted as support for the rolling window method, particularly in situations where visual inspection of data does not reveal a sufficiently long period of constant unconditional volatility.
Overall, it is clear that using simply the GARCH estimates from the entire sample to calibrate the null distribution of the test statistic for the specification of the deterministic component of the volatility is not recommended. For comparison, Table A1 also reports the GARCH estimates from a TV-GARCH model where the TV specification has been completed. The estimated persistence is higher than the ones obtained from the calm period or rolling window variance targeting method, however, as discussed in Silvennoinen and Teräsvirta (2016), underestimation of persistence has a less severe impact on the performance of the TV specification test than does overestimation.
Table A1. Specification stage for the deterministic component in volatilities of each of the four banks. α ˜ and β ˜ are the initial estimates used for calibrating the test statistic distribution. The rolling window method allows the GARCH intercept to adjust to target the unconditional variance in a window of size 400. The ‘calm period’ selects the continuous period from November 2003 to October 2007, which has very little visible variation in the baseline volatility. For comparison, the GARCH estimates from the entire sample period are reported along with the final estimates from the TV-GARCH model.
Table A1. Specification stage for the deterministic component in volatilities of each of the four banks. α ˜ and β ˜ are the initial estimates used for calibrating the test statistic distribution. The rolling window method allows the GARCH intercept to adjust to target the unconditional variance in a window of size 400. The ‘calm period’ selects the continuous period from November 2003 to October 2007, which has very little visible variation in the baseline volatility. For comparison, the GARCH estimates from the entire sample period are reported along with the final estimates from the TV-GARCH model.
α ˜ β ˜ PersistenceKurtosis
Rolling window 400ANZ0.0900.8360.9263.38
CBA0.0870.8500.9373.43
NAB0.0950.8170.9123.36
WBC0.0850.8580.9433.45
Calm periodANZ0.0730.8520.9253.24
CBA0.0810.8420.9233.29
NAB0.0660.8290.8963.14
WBC0.0910.8060.8973.28
Entire period GARCH onlyANZ0.0650.9270.9926.40
CBA0.0890.8900.9794.83
NAB0.1040.8670.9714.85
WBC0.0750.9110.9865.08
Entire period TV-GARCHANZ0.0780.8800.9573.50
CBA0.0910.8600.9503.61
NAB0.1070.8250.9313.62
WBC0.0840.8780.9623.70
Figure A1. Simulated distributions of GARCH estimates and implied persistence and kurtosis measures for a selection of window widths. The baseline g t has a single transition. The dotted vertical lines indicate the true values of the parameters α , β , persistence, and kurtosis.
Figure A1. Simulated distributions of GARCH estimates and implied persistence and kurtosis measures for a selection of window widths. The baseline g t has a single transition. The dotted vertical lines indicate the true values of the parameters α , β , persistence, and kurtosis.
Econometrics 11 00005 g0a1
Figure A2. Simulated distributions of GARCH estimates and implied persistence and kurtosis measures for a selection of window widths. The baseline g t has an asymmetric double transition. The dotted vertical lines indicate the true values of the parameters α , β , persistence, and kurtosis.
Figure A2. Simulated distributions of GARCH estimates and implied persistence and kurtosis measures for a selection of window widths. The baseline g t has an asymmetric double transition. The dotted vertical lines indicate the true values of the parameters α , β , persistence, and kurtosis.
Econometrics 11 00005 g0a2
Figure A3. Simulated distributions of GARCH estimates and implied persistence and kurtosis measures for a selection of window widths. The baseline g t has two double transitions. The dotted vertical lines indicate the true values of the parameters α , β , persistence, and kurtosis.
Figure A3. Simulated distributions of GARCH estimates and implied persistence and kurtosis measures for a selection of window widths. The baseline g t has two double transitions. The dotted vertical lines indicate the true values of the parameters α , β , persistence, and kurtosis.
Econometrics 11 00005 g0a3

Appendix B.2. Evaluation Tests of GARCH Equations

The fact that the evaluation tests discussed in Appendix A.2 are applied to the pre-filtered data P t ^ 1 / 2 ε t is known to potentially alter the distribution of the test statistic. In this section, we present simulation results that show that the size of the tests remains practically unchanged, rendering the tests applicable in the proposed way.
The simulation uses 2000 observations on a bivariate TVGARCH model parametrised as h t = 0.10 + 0.05 ε t 1 2 / g t 1 + 0.85 h t 1 , g t = 1 + 3 ( 1 + exp { e 3 ( t / T 0.5 ) } ) 1 . These are coupled with a CCC model with ρ = 0.5 , and then with an STCC model parametrised as ρ ( 1 ) = 0.3 , ρ ( 2 ) = 0.7 , G t = ( 1 + exp { e 2.5 ( t / T 0.5 ) } ) 1 . The noise terms are iid standard normal. Two estimation procedures were used, a two-step and a multi-step one.
First step
The individual TVGARCH models are estimated, assuming the series are uncorrelated.
Second step
Estimate the correlation model conditional on the volatility model estimates from the previous step. Then, estimate the TVGARCH models conditional on the correlation estimates.
The misspecification tests are then calculated using the TVGARCH estimates from the second step, and the data is pre-filtered with the correlation estimates from the second step. The multi-step continues repeating the procedure of the second step until no further improvements are achieved.
Table A2. Size simulation for the three types of misspecification tests in Amado and Teräsvirta (2017). 2000 replications. T = 2000 , N = 2 . MS1: g t additively misspecified, alternative linearised with a first-order term only; MS2-a: GARCH(1,1) vs. GARCH(1,2); MS2-b: GARCH(1,1) vs. GARCH(2,1); MS3: test for remaining ARCH, lag 1.
Table A2. Size simulation for the three types of misspecification tests in Amado and Teräsvirta (2017). 2000 replications. T = 2000 , N = 2 . MS1: g t additively misspecified, alternative linearised with a first-order term only; MS2-a: GARCH(1,1) vs. GARCH(1,2); MS2-b: GARCH(1,1) vs. GARCH(2,1); MS3: test for remaining ARCH, lag 1.
StandardRobust
10%5%1%10%5%1%
CCC two-stepMS10.1460.0850.0200.1320.0740.016
MS2-a0.1220.0640.0120.1010.0480.013
MS2-b0.1430.0800.0170.1080.0510.008
MS30.1250.0610.0100.1040.0540.010
STCC two-stepMS10.1340.0740.0230.1210.0550.015
MS2-a0.1230.0590.0150.1010.0450.013
MS2-b0.1220.0620.0190.0870.0440.010
MS30.1150.0580.0150.1000.0500.011
CCC multi-stepMS10.1450.0830.0220.1330.0730.014
MS2-a0.1160.0620.0150.0970.0520.009
MS2-b0.1330.0690.0180.1000.0460.010
MS30.1200.0620.0160.1070.0600.014
STCC multi-stepMS10.1470.0840.0230.1350.0680.012
MS2-a0.1300.0590.0110.1030.0460.006
MS2-b0.1200.0670.0160.0900.0390.005
MS30.1120.0550.0120.1040.0470.009
From Table A2, it is evident that the standard form of the tests is slightly oversized. The robust version of the tests, on the other hand, seems to behave well, and there is no need for any adjustments of the test statistics or their distributions. Therefore, the procedure of removing the correlations between the series prior to applying the evaluation tests can be recommended.

Appendix B.3. Tests of Correlations

The simulation experiment investigates the size of the test in an environment where the multivariate model is correctly specified. The number of data series considered in the system is N = 2 , 5 , 10 , 20 . The length varies from T = 25 for the bivariate systems, which is relevant for time series systems in macro applications, up to T = 1000 , which, in turn, is considered to be a fairly small sample size for high frequency returns data. The length of the time series places a constraint on the dimension of the model—that is, the parametric alternative is only feasible if the number of parameters remains comfortably below the amount of available data points. We simulated the test by both assuming that D t I N and that there is conditional heteroskedasticity in the model: D t I N .
When D t I N , we found that the results were fairly independent of the structure of the correlations. We used both equicorrelation and Toeplitz matrices in our simulations, and the results remained the same. Table A3 contains the results of a simulation in which D t I N , and the N × N correlation matrix P = [ ρ i j ] is an equicorrelated one with weak ( ρ = 1 / 3 ) and moderately strong ( ρ = 2 / 3 ) correlation. The table also reports the results from using a Toeplitz correlation matrix such that [ ρ i j ] = ρ | i j | ,   i , j = 1 , , N with ρ = 0.5 representing moderate to weak correlation and ρ = 0.9 representing strong to moderate correlation. It is seen that the empirical size of the test is rather close to the nominal one already when N = 2 and T = 100 . The size holds up across the various correlation patterns.
Table A3. Size-study: Test of constant correlations. Data are generated as an MTV-CCC with an equicorrelation coefficient of 0.33 (CEC33) and 0.67 (CEC67) and a Toeplitz structure with a correlation coefficient of 0.5 (CTC50) and 0.9 (CTC90). Tests are based on the first-order polynomial approximation. A total of 5000 replications.
Table A3. Size-study: Test of constant correlations. Data are generated as an MTV-CCC with an equicorrelation coefficient of 0.33 (CEC33) and 0.67 (CEC67) and a Toeplitz structure with a correlation coefficient of 0.5 (CTC50) and 0.9 (CTC90). Tests are based on the first-order polynomial approximation. A total of 5000 replications.
CEC33CEC67CTC50CTC90
NT1%5%10%1%5%10%1%5%10%1%5%10%
2250.0230.0760.1320.0220.0690.1280.0240.0740.1300.0220.0700.126
500.0150.0630.1160.0160.0640.1150.0160.0640.1150.0150.0620.109
1000.0110.0560.1040.0100.0540.1020.0110.0560.1030.0100.0510.101
2500.0120.0550.1080.0100.0540.1070.0110.0550.1060.0090.0530.108
5000.0100.0510.0970.0090.0490.0970.0100.0500.0960.0090.0500.094
10000.0100.0480.0990.0100.0480.0950.0100.0460.0970.0100.0490.092
51000.0110.0540.1120.0110.0530.1100.0110.0560.1120.0110.0530.111
2500.0140.0540.0990.0120.0510.0990.0130.0530.1000.0120.0510.101
5000.0100.0500.1040.0100.0530.1060.0090.0520.1010.0100.0540.105
10000.0100.0560.1020.0100.0520.1030.0090.0530.1000.0080.0530.103
102500.0130.0550.1120.0130.0570.1120.0130.0570.1100.0120.0540.115
5000.0090.0490.1010.0100.0490.1040.0080.0530.1030.0100.0500.103
10000.0110.0520.1020.0110.0540.1050.0110.0530.0990.0120.0560.103
2010000.0120.0560.1060.0120.0570.1060.0130.0560.1030.0120.0560.107
We next turn to the case D t I N . Table A4 and Table A5 contain results of size simulations where the sensitivity of the test is examined against combinations for the GARCH persistence and kurtosis as well as a selection of strengths of correlations (the equicorrelated and Toepliz ones described above). The test is generally well-sized.
Table A4. Size-study: Test of constant correlations. Data are generated as an MTV-GARCH-CEC with persistence of 0.95 and 0.97, kurtosis of 4 and 6, and an equicorrelation coefficient of 0.33 and 0.67. Tests are based on the first-order polynomial approximation. A total of 2500 replications.
Table A4. Size-study: Test of constant correlations. Data are generated as an MTV-GARCH-CEC with persistence of 0.95 and 0.97, kurtosis of 4 and 6, and an equicorrelation coefficient of 0.33 and 0.67. Tests are based on the first-order polynomial approximation. A total of 2500 replications.
   CEC33CEC67
   kurtosis = 4kurtosis = 6kurtosis = 4kurtosis = 6
PersistenceNT1%5%10%1%5%10%1%5%10%1%5%10%
0.9525000.0120.0560.1080.0160.0560.1060.0160.0700.1220.0160.0920.122
210000.0090.0440.1030.0090.0420.0970.0110.0450.0930.0090.0440.097
220000.0080.0420.0940.0070.0420.0900.0100.0520.0990.0090.0460.092
55000.0060.0620.1180.0060.0700.1140.0180.0760.1400.0180.0820.146
510000.0160.0600.1190.0160.0610.1120.0160.0590.1150.0180.0600.112
520000.0100.0580.1080.0080.0510.1020.0160.0600.1160.0100.0520.098
105000.0160.0580.1180.0200.0640.1140.0200.0680.1160.0240.0800.128
1010000.0180.0530.1040.0150.0510.1010.0140.0610.1110.0170.0630.110
1020000.0140.0720.1260.0120.0600.1120.0180.0820.1420.0130.0620.118
0.9725000.0100.0560.1140.0120.0540.1180.0200.0720.1140.0140.0680.120
210000.0110.0430.1020.0110.0440.1030.0120.0470.1070.0130.0480.103
220000.0090.0460.0940.0070.0420.0890.0100.0560.1080.0120.0500.093
55000.0040.0660.1240.0120.0560.1040.0120.0880.1520.0180.0860.164
510000.0150.0630.1130.0140.0670.1140.0180.0630.1210.0190.0600.125
520000.0100.0600.1100.0080.0500.1000.0150.0600.1180.0120.0500.101
105000.0120.0620.1080.0160.0700.1120.0160.0720.1120.0220.0860.148
1010000.0160.0530.1000.0150.0560.1070.0150.0630.1130.0180.0570.110
1020000.0150.0740.1320.0140.0580.1080.0160.0880.1420.0100.0630.112
Table A5. Size-study: Test of constant correlations. Data are generated as an MTV-GARCH-CTC with persistence of 0.95 and 0.97, kurtosis of 4 and 6, and a correlation matrix with a Toeplitz structure with a correlation coefficient of 0.5 and 0.9. Tests are based on the first-order polynomial approximation. A total of 2500 replications.
Table A5. Size-study: Test of constant correlations. Data are generated as an MTV-GARCH-CTC with persistence of 0.95 and 0.97, kurtosis of 4 and 6, and a correlation matrix with a Toeplitz structure with a correlation coefficient of 0.5 and 0.9. Tests are based on the first-order polynomial approximation. A total of 2500 replications.
   CTC50CTC90
   kurtosis = 4kurtosis = 6kurtosis = 4kurtosis = 6
PersistenceNT1%5%10%1%5%10%1%5%10%1%5%10%
0.9525000.0100.0640.1020.0100.0700.1060.0180.0940.1360.0260.0880.146
210000.0090.0410.0970.0110.0420.1030.0140.0530.0960.0200.0620.104
220000.0080.0440.0960.0090.0440.0900.0170.0660.1200.0140.0480.098
55000.0060.0620.1180.0100.0580.1140.0200.1200.2120.0500.1340.210
510000.0140.0600.1120.0180.0640.1130.0270.0760.1340.0340.0930.144
520000.0110.0570.1100.0080.0520.1050.0200.0750.1420.0180.0580.110
105000.0120.0700.1200.0160.0800.1280.0400.1140.1720.0780.1500.230
1010000.0120.0490.1000.0130.0510.1020.0190.0780.1270.0320.0890.147
1020000.0190.0720.1270.0140.0590.1110.0330.1100.1780.0180.0770.140
0.9725000.0140.0660.1140.0180.0680.1160.0160.0820.1340.0300.1040.164
210000.0090.0440.1010.0080.0420.0990.0160.0510.1120.0220.0630.119
220000.0100.0500.1020.0090.0440.0920.0240.0700.1200.0150.0520.100
55000.0140.0560.1280.0080.0740.1300.0240.1340.2080.0520.1600.256
510000.0130.0590.1120.0160.0660.1230.0220.0820.1570.0370.1020.164
520000.0140.0620.1120.0100.0520.1010.0280.0860.1450.0200.0660.116
105000.0180.0800.1280.0220.0880.1300.0400.1140.1720.1000.1880.278
1010000.0120.0540.1050.0130.0540.1070.0190.0780.1270.0300.1040.181
1020000.0160.0720.1320.0160.0620.1100.0330.1100.1780.0260.0890.150
However, an interesting aspect is that there is slight oversizing when kurtosis decreases (which means shifting the relative weight from α to β in the GARCH equation, while keeping the persistence constant). A change in persistence does not seem to affect the size of the test. As the dimension of the system increases, the test does not perform equally well. Increasing the sample size does not seem to be able to counteract this (the simulations use T = 500 , 1000 , 2000 ).
In yet another simulation (results not reported here), we considered the effects of misspecifying the conditional heteroskedasticity on the correlation test. More specifically, when D t I N but conditional heteroskedasticity is ignored, the test is, as may be expected, heavily oversized. The obvious conclusion is that the constancy of correlations can only be tested after specifying and estimating both D t and S t .

Appendix C. Details of Maximisation by Parts

This appendix describing the outlines of the estimation algorithm derives from Silvennoinen and Teräsvirta (2021). The estimation proceeds as follows.
  • Assume ln h i t ( θ h i , θ g i ) = 0 , i = 1 , , N , and estimate parameters θ g = ( θ g 1 , , θ g N ) , i = 1 , , N , equation by equation, assuming P t ( θ ρ ) = I N . Denote the estimate S t ( θ ^ g ( 1 , 1 ) ) . This means that the deterministic components g i ( t / T , θ g i ) have been estimated once, including the intercept δ i 0 in (2).
  • Estimate P t ( θ ρ ) given θ g = θ ^ g ( 1 , 1 ) . This requires a separate iteration because P t ( θ ρ ) is nonlinear in parameters; see (5) and (6). Denote the estimate P t ( θ ^ ρ ( 1 , 1 ) ) .
  • Re-estimate S t ( θ g ) assuming P t ( θ ρ ) = P t ( θ ^ ρ ( 1 , 1 ) ) . This yields S t ( θ ^ g ( 1 , 2 ) ) . Then, re-estimate P t ( θ ρ ) given θ g = θ ^ g ( 1 , 2 ) . Iterate until convergence. Let the result after R 1 iterations be S t ( θ g ) = S t ( θ ^ g ( 1 , R 1 ) ) and P t ( θ ρ ) = P t ( θ ^ ρ ( 1 , R 1 ) ) . The resulting estimates are maximum likelihood ones under the assumption D t ( θ h , θ g ) = I N .
  • Estimate θ h from D t ( θ h , θ ^ g ( 1 , R 1 ) ) using P t ( θ ρ ) = P t ( θ ^ ρ ( 1 , R 1 ) ) . This is a standard multivariate conditional correlation GARCH estimation step as in Bollerslev (1990), because S t ( θ ^ g ( 1 , R 1 ) ) is fixed and does not affect the maximum and P t ( θ ^ ρ ( 1 , R 1 ) ) is known. In total, steps 1–4 form the first iteration of the maximisation algorithm. Denote the estimate θ ^ h ( 1 ) .
  • Estimate θ g from S t ( θ g ) keeping D t ( θ ^ h ( 1 ) , θ ^ g ( 1 , R 1 ) ) and P t ( θ ^ ρ ( 1 , R 1 ) ) fixed. This step is analogous to the first part of Step 3. The difference is that D t ( θ ^ h ( 1 ) , θ ^ g ( 1 , R 1 ) ) I N . Denote the estimator S t ( θ ^ g ( 2 , 1 ) ) .
  • Estimate P t ( θ ρ ) given θ g = θ ^ g ( 2 , 1 ) and θ h = θ ^ h ( 1 ) . Denote the estimator P t ( θ ^ ρ ( 2 , 1 ) ) . Iterate until convergence, R 2 iterations. The result: S t ( θ g ) = S t ( θ ^ g ( 2 , R 2 ) ) and P t ( θ ρ ) = P t ( θ ^ ρ ( 2 , R 2 ) ) .
  • Estimate θ h from D t ( θ h , θ ^ g ( 2 , R 2 ) ) using P t ( θ ρ ) = P t ( θ ^ ρ ( 2 , R 2 ) ) ( S t ( θ ^ g ( 2 , R 2 ) ) is fixed). The result: θ h = θ ^ h ( 2 ) . This completes the second full iteration.
  • Repeat steps 5–7 and iterate until convergence.
For identification reasons, δ 0 i , i = 1 , , N , is frozen to δ i 0 = δ ^ i 0 ( 1 , R 1 ) . This frees the intercepts in θ h i . Any positive constant would do for δ i 0 ; however, for numerical reasons, the intercepts are fixed to the values they obtain after the first iteration when θ h has not yet been estimated a single time.
In practice, in estimating the slope parameters in transition functions it may be useful to apply the transformation γ i j = exp { η i j } , in which case γ i j need not be restricted when η i j is bounded away from . The motivation for this transformation is that estimating η i j instead of γ i j is numerically convenient in cases where γ i j is large; see Goodwin et al. (2011) or Silvennoinen and Teräsvirta (2016) for discussion.
Another alternative, proposed by Chan and Theoharakis (2011), is to redefine the slope parameter as γ i j = 1 / η i j 2 and estimate η i j . The authors show that this also alleviates the convergence problems sometimes found when γ i j is large. Ekner and Nejstgaard (2013) aim at the same effect by rescaling γ i j to vary between zero and one.

Appendix D. Estimated Transition Equations

This appendix contains the estimated deterministic components in the TV-GARCH Equations (standard deviation estimates in parentheses). Note that the intercept is fixed after the first iteration; hence, it does not have a standard deviation estimate.
ANZ : g ^ 1 t = 2.28 1.234 ( 0.059 ) ( 1 + exp { 5.715 ( 1.223 ) ( t / T 0.404 ( 0.003 ) ) } ) 1 + 12.316 ( 1.518 ) ( 1 + exp { 5.875 ( 0.392 ) ( t / T 0.571 ( 0.002 ) ) } ) 1 11.704 ( 1.514 ) ( 1 + exp { 4.459 ( 0.166 ) ( t / T 0.623 ( 0.004 ) ) } ) 1 .
CBA : g ^ 2 t = 1.35 0.525 ( 0.054 ) ( 1 + exp { 5.638 ( 2.545 ) ( t / T 0.407 ( 0.007 ) ) } ) 1 + 9.257 ( 1.871 ) ( 1 + exp { 5.117 ( 0.374 ) ( t / T 0.574 ( 0.004 ) ) } ) 1 8.944 ( 1.867 ) ( 1 + exp { 4.504 ( 0.252 ) ( t / T 0.621 ( 0.006 ) ) } ) 1 .
NAB : g ^ 3 t = 1.07 + 3.843 ( 1.273 ) ( 1 + exp { 2.518 ( 0.130 ) ( t / T 0.303 ( 0.034 ) ) } ) 1 3.491 ( 1.114 ) ( 1 + exp { 3.787 ( 0.329 ) ( t / T 0.373 ( 0.008 ) ) } ) 1 + 20.026 ( 5.692 ) ( 1 + exp { 4.926 ( 0.229 ) ( t / T 0.576 ( 0.003 ) ) } ) 1 20.039 ( 5.676 ) ( 1 + exp { 4.183 ( 0.123 ) ( t / T 0.609 ( 0.006 ) ) } ) 1 .
WBC : g ^ 4 t = 2.45 3.120 ( 0.554 ) ( 1 + exp { 2.194 ( 0.124 ) ( t / T 0.534 ( 0.034 ) ) } ) 1 + 25.782 ( 12.524 ) ( 1 + exp { 4.569 ( 0.158 ) ( t / T 0.585 ( 0.006 ) ) } ) 1 23.616 ( 12.682 ) ( 1 + exp { 4.767 ( 0.375 ) ( t / T 0.607 ( 0.007 ) ) } ) 1 .
The locations of the transitions are remarkably similar across transitions. The first transition of the WBC equation is very slow. The effect of the transition extends over the whole estimation period and is the reason for the post-crisis decline in the value of g ^ 4 t ; see Figure 5.

Notes

1
2
The operator vecl ( · ) stacks the subdiagonal elements of its argument matrix.
3
See Explanatory Statement, Banking (prudential standard) Determination 2007, Nos 5, 12 and 15. https://www.legislation.gov.au/Details/F2007L04593/ (accessed on 26 January 2023), https://www.legislation.gov.au/Details/F2007L04600/ (accessed on 26 January 2023) and https://www.legislation.gov.au/Details/F2007L04603/ (accessed on 26 January 2023).

References

  1. Amado, Christina, and Timo Teräsvirta. 2008. Modelling Conditional and Unconditional Heteroskedasticity with Smoothly Time-Varying Structure. SSE/EFI Working Paper Series in Economics and Finance 691. Stockholm: Stockholm School of Economics. [Google Scholar]
  2. Amado, Cristina, and Timo Teräsvirta. 2013. Modelling volatility by variance decomposition. Journal of Econometrics 175: 153–65. [Google Scholar] [CrossRef]
  3. Amado, Cristina, and Timo Teräsvirta. 2014. Conditional correlation models of autoregressive conditional heteroscedasticity with nonstationary GARCH equations. Journal of Business and Economic Statistics 32: 69–87. [Google Scholar] [CrossRef]
  4. Amado, Cristina, and Timo Teräsvirta. 2017. Specification and testing of multiplicative time-varying GARCH models with applications. Econometric Reviews 36: 421–46. [Google Scholar] [CrossRef]
  5. Amado, Cristina, Annastiina Silvennoinen, and Timo Teräsvirta. 2017. Modelling and forecasting WIG20 daily returns. Central European Journal of Economic Modelling and Econometrics 9: 173–200. [Google Scholar]
  6. Berben, Robert-Paul, and W. Jos Jansen. 2005. Comovement in international equity markets: A sectoral view. Journal of International Money and Finance 24: 832–57. [Google Scholar] [CrossRef]
  7. Bollerslev, Tim. 1986. Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics 31: 307–27. [Google Scholar] [CrossRef]
  8. Bollerslev, Tim. 1990. Modelling the coherence in short-run nominal exchange rates: A multivariate generalized ARCH model. Review of Economics and Statistics 72: 498–505. [Google Scholar] [CrossRef]
  9. Box, George E. P., and Gwilym M. Jenkins. 1970. Time Series Analysis, Forecasting Furthermore, Control. San Francisco: Holden-Day. [Google Scholar]
  10. Chan, Felix, and Billy Theoharakis. 2011. Estimating m-regimes STAR–GARCH model using QMLE with parameter transformation. Mathematics and Computers in Simulation 81: 1385–96. [Google Scholar] [CrossRef]
  11. Ekner, Line, and Emil Nejstgaard. 2013. Parameter Identification in the Logistic STAR Model. Discussion Paper 13-07. København: Department of Economics, University of Copenhagen. [Google Scholar]
  12. Engle, Robert F. 2002. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. Journal of Business and Economic Statistics 20: 339–50. [Google Scholar] [CrossRef]
  13. Engle, Robert, and Bryan Kelly. 2012. Dynamic equicorrelation. Journal of Business & Economic Statistics 30: 212–28. [Google Scholar]
  14. Feng, Yuanhua. 2004. Simultaneously modeling conditional heteroskedasticity and scale change. Econometric Theory 20: 563–96. [Google Scholar] [CrossRef]
  15. Feng, Yuanhua. 2006. A Local Dynamic Conditional Correlation Model. MPRA Paper 1592. Edinburgh: Maxwell Institute for Mathematical Sciences, Heriot-Watt University. [Google Scholar]
  16. Glosten, Lawrence R., Ravi Jagannathan, and David E. Runkle. 1993. On the relation between the expected value and the volatility of the nominal excess return on stocks. Journal of Finance 48: 1779–801. [Google Scholar] [CrossRef]
  17. Godfrey, Leslie G. 1988. Misspecification Tests in Econometrics. Cambridge: Cambridge University Press. [Google Scholar]
  18. Goodwin, Barry K., Matthew T. Holt, and Jeffrey P. Prestemon. 2011. North American oriented strand board markets, arbitrage activity, and market price dynamics: A smooth transition approach. American Journal of Agricultural Economics 93: 993–1014. [Google Scholar] [CrossRef]
  19. He, Changli, and Timo Teräsvirta. 1999. Properties of moments of a family of GARCH processes. Journal of Econometrics 92: 173–92. [Google Scholar] [CrossRef]
  20. Kang, Jian, Johan Stax Jakobsen, Annastiina Silvennoinen, Timo Teräsvirta, and Glen Wade. 2022. A parsimonious test of constancy of a positive definite correlation matrix in a multivariate time-varying GARCH model. Econometrics 10: 30. [Google Scholar] [CrossRef]
  21. Lütkepohl, Helmut. 1996. Handbook of Matrices. Chichester: John Wiley & Sons. [Google Scholar]
  22. Luukkonen, Ritva, Pentti Saikkonen, and Timo Teräsvirta. 1988. Testing linearity against smooth transition autoregressive models. Biometrika 75: 491–99. [Google Scholar] [CrossRef]
  23. Silvennoinen, Annastiina, and Timo Teräsvirta. 2005. Multivariate Autoregressive Conditional Heteroskedasticity with Smooth Transitions in Conditional Correlations. SSE/EFI Working Paper Series in Economics and Finance No. 577. Stockholm: Stockholm School of Economics. [Google Scholar]
  24. Silvennoinen, Annastiina, and Timo Teräsvirta. 2009. Modelling multivariate autoregressive conditional heteroskedasticity with the double smooth transition conditional correlation GARCH model. Journal of Financial Econometrics 7: 373–411. [Google Scholar] [CrossRef]
  25. Silvennoinen, Annastiina, and Timo Teräsvirta. 2015. Modeling conditional correlations of asset returns: A smooth transition approach. Econometric Reviews 34: 174–97. [Google Scholar] [CrossRef]
  26. Silvennoinen, Annastiina, and Timo Teräsvirta. 2016. Testing constancy of unconditional variance in volatility models by misspecification and specification tests. Studies in Nonlinear Dynamics and Econometrics 20: 347–64. [Google Scholar] [CrossRef]
  27. Silvennoinen, Annastiina, and Timo Teräsvirta. 2021. Consistency and asymptotic normality of maximum likelihood estimators of the multiplicative time-varying smooth transition correlation GARCH model. Econometrics and Statistics. in press. [Google Scholar] [CrossRef]
  28. Song, Peter X.-K., Yanqin Fan, and John D. Kalbfleisch. 2005. Maximization by parts in likelihood inference. Journal of the American Statistical Association 100: 1145–58. [Google Scholar] [CrossRef]
  29. Teräsvirta, Timo. 1994. Specification, estimation, and evaluation of smooth transition autoregressive models. Journal of the American Statistical Association 89: 208–18. [Google Scholar]
  30. Teräsvirta, Timo, Dag Tjøstheim, and Clive W. J. Granger. 2010. Modelling Nonlinear Economic Time Series. Oxford: Oxford University Press. [Google Scholar]
  31. Tse, Yiu Kuen, and Albert K. C. Tsui. 2002. A multivariate generalized autoregressive conditional heteroscedasticity model with time-varying correlations. Journal of Business and Economic Statistics 20: 351–62. [Google Scholar] [CrossRef]
Figure 1. The market capitalisation of the Big Four as percentage of ASX200 (left) and of ASX200 Financials Index (right).
Figure 1. The market capitalisation of the Big Four as percentage of ASX200 (left) and of ASX200 Financials Index (right).
Econometrics 11 00005 g001
Figure 2. Daily returns of the Big Four, from 2 January 1992 to 31 January 2020.
Figure 2. Daily returns of the Big Four, from 2 January 1992 to 31 January 2020.
Econometrics 11 00005 g002
Figure 3. The first 100 autocorrelations of squared returns.
Figure 3. The first 100 autocorrelations of squared returns.
Econometrics 11 00005 g003
Figure 4. The first 100 autocorrelations of squared standardised returns ε i t 2 / g ^ i t .
Figure 4. The first 100 autocorrelations of squared standardised returns ε i t 2 / g ^ i t .
Econometrics 11 00005 g004
Figure 5. Estimated multiplicative component g ^ i t 1 / 2 (solid curve) and the absolute returns | ε i t | (grey area).
Figure 5. Estimated multiplicative component g ^ i t 1 / 2 (solid curve) and the absolute returns | ε i t | (grey area).
Econometrics 11 00005 g005
Figure 6. Estimated conditional variance h ^ i t from the GJR-GARCH model.
Figure 6. Estimated conditional variance h ^ i t from the GJR-GARCH model.
Econometrics 11 00005 g006
Figure 7. Estimated conditional variance h ^ i t from the TV-GJR-GARCH model.
Figure 7. Estimated conditional variance h ^ i t from the TV-GJR-GARCH model.
Econometrics 11 00005 g007
Figure 8. Estimated correlations. Vertical lines correspond to October 2007 and February 2008.
Figure 8. Estimated correlations. Vertical lines correspond to October 2007 and February 2008.
Econometrics 11 00005 g008
Figure 9. Estimated pairwise correlations.
Figure 9. Estimated pairwise correlations.
Econometrics 11 00005 g009
Table 1. Univariate estimation results for the four banks. GJR is the GJR-GARCH(1,1) equation, and TV-GJR is the TV-GJR-GARCH equation; standard errors in parentheses.
Table 1. Univariate estimation results for the four banks. GJR is the GJR-GARCH(1,1) equation, and TV-GJR is the TV-GJR-GARCH equation; standard errors in parentheses.
α i 0 α i 1 κ i 1 β i 1 PersistenceKurtosis
ANZGJR 0.020 ( 0.005 ) 0.039 ( 0.007 ) 0.044 ( 0.008 ) 0.929 ( 0.008 ) 0.9913.76
TV-GJR 0.111 ( 0.016 ) 0.015 ( 0.005 ) 0.046 ( 0.007 ) 0.792 ( 0.027 ) 0.8313.02
CBAGJR 0.035 ( 0.005 ) 0.060 ( 0.008 ) 0.063 ( 0.011 ) 0.886 ( 0.010 ) 0.9773.66
TV-GJR 0.107 ( 0.014 ) 0.021 ( 0.006 ) 0.065 ( 0.010 ) 0.813 ( 0.020 ) 0.8673.06
NABGJR 0.065 ( 0.009 ) 0.077 ( 0.010 ) 0.075 ( 0.014 ) 0.850 ( 0.014 ) 0.9643.68
TV-GJR 0.152 ( 0.019 ) 0.021 ( 0.006 ) 0.058 ( 0.009 ) 0.731 ( 0.030 ) 0.7803.03
WBCGJR 0.031 ( 0.006 ) 0.045 ( 0.007 ) 0.058 ( 0.010 ) 0.910 ( 0.009 ) 0.9853.70
TV-GJR 0.079 ( 0.011 ) 0.015 ( 0.004 ) 0.041 ( 0.006 ) 0.829 ( 0.020 ) 0.8643.02
Table 2. Estimation results for the four banks’ time-varying correlations. A total of 90% of the estimated transition is between the dates 18 October 2006 and 28 February 2008. The centre point of the location corresponds to 28 June 2007 with ± two standard error ranges of 11 May–13 August 2007.
Table 2. Estimation results for the four banks’ time-varying correlations. A total of 90% of the estimated transition is between the dates 18 October 2006 and 28 February 2008. The centre point of the location corresponds to 28 June 2007 with ± two standard error ranges of 11 May–13 August 2007.
P ( 1 ) P ( 2 )
ANZCBANAB ANZCBANAB
CBA 0.485 ( 0.011 ) CBA 0.782 ( 0.006 )
NAB 0.503 ( 0.010 ) 0.525 ( 0.010 ) NAB 0.808 ( 0.005 ) 0.787 ( 0.005 )
WBC 0.606 ( 0.009 ) 0.500 ( 0.011 ) 0.492 ( 0.011 ) WBC 0.830 ( 0.004 ) 0.818 ( 0.005 ) 0.814 ( 0.005 )
Transition parameters:c η
0.552 ( 0.002 ) 5.020 ( 0.162 )
Note: η = ln γ; see Appendix A.2.
Table 3. Estimation results for the four banks’ time-varying bivariate correlations.
Table 3. Estimation results for the four banks’ time-varying bivariate correlations.
P ( 1 ) P ( 2 )
ANZCBANAB ANZCBANAB
CBA 0.484 ( 0.011 ) CBA 0.784 ( 0.006 )
NAB 0.510 ( 0.010 ) 0.518 ( 0.011 ) NAB 0.811 ( 0.005 ) 0.785 ( 0.005 )
WBC 0.607 ( 0.009 ) 0.504 ( 0.011 ) 0.490 ( 0.011 ) WBC 0.831 ( 0.004 ) 0.816 ( 0.005 ) 0.812 ( 0.005 )
Transition parameters:
c η
ANZCBANAB ANZCBANAB
CBA 0.550 ( 0.005 ) CBA 4.764 ( 0.280 )
NAB 0.567 ( 0.001 ) 0.532 ( 0.008 ) NAB 7.000 ( ) 4.514 ( 0.341 )
WBC 0.555 ( 0.005 ) 0.547 ( 0.004 ) 0.549 ( 0.004 ) WBC 5.198 ( 0.308 ) 4.761 ( 0.254 ) 4.873 ( 0.294 )
Note: η = ln γ; see Appendix A.2.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hall, A.D.; Silvennoinen, A.; Teräsvirta, T. Building Multivariate Time-Varying Smooth Transition Correlation GARCH Models, with an Application to the Four Largest Australian Banks. Econometrics 2023, 11, 5. https://doi.org/10.3390/econometrics11010005

AMA Style

Hall AD, Silvennoinen A, Teräsvirta T. Building Multivariate Time-Varying Smooth Transition Correlation GARCH Models, with an Application to the Four Largest Australian Banks. Econometrics. 2023; 11(1):5. https://doi.org/10.3390/econometrics11010005

Chicago/Turabian Style

Hall, Anthony D., Annastiina Silvennoinen, and Timo Teräsvirta. 2023. "Building Multivariate Time-Varying Smooth Transition Correlation GARCH Models, with an Application to the Four Largest Australian Banks" Econometrics 11, no. 1: 5. https://doi.org/10.3390/econometrics11010005

APA Style

Hall, A. D., Silvennoinen, A., & Teräsvirta, T. (2023). Building Multivariate Time-Varying Smooth Transition Correlation GARCH Models, with an Application to the Four Largest Australian Banks. Econometrics, 11(1), 5. https://doi.org/10.3390/econometrics11010005

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop