Pricing vulnerable options under a stochastic volatility model

In this paper, we consider the pricing of vulnerable options when the underlying asset follows a stochastic volatility model. We use multiscale asymptotic analysis to derive an analytic approximation formula for the price of the vulnerable options and study the stochastic volatility effect on the option price. A numerical experiment result is presented to demonstrate our findings graphically.     

An efficient ETD method for pricing American options under stochastic volatility with nonsmooth payoffs

Based on Cox and Matthews Exponential Time Differencing (ETD) approach, a fourth–order strongly–stable method having real distinct poles is developed and applied to solve American options under stochastic volatility with nonsmooth payoffs. A computationally efficient version of the method is constructed using partial fraction splitting technique. This approach requires to solve several backward Euler‐type linear systems at each time step. Numerical experiments are presented to demonstrate the computational efficiency, accuracy, and reliability of the method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Two asset-barrier option under stochastic volatility

Financial products which depend on hitting times for two underlying assets have become very popular in the last decade. Three common examples are double-digital barrier options, two-asset barrier spread options and double lookback options. Analytical expressions for the joint distribution of the endpoints and the maximum and/or minimum values of two assets are essential in order to obtain quasi-closed form solutions for the price of these derivatives. Earlier authors derived quasi-closed form pricing expressions in the context of constant volatility and correlation. More recently solutions were provided in the presence of a common stochastic volatility factor but with restricted correlations due to the use of a method of images. In this article, we generalize this finding by allowing any value for the correlation. In this context, we derive closed-form expressions for some two-asset barrier options.     

Regime-switching stochastic volatility model: estimation and calibration to VIX options

We develop and implement a method for maximum likelihood estimation of a regime-switching stochastic volatility model. Our model uses a continuous time stochastic process for the stock dynamics with the instantaneous variance driven by a Cox–Ingersoll–Ross process and each parameter modulated by a hidden Markov chain. We propose an extension of the EM algorithm through the Baum–Welch implementation to estimate our model and filter the hidden state of the Markov chain while using the VIX index to invert the latent volatility state. Using Monte Carlo simulations, we test the convergence of our algorithm and compare it with an approximate likelihood procedure where the volatility state is replaced by the VIX index. We found that our method is more accurate than the approximate procedure. Then, we apply Fourier methods to derive a semi-analytical expression of S&P500 and VIX option prices, which we calibrate to market data. We show that the model is sufficiently rich to encapsulate important features of the joint dynamics of the stock and the volatility and to consistently fit option market prices.     

Unifying pricing formula for several stochastic volatility models with jumps

In this paper, we introduce a unifying approach to option pricing under continuous‐time stochastic volatility models with jumps. For European style options, a new semi‐closed pricing formula is derived using the generalized complex Fourier transform of the corresponding partial integro‐differential equation. This approach is successfully applied to models with different volatility diffusion and jump processes. We also discuss how to price options with different payoff functions in a similar way. In particular, we focus on a log‐normal and a log‐uniform jump diffusion stochastic volatility model, originally introduced by Bates and Yan and Hanson, respectively. The comparison of existing and newly proposed option pricing formulas with respect to time efficiency and precision is discussed. We also derive a representation of an option price under a new approximative fractional jump diffusion model that differs from the aforementioned models, especially for the out‐of‐the money contracts. Copyright © 2017 John Wiley & Sons, Ltd.     

On the Esscher transforms and other equivalent martingale measures for Barndorff-Nielsen and Shephard stochastic volatility models with jumps

We compute and then discuss the Esscher martingale transform for exponential processes, the Esscher martingale transform for linear processes, the minimal martingale measure, the class of structure preserving martingale measures, and the minimum entropy martingale measure for stochastic volatility models of the Ornstein–Uhlenbeck type as introduced by Barndorff-Nielsen and Shephard. We show that in the model with leverage, with jumps both in the volatility and in the returns, all those measures are different, whereas in the model without leverage, with jumps in the volatility only and a continuous return process, several measures coincide, some simplifications can be made and the results are more explicit. We illustrate our results with parametric examples used in the literature.     

On the explicit evaluation of the Geometric Asian options in stochastic volatility models with jumps

In the present paper we provide a semiexplicit valuation formula for Geometric Asian options, with fixed and floating strike under continuous monitoring, when the underlying stock price process exhibits both stochastic volatility and jumps. More precisely, we shall work in the Barndorff-Nielsen and Shephard (BNS) model framework. We shall provide some numerical illustrations of the results obtained.     

