MGARCH models: Trade-off between feasibility and flexibility

Multivariate GARCH (MGARCH) models need to be restricted so that their estimation is feasible in large systems and so that the covariance stationarity and positive definiteness of conditional covariance matrices are guaranteed. This paper analyzes the limitations of some of the popular restricted parametric MGARCH models that are often used to represent the dynamics observed in real systems of financial returns. These limitations are illustrated using simulated data generated by general VECH models of different dimensions in which volatilities and correlations are interrelated. We show that the restrictions imposed by the BEKK model are very unrealistic, generating potentially misleading forecasts of conditional correlations. On the other hand, models based on the DCC specification provide appropriate forecasts. Alternative estimators of the parameters are important in order to simplify the computations, and do not have implications for the estimates of conditional correlations. The implications of the restrictions imposed by the different specifications of MGARCH models considered are illustrated by forecasting the volatilities and correlations of a five-dimensional system of exchange rate returns.     

A conditionally heteroskedastic independent factor model with an application to financial stock returns

We propose a new conditionally heteroskedastic factor model, the GICA-GARCH model, which combines independent component analysis (ICA) and multivariate GARCH (MGARCH) models. This model assumes that the data are generated by a set of underlying independent components (ICs) that capture the co-movements among the observations, which are assumed to be conditionally heteroskedastic. The GICA-GARCH model separates the estimation of the ICs from their fitting with a univariate ARMA-GARCH model. Here, we will use two ICA approaches to find the ICs: the first estimates the components, maximizing their non-Gaussianity, while the second exploits the temporal structure of the data. After estimating and identifying the common ICs, we fit a univariate GARCH model to each of them in order to estimate their univariate conditional variances. The GICA-GARCH model then provides a new framework for modelling the multivariate conditional heteroskedasticity in which we can explain and forecast the conditional covariances of the observations by modelling the univariate conditional variances of a few common ICs. We report some simulation experiments to show the ability of ICA to discover leading factors in a multivariate vector of financial data. Finally, we present an empirical application to the Madrid stock market, where we evaluate the forecasting performances of the GICA-GARCH and two additional factor GARCH models: the orthogonal GARCH and the conditionally uncorrelated components GARCH.     

REALIZED BETA GARCH: A MULTIVARIATE GARCH MODEL WITH REALIZED MEASURES OF VOLATILITY

We introduce a multivariate generalized autoregressive conditional heteroskedasticity (GARCH) model that incorporates realized measures of variances and covariances. Realized measures extract information about the current levels of volatilities and correlations from high‐frequency data, which is particularly useful for modeling financial returns during periods of rapid changes in the underlying covariance structure. When applied to market returns in conjunction with returns on an individual asset, the model yields a dynamic model specification of the conditional regression coefficient that is known as the beta. We apply the model to a large set of assets and find the conditional betas to be far more variable than usually found with rolling‐window regressions based exclusively on daily returns. In the empirical part of the paper, we examine the cross‐sectional as well as the time variation of the conditional beta series during the financial crises. Copyright © 2014 John Wiley & Sons, Ltd.     

Testing multivariate distributions in GARCH models

In this paper, we consider testing distributional assumptions in multivariate GARCH models based on empirical processes. Using the fact that joint distribution carries the same amount of information as the marginal together with conditional distributions, we first transform the multivariate data into univariate independent data based on the marginal and conditional cumulative distribution functions. We then apply the Khmaladze's martingale transformation (K-transformation) to the empirical process in the presence of estimated parameters. The K-transformation eliminates the effect of parameter estimation, allowing a distribution-free test statistic to be constructed. We show that the K-transformation takes a very simple form for testing multivariate normal and multivariate t-distributions. The procedure is applied to a multivariate financial time series data set.     

Copula-based multivariate GARCH model with uncorrelated dependent errors

