American option pricing with discrete and continuous time models: An empirical comparison

This paper considers discrete time GARCH and continuous time SV models and uses these for American option pricing. We first of all show that with a particular choice of framework the parameters of the SV models can be estimated using simple maximum likelihood techniques. We then perform a Monte Carlo study to examine their differences in terms of option pricing, and we study the convergence of the discrete time option prices to their implied continuous time values. Finally, a large scale empirical analysis using individual stock options and options on an index is performed comparing the estimated prices from discrete time models to the corresponding continuous time model prices. The results show that, while the overall differences in performance are small, for the in the money put options on individual stocks the continuous time SV models do generally perform better than the discrete time GARCH specifications.     

SV—N与SV—M模型测量股市波动性的应用

与传统的GARCH类模型一样,SV模犁(随机波动模型)是用来捕捉股市波动特征的一个较好的模型,该模型在国外得到广泛的应用.实证研究表明:利用SV模型的两个子类,即基于正态分布下的SV模型(SV-N)和均值SV模型(SV-M)来测量我国沪深股市波动性明显优于GARCH类模型,能够更好地描述其统计特征.     

Covariance dependent kernels,a Q-affine GARCH for multi-asset option pricing

This paper introduces a class of multivariate GARCH models that extends the existing literature by explicitly modeling correlation dependent pricing kernels. A large subclass admits closed-form recursive solutions for the moment generating function under the risk-neutral measure, which permits efficient pricing of multi-asset options. We perform a full calibration to three bivariate series of index returns and their corresponding volatility indexes in a joint maximum likelihood estimation. The results empirically confirm the presence of correlation dependance in addition to the well known variance dependance in the pricing kernel. The model improves both the overall likelihood and the VIX-implied likelihoods, with a better fitting of marginal distributions, e.g., 15% less error on one-asset option prices. The new degree of freedom is also shown to significantly impact the shape of marginal and joint pricing kernels, and leads to up to 53% differences for out-of-the-money two-asset correlation option prices.     

Forecasting variance using stochastic volatility and GARCH

This paper estimates the conditional variance of daily Swedish OMX-index returns with stochastic volatility (SV) models and GARCH models and evaluates the in-sample performance as well as the out-of-sample forecasting ability of the models. Asymmetric as well as weekend/holiday effects are allowed for in the variance, and the assumption that errors are Gaussian is released. Evidence is found of a leverage effect and of higher variance during weekends. In both in-sample and out-of-sample comparisons SV models outperform GARCH models. However, while asymmetry, weekend/holiday effects and non-Gaussian errors are important for the in-sample fit, it is found that these factors do not contribute to enhancing the forecasting ability of the SV models.     

Some nonstandard stochastic volatility models and their estimation using structured hidden Markov models

We introduce a number of nonstandard stochastic volatility (SV) models and examine their performance when applied to the series of daily returns on several stocks listed on the New York Stock Exchange. The nonstandard models under investigation extend both the observation process and the volatility-generating process of basic SV models. In particular, we consider dependent as well as independent mixtures of autoregressive components as the log-volatility process, and include in the observation equation a lower bound on the volatility. We also consider an experimental SV model that is based on conditionally gamma-distributed volatilities.Our estimation method is based on the fact that an SV model can be approximated arbitrarily accurately by a hidden Markov model (HMM), whose likelihood is easy to compute and to maximize. The method is close, but not identical, to those of Fridman and Harris (1998), Bartolucci and De Luca (2001, 2003) and Clements et al. (2006), and makes explicit the useful link between HMMs and the methods of those authors. Likelihood-based estimation of the parameters of SV models is usually regarded as challenging because the likelihood is a high-dimensional multiple integral. The HMM approximation is easy to implement and particularly convenient for fitting experimental extensions and variants of SV models such as those we introduce here. In addition, and in contrast to the case of SV models themselves, simple formulae are available for the forecast distributions of HMMs, for computing appropriately defined residuals, and for decoding, i.e. estimating the volatility of the process.     

