Analysis of high dimensional multivariate stochastic volatility models

This paper is concerned with the Bayesian estimation and comparison of flexible, high dimensional multivariate time series models with time varying correlations. The model proposed and considered here combines features of the classical factor model with that of the heavy tailed univariate stochastic volatility model. A unified analysis of the model, and its special cases, is developed that encompasses estimation, filtering and model choice. The centerpieces of the estimation algorithm (which relies on MCMC methods) are: (1) a reduced blocking scheme for sampling the free elements of the loading matrix and the factors and (2) a special method for sampling the parameters of the univariate SV process. The resulting algorithm is scalable in terms of series and factors and simulation-efficient. Methods for estimating the log-likelihood function and the filtered values of the time-varying volatilities and correlations are also provided. The performance and effectiveness of the inferential methods are extensively tested using simulated data where models up to 50 dimensions and 688 parameters are fit and studied. The performance of our model, in relation to various multivariate GARCH models, is also evaluated using a real data set of weekly returns on a set of 10 international stock indices. We consider the performance along two dimensions: the ability to correctly estimate the conditional covariance matrix of future returns and the unconditional and conditional coverage of the 5% and 1% value-at-risk (VaR) measures of four pre-defined portfolios.     

Bayesian semiparametric stochastic volatility modeling

This paper extends the existing fully parametric Bayesian literature on stochastic volatility to allow for more general return distributions. Instead of specifying a particular distribution for the return innovation, nonparametric Bayesian methods are used to flexibly model the skewness and kurtosis of the distribution while the dynamics of volatility continue to be modeled with a parametric structure. Our semiparametric Bayesian approach provides a full characterization of parametric and distributional uncertainty. A Markov chain Monte Carlo sampling approach to estimation is presented with theoretical and computational issues for simulation from the posterior predictive distributions. An empirical example compares the new model to standard parametric stochastic volatility models.     

A threshold stochastic volatility model with explanatory variables

In this paper, we introduce a threshold stochastic volatility model with explanatory variables. The Bayesian method is considered in estimating the parameters of the proposed model via the Markov chain Monte Carlo (MCMC) algorithm. Gibbs sampling and Metropolis–Hastings sampling methods are used for drawing the posterior samples of the parameters and the latent variables. In the simulation study, the accuracy of the MCMC algorithm, the sensitivity of the algorithm for model assumptions, and the robustness of the posterior distribution under different priors are considered. Simulation results indicate that our MCMC algorithm converges fast and that the posterior distribution is robust under different priors and model assumptions. A real data example was analyzed to explain the asymmetric behavior of stock markets.     

Leverage effect on stochastic volatility for option pricing in Hong Kong: A simulation and empirical study

This paper explores the importance of incorporating the financial leverage effect in the stochastic volatility models when pricing options. For the illustrative purpose, we first conduct the simulation experiment by using the Markov Chain Monte Carlo (MCMC) sampling method. We then make an empirical analysis by applying the volatility models to the real return data of the Hang Seng index during the period from January 1, 2013 to December 31, 2017. Our results highlight the accuracy of the stochastic volatility models with leverage in option pricing when leverage is high. In addition, the leverage effect becomes more significant as the maturity of options increases. Moreover, leverage affects the pricing of in-the-money options more than that of at-the-money and out-of-money options. Our study is therefore useful for both asset pricing and portfolio investment in the Hong Kong market where volatility is an inherent nature of the economy.     

Estimation of affine term structure models with spanned or unspanned stochastic volatility

We develop new procedures for maximum likelihood estimation of affine term structure models with spanned or unspanned stochastic volatility. Our approach uses linear regression to reduce the dimension of the numerical optimization problem yet it produces the same estimator as maximizing the likelihood. It improves the numerical behavior of estimation by eliminating parameters from the objective function that cause problems for conventional methods. We find that spanned models capture the cross-section of yields well but not volatility while unspanned models fit volatility at the expense of fitting the cross-section.     

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.     

The Wishart Autoregressive process of multivariate stochastic volatility

The Wishart Autoregressive (WAR) process is a dynamic model for time series of multivariate stochastic volatility. The WAR naturally accommodates the positivity and symmetry of volatility matrices and provides closed-form non-linear forecasts. The estimation of the WAR is straighforward, as it relies on standard methods such as the Method of Moments and Maximum Likelihood. For illustration, the WAR is applied to a sequence of intraday realized volatility–covolatility matrices from the Toronto Stock Market (TSX).     

