Multivariate GARCH DCC Estimation

In this video tutorial, we will demonstrate the process of Multivariate GARCH DCC Estimation using OxMetrics 6. Multivariate Garch can be elaborate in the following written tutorial from: https://vlab.stern.nyu.edu/doc/13?topic=mdls

Definition of Multivariate GARCH

Consider n time series of returns and make the usual assumption that returns are serially uncorrelated. Then, we can define a vector of zero-mean white noises εt=rtμ, where rt is the n1 vector of returns and μ is the vector of expected returns.

Despite of being serially uncorrelated, the returns may present contemporaneous correlation. That is:

t?t1[(rtμ)(rtμ)]

may not be a diagonal matrix. Moreover, this contemporaneous variance may be time-varying, depending on past information.

The GARCH-DCC involves two steps. The first step accounts for the conditional heteroskedasticity. It consists in estimating, for each one of the n series of returns rit, its conditional volatility σit using a GARCH model (see garch documentation). Let Dt be a diagonal matrix with these conditional volatilities, i.e. Di,it=σit and, if ij, Di,jt=0. Then the standardized residuals are:

νtD-1t(rtμ)

and notice that these standardized residuals have unit conditional volatility. Now, define the matrix:

R⎯⎯⎯1Tt=1Tνtνt

This is the Bollerslev's Constant Conditional Correlation (CCC) Estimator (Bollerslev, 1990).

The second step consists in generalizing Bollerslev's CCC to capture dynamics in the correlation, hence the name Dynamic Conditional Correlation (DCC). The DCC correlations are:

Qt=R⎯⎯⎯+α(νt1νt1R⎯⎯⎯)+β(Qt1R⎯⎯⎯)

So, Qi,jt is the correlation between rit and rjt at time t, and that is what is plotted by V-Lab.

Estimation

The estimation of one GARCH model for each of the n time series of returns in the first step is standard. For details on GARCH estimation, see garch documentation.

For the second step, which is the DCC estimation per se, V-Lab estimates both parameters, α and β, simultaneously, by maximizing the log likelihood. The standardized residuals are assumed to be jointly Gaussian. To ease the computation cost of estimating a vast dimensional time-varying correlation model, V-Lab uses a technique called composite likelihood (Engle et al., 2007).

The DCC model captures a stylized facts in financial time series: correlation clustering. The correlation is more likely to be high at time t if it was also high at time t1. Another way of seeing this is noting that a shock at time t1 also impacts the correlation at time t. However, if α+β<1, the correlation itself is mean reverting, and it fluctuates around R⎯⎯⎯, the unconditional correlation.

Usual restrictions on the parameters are α,β>0. Though, it is possible to have α+β=1; the conditional correlation is then an integrated process.

Variance Targeting

Notice that if we had written the DCC model in a fashion similar to the GARCH model:

Qt=Ω+ανt1νt1+βQt1

we would have to estimate the matrix Ω also. That is, instead of estimating only two parameters, we would have to estimate 2+nn+12 parameters (it is not 2+n2 parameters due to the fact that Ω is a symmetric matrix). And then the unconditional correlation implied by the model would have been:

R⎯⎯⎯=Ω1αβ

Instead of estimating Ω, notice that we actually substituted Ω by R⎯⎯⎯(1αβ) in the DCC formula, which is a much more parsimonious way of writing the model. This is called Variance Targeting, introduced by Engle and Mezrich in 1995, and it is a very useful technique when modeling vast dimensional time-varying covariance or correlation models.

DCC(p,q)

The specific model just described can be generalized in two ways.

In the first stage, each GARCH specification used to standardize each one of the n return time series can be generalized to a GARCH(p,q) model (see garch documentation), where p and q can be chosen differently for each return time series, for instance, by Bayesian Information Criterion (BIN), also known as Schwarz Information Criterion (SIC), or by Akaike Information Criterion (AIC). The former tends to be more parsimonious than the latter. V-Lab uses p=1 and q=1 though, because this is usually the option that best fits financial time series.

In the second stage, the DCC model can be generalized to account for more lags in the conditional correlation. A DCC(p,q) model assumes that:

Qt=R⎯⎯⎯+i=1pαi(νtiνtiR⎯⎯⎯)+j=1qβj(QtjR⎯⎯⎯)

where p and q can be chosen, for instance, by information criterion. Again, V-Lab uses p=1 and q=1 though, because this is usually the option that best fits financial time series.

Bibliography

Bollerslev, T., 1990. Modeling The Coherence in Short-Run Nominal Exchange Rates: A Multivariate Generalized ARCH Model. Review of Economics and Statistics 72: 498-505.

Engle, R. F., 2002. Dynamic Conditional Correlation: A Simple Class of Multivariate GARCH Models. Journal of Business and Economic Statistics 20(3).

Engle, R. F., 2009. Anticipating Correlations: A New Paradigm for Risk Management. Princeton University Press.

Engle, R. F. and J. Mezrich, 1995. Grappling with GARCH. Risk: 112-117.

Engle, R. F., N. Shephard, and K. Sheppard, 2007. Fitting and Testing Vast Dimensional Time-Varying Covariance Models. NYU Working Paper FIN-07-046.

Video On Multivariate GARCH DCC Estimation

Leave a Reply

You must be logged in to post a comment.