Pricing Options Under a Generalized Markov-Modulated Jump-Diffusion Model

Abstract

We consider the pricing of options when the dynamics of the risky underlying asset are driven by a Markov-modulated jump-diffusion model. We suppose that the market interest rate, the drift and the volatility of the underlying risky asset switch over time according to the state of an economy, which is modelled by a continuous-time Markov chain. The measure process is defined to be a generalized mixture of Poisson random measure and encompasses a general class of processes, for example, a generalized gamma process, which includes the weighted gamma process and the inverse Gaussian process. Another interesting feature of the measure process is that jump times and jump sizes can be correlated in general. The model considered here can provide market practitioners with flexibility in modelling the dynamics of the underlying risky asset. We employ the generalized regime-switching Esscher transform to determine an equivalent martingale measure in the incomplete market setting. A system of coupled partial-differential-integral equations satisfied by the European option prices is derived. We also derive a decomposition result for an American put option into its European counterpart and early exercise premium. Simulation results of the model have been presented and discussed.     

Application in stochastic volatility models of nonlinear regression with stochastic design

In regression model with stochastic design, the observations have been primarily treated as a simple random sample from a bivariate distribution. It is of enormous practical significance to generalize the situation to stochastic processes. In this paper, estimation and hypothesis testing problems in stochastic volatility model are considered, when the volatility depends on a nonlinear function of the state variable of other stochastic process, but the correlation coefficient |ρ|≠±1. The methods are applied to estimate the volatility of stock returns from Shanghai stock exchange. Copyright © 2009 John Wiley & Sons, Ltd.     

The implied volatility of Forward-Start options: ATM short-time level,skew and curvature

Using Malliavin Calculus techniques, we derive closed-form expressions for the at-the-money behaviour of the forward implied volatility, its skew and its curvature, in general Markovian stochastic volatility models with continuous paths.     

Efficient numerical methods for pricing American options under stochastic volatility

Five numerical methods for pricing American put options under Heston's stochastic volatility model are described and compared. The option prices are obtained as the solution of a two‐dimensional parabolic partial differential inequality. A finite difference discretization on nonuniform grids leading to linear complementarity problems with M‐matrices is proposed. The projected SOR, a projected multigrid method, an operator splitting method, a penalty method, and a componentwise splitting method are considered. The last one is a direct method while all other methods are iterative. The resulting systems of linear equations in the operator splitting method and in the penalty method are solved using a multigrid method. The projected multigrid method and the componentwise splitting method lead to a sequence of linear complementarity problems with one‐dimensional differential operators that are solved using the Brennan and Schwartz algorithm. The numerical experiments compare the accuracy and speed of the considered methods. The accuracies of all methods appear to be similar. Thus, the additional approximations made in the operator splitting method, in the penalty method, and in the componentwise splitting method do not increase the error essentially. The componentwise splitting method is the fastest one. All multigrid‐based methods have similar rapid grid independent convergence rates. They are about two or three times slower that the componentwise splitting method. On the coarsest grid the speed of the projected SOR is comparable with the multigrid methods while on finer grids it is several times slower. ©John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

American option valuation in a stochastic volatility model with transaction costs

In the present paper we analyse the American option valuation problem in a stochastic volatility model when transaction costs are taken into account. We shall show that it can be formulated as a singular stochastic optimal control problem, proving the existence and uniqueness of the viscosity solution for the associated Hamilton–Jacobi–Bellman partial differential equation. Moreover, after performing a dimensionality reduction through a suitable choice of the utility function, we shall provide a numerical example illustrating how American options prices can be computed in the present modelling framework.     

State-space stochastic volatility models: A review of estimation algorithms

Stochastic volatility models (SVMs) represent an important framework for the analysis of financial time series data, together with ARCH-type models; but unlike the latter, the former, at least from the statistical point of view, cannot rely on the possibility of obtaining exact inference, in particular with regard to maximum likelihood estimates for the parameters of interest. For SVMs, usually only approximate results can be obtained, unless particularly sophisticated estimation strategies like exact non-gaussian filtering methods or simulation techniques are employed. In this paper we review SVM and present a new characterization for them, called ‘generalized bilinear stochastic volatility’. © 1996 John Wiley & Sons, Ltd.     

Multifactor Heston's stochastic volatility model for European option pricing

