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).     

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.     

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.     

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.     

Option Pricing under the Double Exponential Jump-Diffusion Model with Stochastic Volatility and Interest Rate

This paper proposes an efficient option pricing model that incorporates stochastic interest rate (SIR), stochastic volatility (SV), and double exponential jump into the jump-diffusion settings. The model comprehensively considers the leptokurtosis and heteroscedasticity of the underlying asset’s returns, rare events, and an SIR. Using the model, we deduce the pricing characteristic function and pricing formula of a European option. Then, we develop the Markov chain Monte Carlo method with latent variable to solve the problem of parameter estimation under the double exponential jump-diffusion model with SIR and SV. For verification purposes, we conduct time efficiency analysis, goodness of fit analysis, and jump/drift term analysis of the proposed model. In addition, we compare the pricing accuracy of the proposed model with those of the Black–Scholes and the Kou (2002) models. The empirical results show that the proposed option pricing model has high time efficiency, and the goodness of fit and pricing accuracy are significantly higher than those of the other two models.     

Impact of volatility jumps in a mean-reverting model: Derivative pricing and empirical evidence

This paper investigates the critical role of volatility jumps under mean reversion models. Based on the empirical tests conducted on the historical prices of commodities, we demonstrate that allowing for the presence of jumps in volatility in addition to price jumps is a crucial factor when confronting non-Gaussian return distributions. By employing the particle filtering method, a comparison of results drawn among several mean-reverting models suggests that incorporating volatility jumps ensures an improved fit to the data. We infer further empirical evidence for the existence of volatility jumps from the possible paths of filtered state variables. Our numerical results indicate that volatility jumps significantly affect the level and shape of implied volatility smiles. Finally, we consider the pricing of options under the mean reversion model, where the underlying asset price and its volatility both have jump components.     

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.     

Filtering and identification of Heston's stochastic volatility model and its market risk

We study the filtering problem for the stochastic volatility model of Heston by using the nonlinear estimation theory. To solve the estimation problem for the stochastic volatility process, we use the random time change method. The derived basic equation for the filtering is the so-called Zakai equation and its numerically realized algorithm is proposed with the aid of the splitting-up method. Regarding the European call option problem, the identification of the market price of the volatility risk is also studied. Some numerical simulation studies are demonstrated to show the advantage of the proposed method.     

Calibration of stochastic volatility models: A Tikhonov regularization approach

We aim to calibrate stochastic volatility models from option prices. We develop a Tikhonov regularization approach with an efficient numerical algorithm to recover the risk neutral drift term of the volatility (or variance) process. In contrast to most existing literature, we do not assume that the drift term has any special structure. As such, our algorithm applies to calibration of general stochastic volatility models. An extensive numerical analysis is presented to demonstrate the efficiency of our approach. Interestingly, our empirical study reveals that the risk neutral variance processes recovered from market prices of options on S&P 500 index and EUR/USD exchange rate are indeed linearly mean-reverting.     

Estimators and hypothesis tests for a stochastic frontier function: A Monte Carlo analysis

This paper uses Monte Carlo experimentation to investigate the finite sample properties of the maximum likelihood (ML) and corrected ordinary least squares (COLS) estimators of the half-normal stochastic frontier production function. Results indicate substantial bias in both ML and COLS when the percentage contribution of inefficiency in the composed error (denoted by *) is small, and also that ML should be used in preference to COLS because of large mean square error advantages when * is greater than 50%. The performance of a number of tests of the existence of technical inefficiency is also investigated. The Wald and likelihood ratio (LR) tests are shown to have incorrect size. A one-sided LR test and a test of the significance of the third moment of the OLS residuals are suggested as alternatives, and are shown to have correct size, with the one-sided LR test having the better power of the two.The author would like to thank Bill Griffiths, George Battese, Howard Doran, Bill Greene and two anonymous referees for valuable comments. Any errors which remain are those of the author.     

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.     

Supply Chain Finance: A supply chain-oriented perspective to mitigate commodity risk and pricing volatility

This study addresses the Supply Chain Finance challenge of Commodity Price Volatility (CPV) by adopting a supply chain-oriented perspective. In particular, the effectiveness of two Supply Chain Risk Management (SCRM) strategies in mitigating CPV, namely, Switching suppliers and Substituting Commodities, and the main factors that may affect their value, are investigated with a simulation analysis. A Real Option Valuation (ROV) model was developed and tested on real cases of CPV mitigation, as experienced by a large multinational company (Fortune 100) leader in the Fast Moving Consumer Goods (FMCG) industry. The results show the effectiveness of Switching suppliers and Substituting Commodities in mitigating CPV, highlighting that the convenience of adopting such strategies is strongly influenced by some specific conditions, like the relative values of the long-term prices of the commodities, the purchasing volume, and the sunk cost needed to build flexibility.     

