Static Hedging of Geometric Average Asian Options with Standard Options

In this article, we first establish a theorem that represents the price of an Asian option in terms of standard European options with a shorter term and different strikes. Then using Gauss–Hermite numerical integration, we discretize our theorem so as to use Monte Carlo simulation to examine the error of the static hedging under the Black–Scholes model and the Merton jump-diffusion model. For ease of comparison, we also provide the error of the dynamic hedging. The numerical results show that the static hedging strategy performs better than the dynamic one under both models.     

A new hybrid Monte Carlo simulation for Asian options pricing

The aim of this paper is to present a new hybrid algorithm for pricing financial derivatives in the arithmetic Asian options. In this paper, two variance reduction techniques are combined, the multiple control variates (MCV) and the antithetic variates (AV). We propose an efficient algorithm for pricing arithmetic Asian options based on the AV and the MCV procedures. A detailed numerical study illustrates the efficiency of the proposed algorithm.     

Pricing arithmetic Asian option under a two-factor stochastic volatility model with jumps

This paper presents an efficient Monte Carlo simulation scheme based on the variance reduction methods to evaluate arithmetic average Asian options in the context of the double Heston's stochastic volatility model with jumps. This paper consists of two essential parts. The first part presents a new flexible stochastic volatility model, namely, the double Heston model with jumps. In the second part, by combining two variance reduction procedures via Monte Carlo simulation, we propose an efficient Monte Carlo simulation scheme for pricing arithmetic average Asian options under the double Heston model with jumps. Numerical results illustrate the efficiency of our method.     

PRICING AUSTRALIAN S&P200 OPTIONS: A BAYESIAN APPROACH BASED ON GENERALIZED DISTRIBUTIONAL FORMS

This paper develops a new class of option price models and applies it to options on the Australian S&P200 Index. The class of models generalizes the traditional Black‐Scholes framework by accommodating time‐varying conditional volatility, skewness and excess kurtosis in the underlying returns process. An important property of these more general pricing models is that the computational requirements are essentially the same as those associated with the Black‐Scholes model, with both methods being based on one‐dimensional integrals. Bayesian inferential methods are used to evaluate a range of models nested in the general framework, using observed market option prices. The evaluation is based on posterior parameter distributions, as well as posterior model probabilities. Various fit and predictive measures, plus implied volatility graphs, are also used to rank the alternative models. The empirical results provide evidence that time‐varying volatility, leptokurtosis and a small degree of negative skewness are priced in Australian stock market options.     

基于BP神经网络的S＆P500指数期权定价

期权定价理论源于影响期权价格的变量和期权价格之间的非线性关系,传统的Black-Scholes期权定价公式过于严格的假设削弱了该公式在现实中的适用性,使其在理论与应用上均存在缺陷。因此,能够以任意精度近似复杂非线性系统的神经网络运用于期权定价。分别利用BP神经网络和Black-Scholes期权定价公式对S＆P 500指数看跌期权进行定价,实证结果表明BP神经网络的定价结果要优于Black-Scholes定价公式。     

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.     

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.     

Valuing equity-indexed annuities with icicled barrier options

Inspired by the recent popularity of autocallable structured products, this paper intends to enhance equity-indexed annuities (EIAs) by introducing a new class of barrier options, termed icicled barrier options. The new class of options has a vertical (icicled) barrier along with the horizontal one of the ordinary barrier options, which may act as an additional knock-in or knock-out trigger. To improve the crediting method of EIAs, we propose a new EIA design, termed autocallable EIA, with payoff structure similar to the autocallable products except for the minimum guarantee, and further investigate the possibility of embedding various icicled barrier options into the plain point-to-point or the ratchet EIAs. Explicit pricing formulas for the proposed EIAs and the icicled barrier options are obtained under the Black–Scholes model. To the purpose, we derive the joint distribution of the logarithmic returns at the icicled time and the maturity, and their running maximum. As an application of the well-known reflection principle, the derivation itself is an interesting probability problem and the joint distribution plays a key role in the subsequent pricing stage. Our option pricing result can be easily transferred to EIAs or other equity-linked products. The pricing formulas for the EIAs and the options are illustrated through numerical examples.     

