Forecasting volatility and co-volatility of crude oil and gold futures: Effects of leverage,jumps, spillovers,and geopolitical risks

To forecast the covariance matrix for the returns of crude oil and gold futures, this paper examines the effects of leverage, jumps, spillovers, and geopolitical risks by using their respective realized covariance matrices. To guarantee the positive definiteness of the forecasts, we consider the full BEKK structure on the conditional Wishart model. By the specification, we can flexibly divide the direct and spillover effects of volatility feedback, negative returns, and jumps. The empirical analysis indicates the benefits of accommodating the spillover effects of negative returns, and the geopolitical risks indicator for modeling and forecasting the covariance matrix.     

The conditional autoregressive Wishart model for multivariate stock market volatility

We propose a Conditional Autoregressive Wishart (CAW) model for the analysis of realized covariance matrices of asset returns. Our model assumes an autoregressive moving average structure for the scale matrix of the Wishart distribution. It accounts for positive definiteness of covariance matrices without imposing parametric restrictions, and can be estimated by Maximum Likelihood. We also propose extensions of the CAW model obtained by including a Mixed Data Sampling (MIDAS) component and Heterogeneous Autoregressive (HAR) dynamics for long-run fluctuations. The CAW models are applied to realized variances and covariances for five New York Stock Exchange stocks.     

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.     

Multivariate rotated ARCH models

This paper introduces a new class of multivariate volatility models which is easy to estimate using covariance targeting, even with rich dynamics. We call them rotated ARCH (RARCH) models. The basic structure is to rotate the returns and then to fit them using a BEKK-type parameterization of the time-varying covariance whose long-run covariance is the identity matrix. This yields the rotated BEKK (RBEKK) model. The extension to DCC-type parameterizations is given, introducing the rotated DCC (RDCC) model. Inference for these models is computationally attractive, and the asymptotics are standard. The techniques are illustrated using data on the DJIA stocks.     

A matrix exponential spatial specification

We introduce the matrix exponential as a way of modelling spatially dependent data. The matrix exponential spatial specification (MESS) simplifies the log-likelihood allowing a closed form solution to the problem of maximum-likelihood estimation, and greatly simplifies the Bayesian estimation of the model. The MESS can produce estimates and inferences similar to those from conventional spatial autoregressive models, but has analytical, computational, and interpretive advantages. We present maximum likelihood and Bayesian approaches to the estimation of this spatial model specification along with methods of model comparisons over different explanatory variables and spatial specifications.     

Modelling security market events in continuous time: Intensity based,multivariate point process models

A continuous time econometric modelling framework for multivariate financial market event (or ‘transactions’) data is developed in which the model is specified via the vector conditional intensity. Generalised Hawkes models are introduced that incorporate inhibitory events and dependence between trading days. Novel omnibus specification tests based on a multivariate random time change theorem are proposed. A bivariate point process model of the timing of trades and mid-quote changes is then presented for a New York Stock Exchange stock and related to the market microstructure literature. The two-way interaction of trades and quote changes in continuous time is found to be important empirically.     

Some further remarks on a superiority problem in misspecified restricted singular linear models

Razzaghi (1987) examined the difference of covariance matrices of competing estimators in misspecified restricted linear models. Further, Gross, Trenkler and Liski (1998) extended Razzaghi's result by asserting his sufficient condition for nonnegative definiteness of the covariance matrix difference to be also necessary. In this paper, when the covariance matrix of the disturbance vector is nonnegative definite, the necessary and sufficient conditions for nonnegative definiteness of the covariance matrix difference are derived. So above results are strengthened. Received: February 1999     

Generalized spectral testing for multivariate continuous-time models

We develop an omnibus specification test for multivariate continuous-time models using the conditional characteristic function, which often has a convenient closed-form or can be accurately approximated for many multivariate continuous-time models in finance and economics. The proposed test fully exploits the information in the joint conditional distribution of underlying economic processes and hence is expected to have good power in a multivariate context. A class of easy-to-interpret diagnostic procedures is supplemented to gauge possible sources of model misspecification. Our tests are also applicable to discrete-time distribution models. Simulation studies show that the tests provide reliable inference in finite samples.     

