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.     

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.     

The tail empirical process for long memory stochastic volatility sequences

This paper describes the limiting behaviour of tail empirical processes associated with long memory stochastic volatility models. We show that such a process has dichotomous behaviour, according to an interplay between the Hurst parameter and the tail index. On the other hand, the tail empirical process with random levels never suffers from long memory. This is very desirable from a practical point of view, since such a process may be used to construct the Hill estimator of the tail index. To prove our results we need to establish new results for regularly varying distributions, which may be of independent interest.     

A semi-analytic pricing formula for lookback options under a general stochastic volatility model

In general, the pricing problems of exotic options in finance do not have analytic solutions under stochastic volatility and so it is hard to compute the option prices or at least it requires much of time to compute them. This paper investigates a semi-analytic pricing method for lookback options in a general stochastic volatility framework. The resultant formula is well connected to the Black–Scholes price that is the first term of a series expansion, which makes computing the option prices relatively efficient. Further, a convergence condition for the expansion is provided with an error bound.     

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.     

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.     

A uniform asymptotic expansion for stochastic volatility model in pricing multi‐asset European options

In this paper, we consider a stochastic volatility model for pricing multi‐asset European options that are widely used in the real world, under the assumption that the volatilities are driven by different OU processes. Using the singular perturbation method for multi‐parameter and the boundary layer theory, we derive a uniform asymptotic expansion for the option prices, as well as the uniform error estimates. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

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.     

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     

Derivative-free Greeks for the Barndorff-Nielsen and Shephard stochastic volatility model

We derive derivative-free formulas for the Delta and other Greeks of options written on an asset modelled by a geometric Brownian motion with stochastic volatility of Barndorff-Nielsen and Shephard type. The method applies the Malliavin calculus in Wiener space which moves differentiation of the payoff function of the option to a random weight function. Our method paves the way for simple Monte Carlo approaches, illustrated by several numerical examples.     

The Heston stochastic volatility model in Hilbert space

We extend the Heston stochastic volatility model to a Hilbert space framework. The tensor Heston stochastic variance process is defined as a tensor product of a Hilbert-valued Ornstein–Uhlenbeck process with itself. The volatility process is then defined by a Cholesky decomposition of the variance process. We define a Hilbert-valued Ornstein–Uhlenbeck process with Wiener noise perturbed by this stochastic volatility, and compute the characteristic functional and covariance operator of this process. This process is then applied to the modeling of forward curves in energy and commodity markets. Finally, we compute the dynamics of the tensor Heston volatility model when the generator is bounded, and study its projection down to the real line for comparison with the classical Heston dynamics.     

Pricing currency options under two-factor Markov-modulated stochastic volatility models

This article investigates the valuation of currency options when the dynamic of the spot Foreign Exchange (FX) rate is governed by a two-factor Markov-modulated stochastic volatility model, with the first stochastic volatility component driven by a lognormal diffusion process and the second independent stochastic volatility component driven by a continuous-time finite-state Markov chain model. The states of the Markov chain can be interpreted as the states of an economy. We employ the regime-switching Esscher transform to determine a martingale pricing measure for valuing currency options under the incomplete market setting. We consider the valuation of the European-style and American-style currency options. In the case of American options, we provide a decomposition result for the American option price into the sum of its European counterpart and the early exercise premium. Numerical results are included.     

Moderate deviations for Euler-Maruyama approximation of Hull-White stochastic volatility model

We consider the Euler-Maruyama discretization of stochastic volatility model dS_t = σ_tS_tdW_t, dσ_t = ωσ_tdZ_t, t ∈ [0, T], which has been widely used in financial practice, where W_t,Z_t, t ∈ [0, T], are two uncorrelated standard Brownian motions. Using asymptotic analysis techniques, the moderate deviation principles for log S_n (or log |S_n| in case S_n is negative) are obtained as n → ∞ under different discretization schemes for the asset price process S_t and the volatility process σ_t. Numerical simulations are presented to compare the convergence speeds in different schemes.     

Analytic solution and estimation of parameters on a stochastic exponential model for a technological diffusion process

In this paper we examine the behaviour of a stochastic model that describes a technological diffusion process (continuously increasing process). Furthermore we obtain a solution for the proposed model through the estimation of the volatility using three different approximations. The adjustment of real data to the final stochastic model confirms its ability of describing and forecasting real cases.     

随机波动率下股价服从O-U过程的期权定价

假设股票价格遵循指数O-U过程,利用随机分析中的鞅方法,得到了具有随机波动率的欧式期权的定价公式,推广了B-S模型.     

An insurance risk model with stochastic volatility

In this paper, we extend the Cramér-Lundberg insurance risk model perturbed by diffusion to incorporate stochastic volatility and study the resulting Gerber-Shiu expected discounted penalty (EDP) function. Under the assumption that volatility is driven by an underlying Ornstein-Uhlenbeck (OU) process, we derive the integro-differential equation which the EDP function satisfies. Not surprisingly, no closed-form solution exists; however, assuming the driving OU process is fast mean-reverting, we apply the singular perturbation theory to obtain an asymptotic expansion of the solution. Two integro-differential equations for the first two terms in this expansion are obtained and explicitly solved. When the claim size distribution is of phase-type, the asymptotic results simplify even further and we succeed in estimating the error of the approximation. Hyper-exponential and mixed-Erlang distributed claims are considered in some detail.     

Valuation of guaranteed annuity options using a stochastic volatility model for equity prices

Guaranteed annuity options are options providing the right to convert a policyholder’s accumulated funds to a life annuity at a fixed rate when the policy matures. These options were a common feature in UK retirement savings contracts issued in the 1970’s and 1980’s when interest rates were high, but caused problems for insurers as the interest rates began to fall in the 1990’s. Currently, these options are frequently sold in the US and Japan as part of variable annuity products. The last decade the literature on pricing and risk management of these options evolved. Until now, for pricing these options generally a geometric Brownian motion for equity prices is assumed. However, given the long maturities of the insurance contracts a stochastic volatility model for equity prices would be more suitable. In this paper explicit expressions are derived for prices of guaranteed annuity options assuming stochastic volatility for equity prices and either a 1-factor or 2-factor Gaussian interest rate model. The results indicate that the impact of ignoring stochastic volatility can be significant.     

Pricing VIX options in a stochastic vol‐of‐vol model

This paper proposes and makes a study of a new model for volatility index option pricing. Factors such as mean‐reversion, jumps, and stochastic volatility are taken into consideration. In particular, the positive volatility skew is addressed by the jump and the stochastic volatility of volatility. Daily calibration is used to check whether the model fits market prices and generates positive volatility skews. Overall, the results show that the mean‐reverting logarithmic jump and stochastic volatility model (called MRLRJSV in the paper) serves as the best model in all the required aspects. Copyright © 2015 John Wiley & Sons, Ltd.