Real-world jump-diffusion term structure models

This paper considers interest rate term structure models in a market attracting both continuous and discrete types of uncertainty. The event-driven noise is modelled by a Poisson random measure. Using as numeraire the growth optimal portfolio, interest rate derivatives are priced under the real-world probability measure. In particular, the real-world dynamics of the forward rates are derived and, for specific volatility structures, finite-dimensional Markovian representations are obtained. Furthermore, allowing for a stochastic short rate in a non-Markovian setting, a class of tractable affine term structures is derived where an equivalent risk-neutral probability measure may not exist.     

Gradient-based simulated maximum likelihood estimation for stochastic volatility models using characteristic functions

Parameter estimation and statistical inference are challenging problems for stochastic volatility (SV) models, especially those driven by pure jump Lévy processes. Maximum likelihood estimation (MLE) is usually preferred when a parametric statistical model is correctly specified, but traditional MLE implementation for SV models is computationally infeasible due to high dimensionality of the integral involved. To overcome this difficulty, we propose a gradient-based simulated MLE method under the hidden Markov structure for SV models, which covers those driven by pure jump Lévy processes. Gradient estimation using characteristic functions and sequential Monte Carlo in the simulation of the hidden states are implemented. Numerical experiments illustrate the efficiency of the proposed method.     

Modeling structural changes in the volatility process

GARCH-type models have been very successful in describing the volatility dynamics of financial return series for short periods of time. However, the time-varying behavior of investors, for example, may cause the structure of volatility to change and the assumption of stationarity is no longer plausible. To deal with this issue, the current paper proposes a conditional volatility model with time-varying coefficients based on a multinomial switching mechanism. By giving more weight to either the persistence or shock term in a GARCH model, conditional on their relative ability to forecast a benchmark volatility measure, the switching reinforces the persistent nature of the GARCH model. The estimation of this benchmark volatility targeting or BVT-GARCH model for Dow 30 stocks indicates that the switching model is able to outperform a number of relevant GARCH setups, both in- and out-of-sample, also without any informational advantages.     

Local volatility function models under a benchmark approach

Without requiring the existence of an equivalent risk-neutral probability measure this paper studies a class of one-factor local volatility function models for stock indices under a benchmark approach. It is assumed that the dynamics for a large diversified index approximates that of the growth optimal portfolio. Fair prices for derivatives when expressed in units of the index are martingales under the real-world probability measure. Different to the classical approach that derives risk-neutral probabilities the paper obtains the transition density for the index with respect to the real-world probability measure. Furthermore, the Dupire formula for the underlying local volatility function is recovered without assuming the existence of an equivalent risk-neutral probability measure. A modification of the constant elasticity of variance model and a version of the minimal market model are discussed as specific examples together with a smoothed local volatility function model that fits a snapshot of S&P500 index options data.     

Estimating GARCH models using support vector machines*

Abstract

Support vector machines (SVMs) are a new nonparametric tool for regression estimation. We will use this tool to estimate the parameters of a GARCH model for predicting the conditional volatility of stock market returns. GARCH models are usually estimated using maximum likelihood (ML) procedures, assuming that the data are normally distributed. In this paper, we will show that GARCH models can be estimated using SVMs and that such estimates have a higher predicting ability than those obtained via common ML methods.     

Modelling conditional heteroscedasticity and jumps in Australian short-term interest rates

The present paper explores a class of jump–diffusion models for the Australian short‐term interest rate. The proposed general model incorporates linear mean‐reverting drift, time‐varying volatility in the form of LEVELS (sensitivity of the volatility to the levels of the short‐rates) and generalized autoregressive conditional heteroscedasticity (GARCH), as well as jumps, to match the salient features of the short‐rate dynamics. Maximum likelihood estimation reveals that pure diffusion models that ignore the jump factor are mis‐specified in the sense that they imply a spuriously high speed of mean‐reversion in the level of short‐rate changes as well as a spuriously high degree of persistence in volatility. Once the jump factor is incorporated, the jump models that can also capture the GARCH‐induced volatility produce reasonable estimates of the speed of mean reversion. The introduction of the jump factor also yields reasonable estimates of the GARCH parameters. Overall, the LEVELS–GARCH–JUMP model fits the data best.     