Multivariate GARCH (MGARCH) models are usually estimated under multivariate normality. In this paper, for non-elliptically distributed financial returns, we propose copula-based multivariate GARCH (C-MGARCH) model with uncorrelated dependent errors, which are generated through a linear combination of dependent random variables. The dependence structure is controlled by a copula function. Our new C-MGARCH model nests a conventional MGARCH model as a special case. The aim of this paper is to model MGARCH for non-normal multivariate distributions using copulas. We model the conditional correlation (by MGARCH) and the remaining dependence (by a copula) separately and simultaneously. We apply this idea to three MGARCH models, namely, the dynamic conditional correlation (DCC) model of Engle [Engle, R.F., 2002. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. Journal of Business and Economic Statistics 20, 339–350], the varying correlation (VC) model of Tse and Tsui [Tse, Y.K., Tsui, A.K., 2002. A multivariate generalized autoregressive conditional heteroscedasticity model with time-varying correlations. Journal of Business and Economic Statistics 20, 351–362], and the BEKK model of Engle and Kroner [Engle, R.F., Kroner, K.F., 1995. Multivariate simultaneous generalized ARCH. Econometric Theory 11, 122–150]. Empirical analysis with three foreign exchange rates indicates that the C-MGARCH models outperform DCC, VC, and BEKK in terms of in-sample model selection and out-of-sample multivariate density forecast, and in terms of these criteria the choice of copula functions is more important than the choice of the volatility models.     

DCC- and DECO-HEAVY: Multivariate GARCH models based on realized variances and correlations

This paper introduces the scalar DCC-HEAVY and DECO-HEAVY models for conditional variances and correlations of daily returns based on measures of realized variances and correlations built from intraday data. Formulas for multi-step forecasts of conditional variances and correlations are provided. Asymmetric versions of the models are developed. An empirical study shows that in terms of forecasts the scalar HEAVY models outperform the scalar BEKK-HEAVY model based on realized covariances and the scalar BEKK, DCC, and DECO multivariate GARCH models based exclusively on daily data.     

Time series copulas for heteroskedastic data

We propose parametric copulas that capture serial dependence in stationary heteroskedastic time series. We suggest copulas for first‐order Markov series, and then extend them to higher orders and multivariate series. We derive the copula of a volatility proxy, based on which we propose new measures of volatility dependence, including co‐movement and spillover in multivariate series. In general, these depend upon the marginal distributions of the series. Using exchange rate returns, we show that the resulting copula models can capture their marginal distributions more accurately than univariate and multivariate generalized autoregressive conditional heteroskedasticity models, and produce more accurate value‐at‐risk forecasts.     

Forecasting the volatility of asset returns: The informational gains from option prices

A new class of forecasting models is proposed that extends the realized GARCH class of models through the inclusion of option prices to forecast the variance of asset returns. The VIX is used to approximate option prices, resulting in a set of cross-equation restrictions on the model’s parameters. The full model is characterized by a nonlinear system of three equations containing asset returns, the realized variance, and the VIX, with estimation of the parameters based on maximum likelihood methods. The forecasting properties of the new class of forecasting models, as well as a number of special cases, are investigated and applied to forecasting the daily S&P500 index realized variance using intra-day and daily data from September 2001 to November 2017. The forecasting results provide strong support for including the realized variance and the VIX to improve variance forecasts, with linear conditional variance models performing well for short-term one-day-ahead forecasts, whereas log-linear conditional variance models tend to perform better for intermediate five-day-ahead forecasts.     

The structure of dynamic correlations in multivariate stochastic volatility models

This paper proposes two types of stochastic correlation structures for Multivariate Stochastic Volatility (MSV) models, namely the constant correlation (CC) MSV and dynamic correlation (DC) MSV models, from which the stochastic covariance structures can easily be obtained. Both structures can be used for purposes of determining optimal portfolio and risk management strategies through the use of correlation matrices, and for calculating Value-at-Risk (VaR) forecasts and optimal capital charges under the Basel Accord through the use of covariance matrices. A technique is developed to estimate the DC MSV model using the Markov Chain Monte Carlo (MCMC) procedure, and simulated data show that the estimation method works well. Various multivariate conditional volatility and MSV models are compared via simulation, including an evaluation of alternative VaR estimators. The DC MSV model is also estimated using three sets of empirical data, namely Nikkei 225 Index, Hang Seng Index and Straits Times Index returns, and significant dynamic correlations are found. The Dynamic Conditional Correlation (DCC) model is also estimated, and is found to be far less sensitive to the covariation in the shocks to the indexes. The correlation process for the DCC model also appears to have a unit root, and hence constant conditional correlations in the long run. In contrast, the estimates arising from the DC MSV model indicate that the dynamic correlation process is stationary.     

