No Arabic abstract
In this paper we consider the pricing of variable annuities (VAs) with guaranteed minimum withdrawal benefits. We consider two pricing approaches, the classical risk-neutral approach and the benchmark approach, and we examine the associated static and optimal behaviors of both the investor and insurer. The first model considered is the so-called minimal market model, where pricing is achieved using the benchmark approach. The benchmark approach was introduced by Platen in 2001 and has received wide acceptance in the finance community. Under this approach, valuing an asset involves determining the minimum-valued replicating portfolio, with reference to the growth optimal portfolio under the real-world probability measure, and it both subsumes classical risk-neutral pricing as a particular case and extends it to situations where risk-neutral pricing is impossible. The second model is the Black-Scholes model for the equity index, where the pricing of contracts is performed within the risk-neutral framework. Crucially, we demonstrate that when the insurer prices and reserves using the Black-Scholes model, while the insured employs a dynamic withdrawal strategy based on the minimal market model, the insurer may be underestimating the value and associated reserves of the contract.
In electricity markets, it is sensible to use a two-factor model with mean reversion for spot prices. One of the factors is an Ornstein-Uhlenbeck (OU) process driven by a Brownian motion and accounts for the small variations. The other factor is an OU process driven by a pure jump Levy process and models the characteristic spikes observed in such markets. When it comes to pricing, a popular choice of pricing measure is given by the Esscher transform that preserves the probabilistic structure of the driving Levy processes, while changing the levels of mean reversion. Using this choice one can generate stochastic risk premiums (in geometric spot models) but with (deterministically) changing sign. In this paper we introduce a pricing change of measure, which is an extension of the Esscher transform. With this new change of measure we also can slow down the speed of mean reversion and generate stochastic risk premiums with stochastic non constant sign, even in arithmetic spot models. In particular, we can generate risk profiles with positive values in the short end of the forward curve and negative values in the long end. Finally, our pricing measure allows us to have a stationary spot dynamics while still having randomly fluctuating forward prices for contracts far from maturity.
This paper focuses on the pricing of continuous geometric Asian options (GAOs) under a multifactor stochastic volatility model. The model considers fast and slow mean reverting factors of volatility, where slow volatility factor is approximated by a quadratic arc. The asymptotic expansion of the price function is assumed, and the first order price approximation is derived using the perturbation techniques for both floating and fixed strike GAOs. Much simplified pricing formulae for the GAOs are obtained in this multifactor stochastic volatility framework. The zeroth order term in the price approximation is the modified Black-Scholes price for the GAOs. This modified price is expressed in terms of the Black-Scholes price for the GAOs. The accuracy of the approximate option pricing formulae is established, and the model parameter is also estimated by capturing the volatility smiles.
We consider option pricing using a discrete-time Markov switching stochastic volatility with co-jump model, which can model volatility clustering and varying mean-reversion speeds of volatility. For pricing European options, we develop a computationally efficient method for obtaining the probability distribution of average integrated variance (AIV), which is key to option pricing under stochastic-volatility-type models. Building upon the efficiency of the European option pricing approach, we are able to price an American-style option, by converting its pricing into the pricing of a portfolio of European options. Our work also provides constructive guidance for analyzing derivatives based on variance, e.g., the variance swap. Numerical results indicate our methods can be implemented very efficiently and accurately.
The goal of this paper is to investigate the method outlined by one of us (PR) in Cherubini et al. (2009) to compute option prices. We name it the SINC approach. While the COS method by Fang and Osterlee (2009) leverages the Fourier-cosine expansion of truncated densities, the SINC approach builds on the Shannon Sampling Theorem revisited for functions with bounded support. We provide several results which were missing in the early derivation: i) a rigorous proof of the convergence of the SINC formula to the correct option price when the support grows and the number of Fourier frequencies increases; ii) ready to implement formulas for put, Cash-or-Nothing, and Asset-or-Nothing options; iii) a systematic comparison with the COS formula for several log-price models; iv) a numerical challenge against alternative Fast Fourier specifications, such as Carr and Madan (1999) and Lewis (2000); v) an extensive pricing exercise under the rough Heston model of Jaisson and Rosenbaum (2015); vi) formulas to evaluate numerically the moments of a truncated density. The advantages of the SINC approach are numerous. When compared to benchmark methodologies, SINC provides the most accurate and fast pricing computation. The method naturally lends itself to price all options in a smile concurrently by means of Fast Fourier techniques, boosting fast calibration. Pricing requires to resort only to odd moments in the Fourier space. A previous version of this manuscript circulated with the title `Rough Heston: The SINC way.
In this paper we investigate price and Greeks computation of a Guaranteed Minimum Withdrawal Benefit (GMWB) Variable Annuity (VA) when both stochastic volatility and stochastic interest rate are considered together in the Heston Hull-White model. We consider a numerical method the solves the dynamic control problem due to the computing of the optimal withdrawal. Moreover, in order to speed up the computation, we employ Gaussian Process Regression (GPR). Starting from observed prices previously computed for some known combinations of model parameters, it is possible to approximate the whole price function on a defined domain. The regression algorithm consists of algorithm training and evaluation. The first step is the most time demanding, but it needs to be performed only once, while the latter is very fast and it requires to be performed only when predicting the target function. The developed method, as well as for the calculation of prices and Greeks, can also be employed to compute the no-arbitrage fee, which is a common practice in the Variable Annuities sector. Numerical experiments show that the accuracy of the values estimated by GPR is high with very low computational cost. Finally, we stress out that the analysis is carried out for a GMWB annuity but it could be generalized to other insurance products.