In this study, we extend the multiscale stochastic volatility model of [Fouque J‐P, Lorig MJ, SIAM J Financial Math. 2011;2(1):221‐254] by incorporating a slow varying factor of volatility. The resulting model can be viewed as a multifactor extension of the Heston model with two additional factors driving the volatility levels. An asymptotic analysis consisting of singular and regular perturbation expansions is developed to obtain an approximation to European option prices. We also find explicit expressions for some essential functions that are available only in integral formulas in the work of [Fouque J‐P, Lorig MJ, SIAM J Financial Math. 2011;2(1):221‐254]. This finding basically leads to considerable reduction in computational time for numerical calculation as well as calibration problems. An accuracy result of the asymptotic approximation is also provided. For numerical illustration, the multifactor Heston model is calibrated to index options on the market, and we find that the resulting implied volatility surfaces fit the market data better than those produced by the multiscale stochastic volatility model of [Fouque J‐P, Lorig MJ, SIAM J Financial Math. 2011;2(1):221‐254], particularly for long‐maturity call options.     

Efficient L-stable method for parabolic problems with application to pricing American options under stochastic volatility

Efficient L-stable numerical method for semilinear parabolic problems with nonsmooth initial data is proposed and implemented to solve Heston’s stochastic volatility model based PDE for pricing American options under stochastic volatility. The proposed new method is also used to solve two asset American options pricing problem. Cox and Matthews [S.M. Cox, P.C. Matthews, Exponential time differencing for stiff systems, Journal of Computational Physics 176 (2002) 430-455] developed a class of exponential time differencing Runge-Kutta schemes (ETDRK) for nonlinear parabolic problems. Kassam and Trefethen [A.K. Kassam, L.N. Trefethen, Fourth-order time stepping for stiff PDEs, SIAM Journal on Scientific Computing 26 (4) (2005) 1214-1233] showed that while computing certain functions involved in the Cox-Matthews schemes, severe cancelation errors can occur which affect the accuracy and stability of the schemes. Kassam and Trefethen proposed complex contour integration technique to implement these schemes in a way that avoids these cancelation errors. But this approach creates new difficulties in choosing and evaluating the contour integrals for larger problems. We modify the ETDRK schemes using positivity preserving Padé approximations of the matrix exponential functions and construct computationally efficient parallel version using splitting technique. As a result of this approach it is required only to solve several backward Euler linear problems in serial or parallel.     

Estimation of integrated volatility in stochastic volatility models

In the framework of stochastic volatility models we examine estimators for the integrated volatility based on the pth power variation (i.e. the sum of pth absolute powers of the log‐returns). We derive consistency and distributional results for the estimators given high‐frequency data, especially taking into account what kind of process we may add to our model without affecting the estimate of the integrated volatility. This may on the one hand be interpreted as a possible flexibility in modelling, for example adding jumps or even leaving the framework of semimartingales by adding a fractional Brownian motion, or on the other hand as robustness against model misspecification. We will discuss possible choices of p under different model assumptions and irregularly spaced data. Copyright © 2005 John Wiley & Sons, Ltd.     

Pricing Volatility Swaps Under Heston's Stochastic Volatility Model with Regime Switching

A model is developed for pricing volatility derivatives, such as variance swaps and volatility swaps under a continuous‐time Markov‐modulated version of the stochastic volatility (SV) model developed by Heston. In particular, it is supposed that the parameters of this version of Heston's SV model depend on the states of a continuous‐time observable Markov chain process, which can be interpreted as the states of an observable macroeconomic factor. The market considered is incomplete in general, and hence, there is more than one equivalent martingale pricing measure. The regime switching Esscher transform used by Elliott et al. is adopted to determine a martingale pricing measure for the valuation of variance and volatility swaps in this incomplete market. Both probabilistic and partial differential equation (PDE) approaches are considered for the valuation of volatility derivatives.     

基于非对称漂移双gamma跳-扩散过程的创新幂型期权定价模型

针对假设股价的对数收益率布朗运动在期权定价时产生的无法解释股价对数收益率的尖峰厚尾和序列相关性的缺陷,采用了Zhang提出的非对称漂移双gamma跳-扩散过程来描述资产(股价)的对数收益率运动形态(该过程是kou提出的双指数跳-扩散过程的推广),并利用Esscher风险中性变换,研究了幂型期权的定价公式.还设计了两种创新的幂型期权,在非对称漂移双gamma跳-扩散过程下给出了相应的定价公式.