Extreme returns and intensity of trading

Asymmetric information models of market microstructure claim that variables such as trading intensity are proxies for latent information on the value of financial assets. We consider the interval‐valued time series (ITS) of low/high returns and explore the relationship between these extreme returns and the intensity of trading. We assume that the returns (or prices) are generated by a latent process with some unknown conditional density. At each period of time, from this density, we have some random draws (trades) and the lowest and highest returns are the realized extreme observations of the latent process over the sample of draws. In this context, we propose a semiparametric model of extreme returns that exploits the results provided by extreme value theory. If properly centered and standardized extremes have well‐defined limiting distributions, the conditional mean of extreme returns is a nonlinear function of the conditional moments of the latent process and of the conditional intensity of the process that governs the number of draws. We implement a two‐step estimation procedure. First, we estimate parametrically the regressors that will enter into the nonlinear function, and in a second step we estimate nonparametrically the conditional mean of extreme returns as a function of the generated regressors. Unlike current models for ITS, the proposed semiparametric model is robust to misspecification of the conditional density of the latent process. We fit several nonlinear and linear models to the 5‐minute and 1‐minute low/high returns to seven major banks and technology stocks, and find that the nonlinear specification is superior to the current linear models and that the conditional volatility of the latent process and the conditional intensity of the trading process are major drivers of the dynamics of extreme returns.     

Order-invariant tests for proper calibration of multivariate density forecasts

Established tests for proper calibration of multivariate density forecasts based on Rosenblatt probability integral transforms can be manipulated by changing the order of variables in the forecasting model. We derive order-invariant tests. The new tests are applicable to densities of arbitrary dimensions and can deal with parameter estimation uncertainty and dynamic misspecification. Monte Carlo simulations show that they often have superior power relative to established approaches. We use the tests to evaluate generalized autoregressive conditional heteroskedasticity-based multivariate density forecasts for a vector of stock market returns and macroeconomic forecasts from a Bayesian vector autoregression with time-varying parameters.     

Estimating and predicting multivariate volatility thresholds in global stock markets

We propose a general double tree structured AR‐GARCH model for the analysis of global equity index returns. The model extends previous approaches by incorporating (i) several multivariate thresholds in conditional means and volatilities of index returns and (ii) a richer specification for the impact of lagged foreign (US) index returns in each threshold. We evaluate the out‐of‐sample forecasting power of our model for eight major equity indices in comparison to some existing volatility models in the literature. We find strong evidence for more than one multivariate threshold (more than two regimes) in conditional means and variances of global equity index returns. Such multivariate thresholds are affected by foreign (US) lagged index returns and yield a higher out‐of‐sample predictive power for our tree structured model setting. Copyright © 2006 John Wiley & Sons, Ltd.     

The value of multivariate model sophistication: An application to pricing Dow Jones Industrial Average options

We assess the predictive accuracies of a large number of multivariate volatility models in terms of pricing options on the Dow Jones Industrial Average. We measure the value of model sophistication in terms of dollar losses by considering a set of 444 multivariate models that differ in their specification of the conditional variance, conditional correlation, innovation distribution, and estimation approach. All of the models belong to the dynamic conditional correlation class, which is particularly suitable because it allows consistent estimations of the risk neutral dynamics with a manageable amount of computational effort for relatively large scale problems. It turns out that increasing the sophistication in the marginal variance processes (i.e., nonlinearity, asymmetry and component structure) leads to important gains in pricing accuracy. Enriching the model with more complex existing correlation specifications does not improve the performance significantly. Estimating the standard dynamic conditional correlation model by composite likelihood, in order to take into account potential biases in the parameter estimates, generates only slightly better results. To enhance this poor performance of correlation models, we propose a new model that allows for correlation spillovers without too many parameters. This model performs about 60% better than the existing correlation models we consider. Relaxing a Gaussian innovation for a Laplace innovation assumption improves the pricing in a more minor way. In addition to investigating the value of model sophistication in terms of dollar losses directly, we also use the model confidence set approach to statistically infer the set of models that delivers the best pricing performances.     

Forecasting cryptocurrencies under model and parameter instability