Dynamic correlation multivariate stochastic volatility with latent factors

Modeling the correlation structure of returns is essential in many financial applications. Considerable evidence from empirical studies has shown that the correlation among asset returns is not stable over time. A recent development in the multivariate stochastic volatility literature is the application of inverse Wishart processes to characterize the evolution of return correlation matrices. Within the inverse Wishart multivariate stochastic volatility framework, we propose a flexible correlated latent factor model to achieve dimension reduction and capture the stylized fact of ‘correlation breakdown’ simultaneously. The parameter estimation is based on existing Markov chain Monte Carlo methods. We illustrate the proposed model with several empirical studies. In particular, we use high‐dimensional stock return data to compare our model with competing models based on multiple performance metrics and tests. The results show that the proposed model not only describes historic stylized facts reasonably but also provides the best overall performance.     

Stochastic volatility with leverage: Fast and efficient likelihood inference

This paper is concerned with the Bayesian analysis of stochastic volatility (SV) models with leverage. Specifically, the paper shows how the often used Kim et al. [1998. Stochastic volatility: likelihood inference and comparison with ARCH models. Review of Economic Studies 65, 361–393] method that was developed for SV models without leverage can be extended to models with leverage. The approach relies on the novel idea of approximating the joint distribution of the outcome and volatility innovations by a suitably constructed ten-component mixture of bivariate normal distributions. The resulting posterior distribution is summarized by MCMC methods and the small approximation error in working with the mixture approximation is corrected by a reweighting procedure. The overall procedure is fast and highly efficient. We illustrate the ideas on daily returns of the Tokyo Stock Price Index. Finally, extensions of the method are described for superposition models (where the log-volatility is made up of a linear combination of heterogenous and independent autoregressions) and heavy-tailed error distributions (student and log-normal).     

Estimation methods for stochastic volatility models: a survey

Abstract.  Although stochastic volatility (SV) models have an intuitive appeal, their empirical application has been limited mainly due to difficulties involved in their estimation. The main problem is that the likelihood function is hard to evaluate. However, recently, several new estimation methods have been introduced and the literature on SV models has grown substantially. In this article, we review this literature. We describe the main estimators of the parameters and the underlying volatilities focusing on their advantages and limitations both from the theoretical and empirical point of view. We complete the survey with an application of the most important procedures to the S&P 500 stock price index.     

Bayesian estimation of an extended local scale stochastic volatility model

A new version of the local scale model of Shephard (1994) is presented. Its features are identically distributed evolution equation disturbances, the incorporation of in-the-mean effects, and the incorporation of variance regressors. A Bayesian posterior simulator and a new simulation smoother are presented. The model is applied to publicly available daily exchange rate and asset return series, and is compared with t-GARCH and Lognormal stochastic volatility formulations using Bayes factors.     

Analytical quasi maximum likelihood inference in multivariate volatility models

Quasi maximum likelihood estimation and inference in multivariate volatility models remains a challenging computational task if, for example, the dimension of the parameter space is high. One of the reasons is that typically numerical procedures are used to compute the score and the Hessian, and often they are numerically unstable. We provide analytical formulae for the score and the Hessian for a variety of multivariate GARCH models including the Vec and BEKK specifications as well as the recent dynamic conditional correlation model. By means of a Monte Carlo investigation of the BEKK–GARCH model we illustrate that employing analytical derivatives for inference is clearly preferable to numerical methods.     

Spatial stochastic frontier models: accounting for unobserved local determinants of inefficiency

This paper analyzes the productivity of farms across 370 municipalities in the Center-West region of Brazil. A stochastic frontier model with a latent spatial structure is proposed to account for possible unknown geographical variation of the outputs. The paper compares versions of the model that include the latent spatial effect in the mean of output or as a variable that conditions the distribution of inefficiency, include or not observed municipal variables, and specify independent normal or conditional autoregressive priors for the spatial effects. The Bayesian paradigm is used to estimate the proposed models. As the resultant posterior distributions do not have a closed form, stochastic simulation techniques are used to obtain samples from them. Two model comparison criteria provide support for including the latent spatial effects, even after considering covariates at the municipal level. Models that ignore the latent spatial effects produce significantly different rankings of inefficiencies across agents.

Bayesian stochastic frontier analysis using WinBUGS