Forecasting Value-at-Risk using functional volatility incorporating an exogenous effect

This paper proposes a novel extension of log and exponential GARCH models, where time-varying parameters are approximated by orthogonal polynomial systems. These expansions enable us to add and study the effects of market-wide and external international shocks on the volatility forecasts and provide a flexible mechanism to capture various dynamics of the parameters. We examine the performance of the new model in both theoretical and empirical analysis. We investigate the asymptotic properties of the quasi-maximum likelihood estimators under mild conditions. The small-sample behavior of the estimators is studied via Monte Carlo simulation. The performance of the proposed models, in terms of accuracy of both volatility estimation and Value-at-Risk forecasts, is assessed in an empirical study of a set of major stock market indices. The results support the proposed specifications with respect to the corresponding constant-parameters models and to other time-varying parameter models.     

Asian Pacific Stock Market Volatility Modeling and Value at Risk Analysis

The potential for stock market growth in Asian Pacific countries has attracted foreign investors. However, higher growth rates come with higher risk. We apply value at risk (VaR) analysis to measure and analyze stock market index risks in Asian Pacific countries, exposing and detailing both the unique risks and system risks embedded in those markets. To implement the VaR measure, it is necessary to perform "volatility modeling" by mixture switch, exponentially weighted moving average (EWMA), or generalized autoregressive conditional heteroskedasticity (GARCH) models. After estimating the volatility parameters, we can calibrate the VaR values of individual and system risks. Empirically, we find that, on average, Indonesia and Korea exhibit the highest VaRs and VaR sensitivity, and currently, Australia exhibits relatively low values. Taiwan is liable to be in high-state volatility. In addition, the Kupiec test indicates that the mixture switch VaR is superior to delta normal VaR; the quadratic probability score (QPS) shows that the EWMA is inclined to underestimate the VaR for a single series, and GARCH shows no difference from GARCH t and GARCH generalized error distribution (GED) for a multivariate VaR estimate with more assets.     

Estimating stochastic volatility models using integrated nested Laplace approximations

Volatility in financial time series is mainly analysed through two classes of models; the generalized autoregressive conditional heteroscedasticity (GARCH) models and the stochastic volatility (SV) ones. GARCH models are straightforward to estimate using maximum-likelihood techniques, while SV models require more complex inferential and computational tools, such as Markov Chain Monte Carlo (MCMC). Hence, although provided with a series of theoretical advantages, SV models are in practice much less popular than GARCH ones. In this paper, we solve the problem of inference for some SV models by applying a new inferential tool, integrated nested Laplace approximations (INLAs). INLA substitutes MCMC simulations with accurate deterministic approximations, making a full Bayesian analysis of many kinds of SV models extremely fast and accurate. Our hope is that the use of INLA will help SV models to become more appealing to the financial industry, where, due to their complexity, they are rarely used in practice.     

The GARCH Option Pricing Model: A Modification of Lattice Approach

Ritchken and Trevor (1999) proposed a lattice approach for pricing American options under discrete time-varying volatility GARCH frameworks. Even though the lattice approach worked well for the pricing of the GARCH options, it was inappropriate when the option price was computed on the lattice using standard backward recursive procedures, even if the concepts of Cakici and Topyan (2000) were incorporated. This paper shows how to correct the deficiency and that with our adjustment, the lattice method performs properly for option pricing under the GARCH process. JEL Classification: C10, C32, C51, F37, G12     

Bayesian comparison of several continuous time models of the Australian short rate