Generalizing weak instrument robust IV statistics towards multiple parameters,unrestricted covariance matrices and identification statistics

We generalize the weak instrument robust score or Lagrange multiplier and likelihood ratio instrumental variables (IV) statistics towards multiple parameters and a general covariance matrix so they can be used in the generalized method of moments (GMM). The GMM extension of Moreira's [2003. A conditional likelihood ratio test for structural models. Econometrica 71, 1027–1048] conditional likelihood ratio statistic towards GMM preserves its expression except that it becomes conditional on a statistic that tests the rank of a matrix. We analyze the spurious power decline of Kleibergen's [2002. Pivotal statistics for testing structural parameters in instrumental variables regression. Econometrica 70, 1781–1803, 2005. Testing parameters in GMM without assuming that they are identified. Econometrica 73, 1103–1124] score statistic and show that an independent misspecification pre-test overcomes it. We construct identification statistics that reflect if the confidence sets of the parameters are bounded. A power study and the possible shapes of confidence sets illustrate the analysis.     

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.     

A multivariate GARCH in mean approach to testing uncovered interest parity: evidence from Asia-Pacific foreign exchange markets

The existence of time-varying risk premia in deviations from uncovered interest parity (UIP) is investigated based on a conditional capital asset pricing model (CAPM) using data from four Asia-Pacific foreign exchange markets. A parsimonious multivariate generalized autoregressive conditional heteroskedasticity in mean (GARCH-M) parameterization is employed to model the conditional covariance matrix of excess returns. The empirical results indicate that when each currency is estimated separately with an univariate GARCH-M parameterization, no evidence of time-varying risk premia is found except Malaysian ringgit. However, when all currencies are estimated simultaneously with the multivariate GARCH-M parameterization, strong evidence of time-varying risk premia is detected. As a result, the evidence supports the idea that deviations from UIP are due to a risk premium and not to irrationality among market participants. In addition, the empirical evidence found in this study points out that simply modeling the conditional second moments is not sufficient enough to explain the dynamics of the risk premia. A time-varying price of risk is still needed in addition to the conditional volatility. Finally, significant asymmetric world market volatility shocks are found in Asia-Pacific foreign exchange markets.     

Forecasting multivariate realized stock market volatility

We present a new matrix-logarithm model of the realized covariance matrix of stock returns. The model uses latent factors which are functions of lagged volatility, lagged returns and other forecasting variables. The model has several advantages: it is parsimonious; it does not require imposing parameter restrictions; and, it results in a positive-definite estimated covariance matrix. We apply the model to the covariance matrix of size-sorted stock returns and find that two factors are sufficient to capture most of the dynamics.     

Pre-averaging estimators of the ex-post covariance matrix in noisy diffusion models with non-synchronous data

We show how pre-averaging can be applied to the problem of measuring the ex-post covariance of financial asset returns under microstructure noise and non-synchronous trading. A pre-averaged realised covariance is proposed, and we present an asymptotic theory for this new estimator, which can be configured to possess an optimal convergence rate or to ensure positive semi-definite covariance matrix estimates. We also derive a noise-robust Hayashi–Yoshida estimator that can be implemented on the original data without prior alignment of prices. We uncover the finite sample properties of our estimators with simulations and illustrate their practical use on high-frequency equity data.     

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.     

PRACTITIONERS CORNER A POSITIVE SEMI-DEFINITE COVARIANCE MATRIX FOR HAUSMAN SPECIFICATION TESTS OF CONDITIONAL AND MARGINAL DENSITIES

ABSTRACT When the joint density of data can be factorized into a conditional and marginal densities, Hausman test can be used for diagnosing misspecifications of these densities. However, since common covariance estimates of the difference of the two estimators used in Hausman test need not be positive semi-definite in finite samples, the test statistic may be negative. This paper presents a simple and consistent covariance matrix which is positive semi-definite in any finite sample.     