Inference for a Class of Stochastic Volatility Models Using Option and Spot Prices: Application of a Bivariate Kalman Filter

In this paper Bayesian methods are applied to a stochastic volatility model using both the prices of the asset and the prices of options written on the asset. Posterior densities for all model parameters, latent volatilities and the market price of volatility risk are produced via a Markov Chain Monte Carlo (MCMC) sampling algorithm. Candidate draws for the unobserved volatilities are obtained in blocks by applying the Kalman filter and simulation smoother to a linearization of a nonlinear state space representation of the model. Crucially, information from both the spot and option prices affects the draws via the specification of a bivariate measurement equation, with implied Black–Scholes volatilities used to proxy observed option prices in the candidate model. Alternative models nested within the Heston (1993) framework are ranked via posterior odds ratios, as well as via fit, predictive and hedging performance. The method is illustrated using Australian News Corporation spot and option price data.     

Derivation of Kurtosis and Option Pricing Formulas for Popular Volatility Models with Applications in Finance

This article discusses some topics relevant to financial modeling. The kurtosis of a distribution plays an important role in controlling tail-behavior and is used in edgeworth expansion of the call prices. We present derivations of the kurtosis for a number of popular volatility models useful in financial applications, including the class of random coefficient GARCH models. Option pricing formulas for various classes of volatility models are also derived and a simple proof of the option pricing formula under the Black–Scholes model is given.     

Examining the interrelation dynamics between option and stock markets using the Markov-switching vector error correction model

This study examines the dynamics of the interrelation between option and stock markets using the Markov-switching vector error correction model. Specifically, we calculate the implied stock prices from the Black–Scholes 6 Black, F. and Scholes, M. 1973. The pricing of options and corporate liabilities. J. Polit. Econ., 81: 637–659. [Crossref], [Web of Science ®] [Google Scholar] model and establish a statistic framework in which the parameter of the price discrepancy between the observed and implied prices switches according to the phase of the volatility regime. The model is tested in the US S&P 500 stock market. The empirical findings of this work are consistent with the following notions. First, while option markets react more quickly to the newest stock–option disequilibrium shocks than spot markets, as found by earlier studies, we further indicate that the price adjustment process occurring in option markets is pronounced when the high variance condition is concerned, but less so during the stable period. Second, the degree of the co-movement between the observed and implied prices is significantly reduced during the high variance state. Last, the lagged price deviation between the observed and implied prices functions as an indicator of the variance-turning process.     

Pricing power options with a generalized jump diffusion

In this article, the valuation of power option is investigated when the dynamic of the stock price is governed by a generalized jump-diffusion Markov-modulated model. The systematic risk is characterized by the diffusion part, and the non systematic risk is characterized by the pure jump process. The jumps are described by a generalized renewal process with generalized jump amplitude. By introducing NASDAQ Index Model, their risk premium is identified respectively. A risk-neutral measure is identified by employing Esscher transform with two families of parameters, which represent the two parts risk premium. In this article, the non systematic risk premium is considered, based on which the price of power option is studied under the generalized jump-diffusion Markov-modulated model. In the case of a special renewal process with log double exponential jump amplitude, the accurate expressions for the Esscher parameters and the pricing formula are provided. By numerical simulation, the influence of the non systematic risk’s price and the index of the power options on the price of the option is depicted.     

Option pricing for generalized distributions

The Black Scholes formula has been widely used to price financial instruments. The derivation of this formula is based on the assumption of lognormally distributed returns which is often in poor agreement with actual data. An option pricing formula based on the generalized beta of the second kind (GB2) is presented. This formula includes the Black Scholes formula as a special case and accommodates a wide variety of nonlognormally distributed returns. The sensitivity of option values to departures from the skewness and kurtosis associated with the lognormal distribution is investigated.     

An FFT approach for option pricing under a regime-switching stochastic interest rate model