This paper provides an empirical analysis of a range of alternative single‐factor continuous time models for the Australian short‐term interest rate. The models are nested in a general single‐factor diffusion process for the short rate, with each alternative model indexed by the level effect parameter for the volatility. The inferential approach adopted is Bayesian, with estimation of the models proceeding through a Markov chain Monte Carlo simulation scheme. Discrimination between the alternative models is based on Bayes factors. A data augmentation approach is used to improve the accuracy of the discrete time approximation of the continuous time models. An empirical investigation is conducted using weekly observations on the Australian 90 day interest rate from January 1990 to July 2000. The Bayes factors indicate that the square root diffusion model has the highest posterior probability of all models considered.     

Forecasting volatility of crude oil futures using a GARCH–RNN hybrid approach

Volatility is an important element for various financial instruments owing to its ability to measure the risk and reward value of a given financial asset. Owing to its importance, forecasting volatility has become a critical task in financial forecasting. In this paper, we propose a suite of hybrid models for forecasting volatility of crude oil under different forecasting horizons. Specifically, we combine the parameters of generalized autoregressive conditional heteroscedasticity (GARCH) and Glosten–Jagannathan–Runkle (GJR)-GARCH with long short-term memory (LSTM) to create three new forecasting models named GARCH–LSTM, GJR-LSTM, and GARCH-GJRGARCH LSTM in order to forecast crude oil volatility of West Texas Intermediate on different forecasting horizons and compare their performance with the classical volatility forecasting models. Specifically, we compare the performances against existing methodologies of forecasting volatility such as GARCH and found that the proposed hybrid models improve upon the forecasting accuracy of Crude Oil: West Texas Intermediate under various forecasting horizons and perform better than GARCH and GJR-GARCH, with GG–LSTM performing the best of the three proposed models at 7-, 14-, and 21-day-ahead forecasts in terms of heteroscedasticity-adjusted mean square error and heteroscedasticity-adjusted mean absolute error, with significance testing conducted through the model confidence set showing that GG–LSTM is a strong contender for forecasting crude oil volatility under different forecasting regimes and rolling-window schemes. The contribution of the paper is that it enhances the forecasting ability of crude oil futures volatility, which is essential for trading, hedging, and purposes of arbitrage, and that the proposed model dwells upon existing literature and enhances the forecasting accuracy of crude oil volatility by fusing a neural network model with multiple econometric models.     

Pricing Options under Generalized GARCH and Stochastic Volatility Processes

In this paper, we develop an efficient lattice algorithm to price European and American options under discrete time GARCH processes. We show that this algorithm is easily extended to price options under generalized GARCH processes, with many of the existing stochastic volatility bivariate diffusion models appearing as limiting cases. We establish one unifying algorithm that can price options under almost all existing GARCH specifications as well as under a large family of bivariate diffusions in which volatility follows its own, perhaps correlated, process.     

Identifying Term Structure Volatility from the LIBOR-Swap Curve

This paper proposes a new family of specification tests andapplies them to affine term structure models of the London InterbankOffered Rate (LIBOR)-swap curve. Contrary to Dai and Singleton(2000), the tests show that when standard estimation techniquesare used, affine models do a poor job of forecasting volatilityat the short end of the term structure. Improving the volatilityforecast does not require different models; rather, it requiresa different estimation technique. The paper distinguishes betweentwo econometric procedures for identifying volatility. The "cross-sectional"approach backs out volatility from a cross section of bond yields,and the "time-series" approach imputes volatility from time-seriesvariation in yields. For an affine model, the volatility impliedby the time-series procedure passes the specification tests,while the cross-sectionally identified volatility does not.This is surprising, since under correct specification, the "cross-sectional"approach is maximum likelihood. One explanation is that affinemodels are slightly misspecified; another is that bond yieldsdo not span volatility, as in Collin-Dufresne and Goldstein(2002).