This paper studies the predictability of cryptocurrency time series. We compare several alternative univariate and multivariate models for point and density forecasting of four of the most capitalized series: Bitcoin, Litecoin, Ripple and Ethereum. We apply a set of crypto-predictors and rely on dynamic model averaging to combine a large set of univariate dynamic linear models and several multivariate vector autoregressive models with different forms of time variation. We find statistically significant improvements in point forecasting when using combinations of univariate models, and in density forecasting when relying on the selection of multivariate models. Both schemes deliver sizable directional predictability.     

Conditional correlations and volatility spillovers between crude oil and stock index returns

This paper investigates the conditional correlations and volatility spillovers between the crude oil and financial markets, based on crude oil returns and stock index returns. Daily returns from 2 January 1998 to 4 November 2009 of the crude oil spot, forward and futures prices from the WTI and Brent markets, and the FTSE100, NYSE, Dow Jones and S&P500 stock index returns, are analysed using the CCC model of Bollerslev (1990), VARMA-GARCH model of Ling and McAleer (2003), VARMA-AGARCH model of McAleer, Hoti, and Chan (2008), and DCC model of Engle (2002). Based on the CCC model, the estimates of conditional correlations for returns across markets are very low, and some are not statistically significant, which means the conditional shocks are correlated only in the same market and not across markets. However, the DCC estimates of the conditional correlations are always significant. This result makes it clear that the assumption of constant conditional correlations is not supported empirically. Surprisingly, the empirical results from the VARMA-GARCH and VARMA-AGARCH models provide little evidence of volatility spillovers between the crude oil and financial markets. The evidence of asymmetric effects of negative and positive shocks of equal magnitude on the conditional variances suggests that VARMA-AGARCH is superior to VARMA-GARCH and CCC.     

Econometric modelling in finance and risk management: An overview

This paper gives an overview about the sixteen papers included in this special issue. The papers in this special issue cover a wide range of topics. Such topics include discussing a class of tests for correlation, estimation of realized volatility, modeling time series and continuous-time models with long-range dependence, estimation and specification testing of time series models, estimation in a factor model with high-dimensional problems, finite-sample examination of quasi-maximum likelihood estimation in an autoregressive conditional duration model, and estimation in a dynamic additive quantile model.     

On a double-threshold autoregressive heteroscedastic time series model

Tong's threshold models have been found useful in modelling nonlinearities in the conditional mean of a time series. The threshold model is extended to the so-called double-threshold ARCH(DTARCH) model, which can handle the situation where both the conditional mean and the conditional variance specifications are piecewise linear given previous information. Potential applications of such models include financial data with different (asymmetric) behaviour in a rising versus a falling market and business cycle modelling. Model identification, estimation and diagnostic checking techniques are developed. Maximum likelihood estimation can be achieved via an easy-to-use iteratively weighted least squares algorithm. Portmanteau-type statistics are also derived for checking model adequacy. An illustrative example demonstrates that asymmetric behaviour in the mean and the variance could be present in financial series and that the DTARCH model is capable of capturing these phenomena.     

An efficient Bayesian approach to multiple structural change in multivariate time series

This paper provides a feasible approach to estimation and forecasting of multiple structural breaks for vector autoregressions and other multivariate models. Owing to conjugate prior assumptions we obtain a very efficient sampler for the regime allocation variable. A new hierarchical prior is introduced to allow for learning over different structural breaks. The model is extended to independent breaks in regression coefficients and the volatility parameters. Two empirical applications show the improvements the model has over benchmarks. In a macro application with seven variables we empirically demonstrate the benefits from moving from a multivariate structural break model to a set of univariate structural break models to account for heterogeneous break patterns across data series.     

基于极值分布理论的VaR与ES度量

本文应用极值分布理论对金融收益序列的尾部进行估计,计算收益序列的在险价值VaR和预期不足ES来度量市场风险。通过伪最大似然估计方法估计的GARCH模型对收益数据进行拟合,应用极值理论中的GPD对新息分布的尾部建模,得到了基于尾部估计产生收益序列的VaR和ES值。采用上证指数日对数收益数据为样本,得到了度量条件极值和无条件极值下VaR和ES的结果。实证研究表明:在置信水平很高(如99%)的条件下,采用极值方法度量风险值效果更好。而置信水平在95%下,其他方法和极值方法结合效果会很好。用ES度量风险能够使我们了解不利情况发生时风险的可能情况。