A model for identifying and ranking dangerous accident locations: a case study in Flanders

These days, road safety has become a major concern in most modern societies. In this respect, the determination of road locations that are more dangerous than others (black spots or also called sites with promise) can help in better scheduling road safety policies. The present paper proposes a multivariate model to identify and rank sites according to their total expected cost to the society. Bayesian estimation of the model via a Markov Chain Monte Carlo approach is discussed in this paper. To illustrate the proposed model, accident data from 23,184 accident locations in Flanders (Belgium) are used and a cost function proposed by the European Transport Safety Council is adopted to illustrate the model. It is shown in the paper that the model produces insightful results that can help policy makers in prioritizing road infrastructure investments.     

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.

Forecasting volatility: A reality check based on option pricing, utility function, value-at-risk, and predictive likelihood

We analyze the predictive performance of various volatility models for stock returns. To compare their performance, we choose loss functions for which volatility estimation is of paramount importance. We deal with two economic loss functions (an option pricing function and an utility function) and two statistical loss functions (a goodness-of-fit measure for a value-at-risk (VaR) calculation and a predictive likelihood function). We implement the tests for superior predictive ability of White [Econometrica 68 (5) (2000) 1097] and Hansen [Hansen, P. R. (2001). An unbiased and powerful test for superior predictive ability. Brown University]. We find that, for option pricing, simple models like the Riskmetrics exponentially weighted moving average (EWMA) or a simple moving average, which do not require estimation, perform as well as other more sophisticated specifications. For a utility-based loss function, an asymmetric quadratic GARCH seems to dominate, and this result is robust to different degrees of risk aversion. For a VaR-based loss function, a stochastic volatility model is preferred. Interestingly, the Riskmetrics EWMA model, proposed to calculate VaR, seems to be the worst performer. For the predictive likelihood-based loss function, modeling the conditional standard deviation instead of the variance seems to be a dominant modeling strategy.     

The effect of market sentiment and information asymmetry on option pricing

This work addresses the impact of imperfections, such as information asymmetry and market sentiment, on the performance of option pricing models. More precisely, this work compares the option pricing model of Black and Scholes and the same model in the presence of imperfections. This study is based on S&P 500 options that cover the period between 17/03/2000 and 14/06/2013. The achieved results show that, in general, in the presence of imperfections, the model is more effective than the Black and Scholes model. This research appears to be promising for the incorporation of imperfections into the assessment of options.     

Modeling node heterogeneity in latent space models for multidimensional networks

Multidimensional network data can have different levels of complexity, as nodes may be characterized by heterogeneous individual-specific features, which may vary across the networks. This article introduces a class of models for multidimensional network data, where different levels of heterogeneity within and between networks can be considered. The proposed framework is developed in the family of latent space models, and it aims to distinguish symmetric relations between the nodes and node-specific features. Model parameters are estimated via a Markov Chain Monte Carlo algorithm. Simulated data and an application to a real example, on fruits import/export data, are used to illustrate and comment on the performance of the proposed models.     

Ordinal-response GARCH models for transaction data: A forecasting exercise

We use numerous high-frequency transaction data sets to evaluate the forecasting performances of several dynamic ordinal-response time series models with generalized autoregressive conditional heteroscedasticity (GARCH). The specifications account for three components: leverage effects, in-mean effects and moving average error terms. We estimate the model parameters by developing Markov chain Monte Carlo algorithms. Our empirical analysis shows that the proposed ordinal-response GARCH models achieve better point and density forecasts than standard benchmarks.     

Efficient learning via simulation: A marginalized resample-move approach

In state–space models, parameter learning is practically difficult and is still an open issue. This paper proposes an efficient simulation-based parameter learning method. First, the approach breaks up the interdependence of the hidden states and the static parameters by marginalizing out the states using a particle filter. Second, it applies a Bayesian resample-move approach to this marginalized system. The methodology is generic and needs little design effort. Different from batch estimation methods, it provides posterior quantities necessary for full sequential inference and recursive model monitoring. The algorithm is implemented both on simulated data in a linear Gaussian model for illustration and comparison and on real data in a Lévy jump stochastic volatility model and a structural credit risk model.     

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.