A martingale approach for testing diffusion models based on infinitesimal operator

I develop an omnibus specification test for diffusion models based on the infinitesimal operator. The infinitesimal operator based identification of the diffusion process is equivalent to a “martingale hypothesis” for the processes obtained by a transformation of the original diffusion model. My test procedure is then constructed by checking the “martingale hypothesis” via a multivariate generalized spectral derivative based approach that delivers a N(0,1) asymptotical null distribution for the test statistic. The infinitesimal operator of the diffusion process is a closed-form function of drift and diffusion terms. Consequently, my test procedure covers both univariate and multivariate diffusion models in a unified framework and is particularly convenient for the multivariate case. Moreover, different transformed martingale processes contain separate information about the drift and diffusion specifications. This motivates me to propose a separate inferential test procedure to explore the sources of rejection when a parametric form is rejected. Simulation studies show that the proposed tests have reasonable size and excellent power performance. An empirical application of my test procedure using Eurodollar interest rates finds that most popular short-rate models are rejected and the drift misspecification plays an important role in such rejections.     

Bootstrapping realized multivariate volatility measures

We propose a bootstrap method for statistics that are a function of multivariate high frequency returns such as realized regression, covariance and correlation coefficients. We show that the finite sample performance of the bootstrap is superior to the existing first-order asymptotic theory. Nevertheless, and contrary to the existing results in the bootstrap literature for regression models subject to error heteroskedasticity, the Edgeworth expansion for the pairs bootstrap that we develop here shows that this method is not second-order accurate. We argue that this is due to the fact that the conditional mean parameters of realized regression models are heterogeneous under stochastic volatility.     

Missing in Asynchronicity: A Kalman‐em Approach for Multivariate Realized Covariance Estimation

Motivated by the need for a positive‐semidefinite estimator of multivariate realized covariance matrices, we model noisy and asynchronous ultra‐high‐frequency asset prices in a state‐space framework with missing data. We then estimate the covariance matrix of the latent states through a Kalman smoother and expectation maximization (KEM) algorithm. Iterating between the two EM steps, we obtain a covariance matrix estimate which is robust to both asynchronicity and microstructure noise, and positive‐semidefinite by construction. We show the performance of the KEM estimator using extensive Monte Carlo simulations that mimic the liquidity and market microstructure characteristics of the S&P 500 universe as well as in a high‐dimensional application on US stocks. KEM provides very accurate covariance matrix estimates and significantly outperforms alternative approaches recently introduced in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

Sequential conditional correlations: Inference and evaluation

This paper presents a new approach to the modeling of conditional correlation matrices within the multivariate GARCH framework. The procedure, which consists of breaking the matrix into the product of a sequence of matrices with desirable characteristics, in effect converts a highly dimensional and intractable optimization problem into a series of simple and feasible estimations. This in turn allows for richer parameterizations and complex functional forms for the single components. An empirical application involving the conditional second moments of 69 selected stocks from the NASDAQ100 shows how the new procedure results in strikingly accurate measures of the conditional correlations.     

Asymptotic distribution of the cointegrating vector estimator in error correction models with conditional heteroskedasticity

This paper explores the asymptotic distribution of the cointegrating vector estimator in error correction models with conditionally heteroskedastic errors. Asymptotic properties of the maximum likelihood estimator (MLE) of the cointegrating vector, which estimates the cointegrating vector and the multivariate GARCH process jointly, are provided. The MLE of the cointegrating vector follows mixture normal, and its asymptotic distribution depends on the conditional heteroskedasticity and the kurtosis of standardized innovations. The reduced rank regression (RRR) estimator and the regression-based cointegrating vector estimators do not consider conditional heteroskedasticity, and thus the efficiency gain of the MLE emerges as the magnitude of conditional heteroskedasticity increases. The simulation results indicate that the relative power of the t-statistics based on the MLE improves significantly as the GARCH effect increases.