In this article, we investigate the pricing of European-style options under a Markovian regime-switching Hull–White interest rate model. The parameters of this model, including the mean-reversion level, the volatility of the stochastic interest rate, and the volatility of an asset’s value, are modulated by an observable, continuous-time, finite-state Markov chain. A closed-form expression for the characteristic function of the logarithmic terminal asset price is derived. Then, using the fast Fourier transform, a price of a European-style option is computed. In a two-state Markov chain case, numerical examples and empirical studies are presented to illustrate the practical implementation of the model.     

A discrete/continuous choice model on a nonconvex budget set

Decreasing block rate pricing is a nonlinear price system often used for public utility services. Residential gas services in Japan and the United Kingdom are provided under this price schedule. The discrete/continuous choice approach is used to analyze the demand under decreasing block rate pricing. However, the nonlinearity problem, which has not been examined in previous studies, arises because a consumer’s budget set (a set of affordable consumption amounts) is nonconvex, and hence, the resulting model includes highly nonlinear functions. To address this problem, we propose a feasible, efficient method of demand estimation on the nonconvex budget. The advantages of our method are as follows: (i) the construction of an Markov chain Monte Carlo algorithm with an efficient blanket based on the Hermite–Hadamard integral inequality and the power-mean inequality, (ii) the explicit consideration of the (highly nonlinear) separability condition, which often makes numerical likelihood maximization difficult, and (iii) the introduction of normal disturbance into the discrete/continuous choice model on the nonconvex budget set. The proposed method is applied to estimate the Japanese residential gas demand function and evaluate the effect of price schedule changes as a policy experiment.     

The homotopy perturbation method for the Black–Scholes equation

The homotopy perturbation method is designed to obtain a quick and accurate solution to the Black–Scholes equation and boundary conditions for a European option pricing problem. The problem of pricing a European option can be cast a partial differential equation. The analytical solution of the equation is calculated in the form of a convergent power series with easily computable components.     

Exact Estimation of Demand Functions under Block-Rate Pricing

This article proposes an exact estimation of demand functions under block-rate pricing by focusing on increasing block-rate pricing. This is the first study that explicitly considers the separability condition which has been ignored in previous literature. Under this pricing structure, the price changes when consumption exceeds a certain threshold and the consumer faces a utility maximization problem subject to a piecewise-linear budget constraint. Solving this maximization problem leads to a statistical model in which model parameters are strongly restricted by the separability condition. In this article, by taking a hierarchical Bayesian approach, we implement a Markov chain Monte Carlo simulation to properly estimate the demand function. We find, however, that the convergence of the distribution of simulated samples to the posterior distribution is slow, requiring an additional scale transformation step for parameters to the Gibbs sampler. These proposed methods are then applied to estimate the Japanese residential water demand function.     

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.     

An analytical approximation method for pricing barrier options under the double Heston model

Abstract

The purpose of the paper is to provide an efficient pricing method for single barrier options under the double Heston model. By rewriting the model as a singular and regular perturbed BS model, the double Heston model can separately mimic a fast time-scale and a slow time-scale. With the singular and regular perturbation techniques, we analytically derive the first-order asymptotic expansion of the solution to a barrier option pricing partial differential equation. The convergence and efficiency of the approximate method is verified by Monte Carlo simulation. Numerical results show that the presented asymptotic pricing method is fast and accurate.     

Efficient and flexible model-based clustering of jumps in diffusion processes

Jump–diffusion processes involving diffusion processes with discontinuous movements, called jumps, are widely used to model time-series data that commonly exhibit discontinuity in their sample paths. The existing jump–diffusion models have been recently extended to multivariate time-series data. The models are, however, still limited by a single parametric jump-size distribution that is common across different subjects. Such strong parametric assumptions for the shape and structure of a jump-size distribution may be too restrictive and unrealistic for multiple subjects with different characteristics. This paper thus proposes an efficient Bayesian nonparametric method to flexibly model a jump-size distribution while borrowing information across subjects in a clustering procedure using a nested Dirichlet process. For efficient posterior computation, a partially collapsed Gibbs sampler is devised to fit the proposed model. The proposed methodology is illustrated through a simulation study and an application to daily stock price data for companies in the S&P 100 index from June 2007 to June 2017.