Markov chain Monte Carlo (MCMC) methods have become a ubiquitous tool in Bayesian analysis. This paper implements MCMC methods for Bayesian analysis of stochastic frontier models using the WinBUGS package, a freely available software. General code for cross-sectional and panel data are presented and various ways of summarizing posterior inference are discussed. Several examples illustrate that analyses with models of genuine practical interest can be performed straightforwardly and model changes are easily implemented. Although WinBUGS may not be that efficient for more complicated models, it does make Bayesian inference with stochastic frontier models easily accessible for applied researchers and its generic structure allows for a lot of flexibility in model specification.      

Cross-sectional implications of dynamic asset pricing with stochastic volatility and ambiguity aversion

Building on recent research that highlights the importance of macroeconomic volatility and ambiguity aversion in explaining the dynamics of stock returns, in this paper we propose a dynamic asset pricing model that simultaneously accounts for stochastic macroeconomic volatility and ambiguity, assuming that investors deal with uncertainty about the mechanics of macroeconomic fluctuations using first-release consumption and revisions to aggregate consumption on vintage data. Our results show that the proposed model captures a large fraction of the cross-sectional variation of excess returns for a wide range of market anomaly portfolios. Furthermore, while the price of risk for ambiguity is positive and significant for the vast majority of assets under study, macroeconomic volatility yields ambiguous outcomes, although it significantly increases the explanatory power of the model for specific assets. Our results suggest that macroeconomic volatility and ambiguity complement each other in explaining the cross-sectional behavior of stock returns.     

Semi-parametric dynamic asymmetric Laplace models for tail risk forecasting,incorporating realized measures

This paper extends the joint Value-at-Risk (VaR) and expected shortfall (ES) quantile regression model of Taylor (2019), by incorporating a realized measure to drive the tail risk dynamics, as a potentially more efficient driver than daily returns. Furthermore, we propose and test a new model for the dynamics of the ES component. Both a maximum likelihood and an adaptive Bayesian Markov chain Monte Carlo method are employed for estimation, the properties of which are compared in a simulation study. The results favour the Bayesian approach, which is employed subsequently in a forecasting study of seven financial market indices. The proposed models are compared to a range of parametric, non-parametric and semi-parametric competitors, including GARCH, realized GARCH, the extreme value theory method and the joint VaR and ES models of Taylor (2019), in terms of the accuracy of one-day-ahead VaR and ES forecasts, over a long forecast sample period that includes the global financial crisis in 2007–2008. The results are favorable for the proposed models incorporating a realized measure, especially when employing the sub-sampled realized variance and the sub-sampled realized range.     

Bayesian modification indices for IRT models

Bayesian modification indices are presented that provide information for the process of model evaluation and model modification. These indices can be used to investigate the improvement in a model if fixed parameters are re-specified as free parameters. The indices can be seen as a Bayesian analogue to the modification indices commonly used in a frequentist framework. The aim is to provide diagnostic information for multi-parameter models where the number of possible model violations and the related number of alternative models is too large to render estimation of each alternative practical. As an example, the method is applied to an item response theory (IRT) model, that is, to the two-parameter model. The method is used to investigate differential item functioning and violations of the assumption of local independence.     

A stochastic frontier analysis to estimate the relative efficiency of Spanish airports

There exists a common belief among researchers and regional policy makers that the actual central system of Aeropuertos Españoles y Navegación Aérea (AENA) should be changed to one more decentralized where airport managers could have more autonomy. The main objective of this article is to evaluate the efficiency of the Spanish airports using Markov chain Monte Carlo (MCMC) simulation to estimate a stochastic frontier analysis (SFA) model. Our results show the existence of a significant level of inefficiency in airport operations. Additionally, we provide efficient marginal cost estimates for each airport which also cast some doubts about the current pricing practices.     

A Model Confidence Set approach to the combination of multivariate volatility forecasts

In predicting conditional covariance matrices of financial portfolios, practitioners are required to choose among several alternative options, facing a number of different sources of uncertainty. A first source is related to the frequency at which prices are observed, either daily or intradaily. Using prices sampled at higher frequency inevitably poses additional sources of uncertainty related to the selection of the optimal intradaily sampling frequency and to the construction of the best realized estimator. Likewise, the choices of model structure and estimation method also have a critical role. In order to alleviate the impact of these sources of uncertainty, we propose a forecast combination strategy based on the Model Confidence Set [MCS] to adaptively identify the set of most accurate predictors. The combined predictor is shown to achieve superior performance with respect to the whole model universe plus three additional competitors, independently of the MCS or portfolio settings.