Analytical and numerical studies on the second-order asymptotic expansion method for European option pricing under two-factor stochastic volatilities

The celebrated Black–Scholes model made the assumption of constant volatility but empirical studies on implied volatility and asset dynamics motivated the use of stochastic volatilities. Christoffersen in 2009 showed that multi-factor stochastic volatilities models capture the asset dynamics more realistically. Fouque in 2012 used it to price European options. In 2013, Chiarella and Ziveyi considered Christoffersen’s ideas and introduced an asset dynamics where the two volatilities of the Heston type act separately and independently on the asset price, and using Fourier transform for the asset price process and double Laplace transform for the two volatilities processes, solved a pricing problem for American options. This paper considers the Chiarella and Ziveyi model and parameterizes it so that the volatilities revert to the long-run-mean with reversion rates that mimic fast (for example daily) and slow (for example seasonal) random effects. Applying asymptotic expansion method presented by Fouque in 2012, we make an extensive and detailed derivation of the approximation prices for European options. We also present numerical studies on the behavior and accuracy of our first- and second-order asymptotic expansion formulas.     

On the arbitrage price of European call options

We show that in a discrete price and discrete time model for option pricing, specifically that given by the Cox–Ross–Rubinstein model, the arbitrage price of a European call option can depend on parameters other than volatility (the standard deviation of the log asset price). We provide two theorems to illustrate this phenomenon. Our first theorem considers two securities with the same volatility so that at a specified time n₀, with probability near 1, the two securities are equal. If their call options differ, both the discounted securities will be martingales. Our second theorem considers two securities with the same volatility so that at times n = 0, ..., N ? 1 the securities are equal with probability near 1. If their call options differ, one of the discounted securities will be a martingale and the other discounted security will be a supermartingale.     

Option pricing under the Merton model of the short rate in subdiffusive Brownian motion regime

In this paper, we propose an extension of the Merton short rate model, which reflects the subdiffusive nature of the short rate dynamics. The subdiffusive property is manifested by the random periods of time, during which the asset price does not change. We derive explicit formulas for European call and put options and present some simulation results for the case of α stable. Moreover, we discuss the implied volatility of this model.     

Forecasting energy futures volatility with threshold augmented heterogeneous autoregressive jump models

Abstract

This study forecasts the volatility of two energy futures markets (oil and gas), using high-frequency data. We, first, disentangle volatility into continuous volatility and jumps. Second, we apply wavelet analysis to study the relationship between volume and the volatility measures for different horizons. Third, we augment the heterogeneous autoregressive (HAR) model by nonlinearly including both jumps and volume. We then propose different empirical extensions of the HAR model. Our study shows that oil and gas volatilities nonlinearly depend on public information (jumps), private information (continuous volatility), and trading volume. Moreover, our threshold augmented HAR model with heterogeneous jumps and continuous volatility outperforms HAR model in forecasting volatility.     

对上海期货交易所金属铜量价关系的实证分析

一、引言对商品价格变动与交易量之间相互关系的分析一直是金融市场讨论的热门话题。Ｋａｒｐｏｆｆ(1987)指出 ,研究商品价格波动与交易量之间关系有助于了解市场的结构、市场信息传播的方式和速度 ,以及市场价格如何对信息作出反映。通常价格的变化可以理解为市场对新的信息的反映 ,而交易量则可解释为投资者新信息认识的差异程度 ,如果商品价格的变动与交易量之间存在关联 ,则可利用价格变动与交易量之间的联合分布信息进行有益的分析和推断。Ｇｒａｎｇｅｒ(196 3)等最早利用周数据资料 ,采用谱分析方法对美国证券交易综合价格指数 (Ｓ…     

New HEAVY Models for Fat-Tailed Realized Covariances and Returns

ABSTRACT

We develop a new score-driven model for the joint dynamics of fat-tailed realized covariance matrix observations and daily returns. The score dynamics for the unobserved true covariance matrix are robust to outliers and incidental large observations in both types of data by assuming a matrix-F distribution for the realized covariance measures and a multivariate Student's t distribution for the daily returns. The filter for the unknown covariance matrix has a computationally efficient matrix formulation, which proves beneficial for estimation and simulation purposes. We formulate parameter restrictions for stationarity and positive definiteness. Our simulation study shows that the new model is able to deal with high-dimensional settings (50 or more) and captures unobserved volatility dynamics even if the model is misspecified. We provide an empirical application to daily equity returns and realized covariance matrices up to 30 dimensions. The model statistically and economically outperforms competing multivariate volatility models out-of-sample. Supplementary materials for this article are available online.     

A wavelet Whittle estimator of generalized long-memory stochastic volatility

We consider a k-GARMA generalization of the long-memory stochastic volatility model, discuss the properties of the model and propose a wavelet-based Whittle estimator for its parameters. Its consistency is shown. Monte Carlo experiments show that the small sample properties are essentially indistinguishable from those of the Whittle estimator, but are favorable with respect to a wavelet-based approximate maximum likelihood estimator. An application is given for the Microsoft Corporation stock, modeling the intraday seasonal patterns of its realized volatility.     

Markov switch smooth transition HYGARCH model: Stability and estimation

Abstract

HYGARCH model is basically used to model long-range dependence in volatility. We propose Markov switch smooth-transition HYGARCH model, where the volatility in each state is a time-dependent convex combination of GARCH and FIGARCH. This model provides a flexible structure to capture different levels of volatilities and also short and long memory effects. The necessary and sufficient condition for the asymptotic stability is derived. Forecast of conditional variance is studied by using all past information through a parsimonious way. Bayesian estimations based on Gibbs sampling are provided. A simulation study has been given to evaluate the estimations and model stability. The competitive performance of the proposed model is shown by comparing it with the HYGARCH and smooth-transition HYGARCH models for some period of the S&P500 and Dow Jones industrial average indices based on volatility and value-at-risk forecasts.     

FORECASTING VOLATILITY IN THE PRESENCE OF MODEL INSTABILITY

Recent advances in financial econometrics have allowed for the construction of efficient ex post measures of daily volatility. This paper investigates the importance of instability in models of realised volatility and their corresponding forecasts. Testing for model instability is conducted with a subsampling method. We show that removing structurally unstable data of a short duration has a negligible impact on the accuracy of conditional mean forecasts of volatility. In contrast, it does provide a substantial improvement in a model's forecast density of volatility. In addition, the forecasting performance improves, often dramatically, when we evaluate models on structurally stable data.     

Volatility,Momentum, and Time-Varying Skewness in Foreign Exchange Returns

This article tests a stochastic volatility model of exchange rates that links both the level of volatility and its instantaneous covariance with returns to pathwise properties of the currency. In particular, the model implies that the return–volatility covariance behaves like a weighted average of recent returns and hence switches signs according to the direction of trends in the data. This implies that the skewness of the finite-horizon return distribution likewise switches sign, leading to time-varying implied volatility “smiles” in options prices. The model is fit and assessed using Bayesian techniques. Some previously reported volatility results are accounted for by the fitted models. The predicted pattern of skewness dynamics accords well with that found in historical options prices.     

Conditional expectation determination based on the J-process using Malliavin calculus applied to pricing American options

The aim of our paper is to elaborate a theoretical methodology based on the Malliavin calculus to calculate the following conditional expectation (P_t(X_t)|(X_s)) for s≤t where the only state variable follows a J-process [Jerbi Y. A new closed-form solution as an extension of the Black—Scholes formula allowing smile curve plotting. Quant Finance. 2013; Online First Article. doi:10.1080/14697688.2012.762458]. The theoretical results are applied to the American option pricing, consisting of an extension of the work of Bally et al. [Pricing and hedging American options by Monte Carlo methods using a Malliavin calculus approach. Monte Carlo Methods Appl. 2005;11-2:97–133], as well as the J-process (with additional parameters λ and θ) is an extension of the Wiener process. The introduction of the aforesaid parameters induces skewness and kurtosis effects, i.e. smile curve allowing to fit with the reality of financial market. In his work Jerbi [Jerbi Y. A new closed-form solution as an extension of the Black–-Scholes formula allowing smile curve plotting. Quant Finance. 2013; Online First Article. doi:10.1080/14697688.2012.762458] showed that the use of the J-process is equivalent to the use of a stochastic volatility model based on the Wiener process as in Heston's. The present work consists on extending this result to the American options. We studied the influence of the parameters λ and θ on the American option price and we find empirical results fitting with the options theory.     

Books on probability ideas are reviewed

In this paper, semiparametric methods are applied to estimate multivariate volatility functions, using a residual approach as in [J. Fan and Q. Yao, Efficient estimation of conditional variance functions in stochastic regression, Biometrika 85 (1998), pp. 645–660; F.A. Ziegelmann, Nonparametric estimation of volatility functions: The local exponential estimator, Econometric Theory 18 (2002), pp. 985–991; F.A. Ziegelmann, A local linear least-absolute-deviations estimator of volatility, Comm. Statist. Simulation Comput. 37 (2008), pp. 1543–1564], among others. Our main goal here is two-fold: (1) describe and implement a number of semiparametric models, such as additive, single-index and (adaptive) functional-coefficient, in volatility estimation, all motivated as alternatives to deal with the curse of dimensionality present in fully nonparametric models; and (2) propose the use of a variation of the traditional cross-validation method to deal with model choice in the class of adaptive functional-coefficient models, choosing simultaneously the bandwidth, the number of covariates in the model and also the single-index smoothing variable. The modified cross-validation algorithm is able to tackle the computational burden caused by the model complexity, providing an important tool in semiparametric volatility estimation. We briefly discuss model identifiability when estimating volatility as well as nonnegativity of the resulting estimators. Furthermore, Monte Carlo simulations for several underlying generating models are implemented and applications to real data are provided.     

Capturing the Spillover Effect With Multiplicative Error Models

Recent statistical models for the analysis of volatility in financial markets serve the purpose of incorporating the effect of other markets in their structure, in order to study the spillover or the contagion phenomena. Extending the Multiplicative Error Model we are able to capture these characteristics, under the assumption that the conditional mean of the volatility can be decomposed into the sum of one component representing the proper volatility of the time series analyzed, and other components, each representing the volatility transmitted from one other market. Each component follows a proper dynamics with elements that can be usefully interpreted. This particular decomposition allows to establish, each time, the contribution brought by each individual market to the global volatility of the market object of the analysis. We experiment this model with four stock indices.     

中国股指期货与股票市场波动性关系的实证分析

使用修正的EGARCH模型与VaR方法检验股指期货的推出对中国股票市场波动性所产生的影响。采用的数据为沪深300指数,样本数据分为股指期货推出前,股指期货推出后的短期、中期和长期与样本数据全体五个时间段。研究表明,从股指期货推出的短期与中期来看,市场对信息的反应比较混乱。从长期来看,股指期货的推出加速了信息的传递速度并且弱化了非对称效应,并没有加大股市的波动性。VaR方法检验表明,股指期货的推出有效降低了股市风险,使A股市场更加成熟和完善。     

A Stochastic Volatility Model With Realized Measures for Option Pricing

Abstract

Based on the fact that realized measures of volatility are affected by measurement errors, we introduce a new family of discrete-time stochastic volatility models having two measurement equations relating both observed returns and realized measures to the latent conditional variance. A semi-analytical option pricing framework is developed for this class of models. In addition, we provide analytical filtering and smoothing recursions for the basic specification of the model, and an effective MCMC algorithm for its richer variants. The empirical analysis shows the effectiveness of filtering and smoothing realized measures in inflating the latent volatility persistence—the crucial parameter in pricing Standard and Poor’s 500 Index options.     

Markov-switching BILINEAR − GARCH models: Structure and estimation

Markov-switching (MS) models are becoming increasingly popular as efficient tools of modeling various phenomena in different disciplines, in particular for non Gaussian time series. In this articlept", we propose a broad class of Markov-switching BILINEAR–GARCH processes (MS ? BLGARCH hereafter) obtained by adding to a MS ? GARCH model one or more interaction components between the observed series and its volatility process. This parameterization offers remarkably rich dynamics and complex behavior for modeling and forecasting financial time-series data which exhibit structural changes. In these models, the parameters of conditional variance are allowed to vary according to some latent time-homogeneous Markov chain with finite state space or “regimes.” The main aim of this new model is to capture asymmetric and hence purported to be able to capture leverage effect characterized by the negativity of the correlation between returns shocks and subsequent shocks in volatility patterns in different regimes. So, first, some basic structural properties of this new model including sufficient conditions ensuring the existence of stationary, causal, ergodic solutions, and moments properties are given. Second, since the second-order structure provides a useful information to identify an appropriate time-series model, we derive the expression of the covariance function of for MS ? BLGARCH and for its powers. As a consequence, we find that the second (resp. higher)-order structure is similar to some linear processes, and hence MS ? BLGARCH (resp. its powers) admit an ARMA representation. This finding allows us for parameter estimation via GMM procedure proved by a Monte Carlo study and applied to foreign exchange rate of the Algerian Dinar against the single European currency.     

Volatility forecasting of financial time series using wavelet based exponential generalized autoregressive conditional heteroscedasticity model

Abstract

An improved forecasting model by merging two different computational models in predicting future volatility was proposed. The model integrates wavelet and EGARCH model where the pre-processing activity based on wavelet transform is performed with de-noising technique to eliminate noise in observed signal. The denoised signal is then feed into EGARCH model to forecast the volatility. The predictive capability of the proposed model is compared with the existing EGARCH model. The results show that the hybrid model has increased the accuracy of forecasting future volatility.     

Empirical Performance and Asset Pricing in Hidden Markov Models

Abstract

To improve the empirical performance of the Black-Scholes model, many alternative models have been proposed to address leptokurtic feature, volatility smile, and volatility clustering effects of the asset return distributions. However, analytical tractability remains a problem for most alternative models. In this article, we study a class of hidden Markov models including Markov switching models and stochastic volatility models, that can incorporate leptokurtic feature, volatility clustering effects, as well as provide analytical solutions to option pricing. We show that these models can generate long memory phenomena when the transition probabilities depend on the time scale. We also provide an explicit analytic formula for the arbitrage-free price of the European options under these models. The issues of statistical estimation and errors in option pricing are also discussed in the Markov switching models.     

Fourier-cosine method for pricing forward starting options with stochastic volatility and jumps

This article provides an efficient method for pricing forward starting options under stochastic volatility model with double exponential jumps. The forward characteristic function of the log asset price is derived and thereby forward starting options are well evaluated by Fourier-cosine technique. Based on adaptive simulated annealing algorithm, the model is calibrated to obtain the estimated parameters. Numerical results show that the pricing method is accurate and fast. Double exponential jumps have pronounced impacts on long-term forward starting options prices. Stochastic volatility model with double exponential jumps fits forward implied volatility smile pretty well in contrast to stochastic volatility model.     

On the invertibility of EGARCH(p,q)

Of the two most widely estimated univariate asymmetric conditional volatility models, the exponential GARCH (or EGARCH) specification is said to be able to capture asymmetry, which refers to the different effects on conditional volatility of positive and negative effects of equal magnitude, and leverage, which refers to the negative correlation between the returns shocks and subsequent shocks to volatility. However, the statistical properties of the (quasi-)maximum likelihood estimator (QMLE) of the EGARCH(p, q) parameters are not available under general conditions, but only for special cases under highly restrictive and unverifiable sufficient conditions, such as EGARCH(1,0) or EGARCH(1,1), and possibly only under simulation. A limitation in the development of asymptotic properties of the QMLE for the EGARCH(p, q) model is the lack of an invertibility condition for the returns shocks underlying the model. It is shown in this article that the EGARCH(p, q) model can be derived from a stochastic process, for which sufficient invertibility conditions can be stated simply and explicitly when the parameters respect a simple condition.¹¹Using the notation introduced in part 2, this refers to the cases where α ≥ |γ| or α ≤ ? |γ|. The first inequality is generally assumed in the literature related to the invertibility of EGARCH. This article provides (in the Appendix) an argument for the possible lack of invertibility when these conditions are not met. This will be useful in reinterpreting the existing properties of the QMLE of the EGARCH(p, q) parameters.