No Arabic abstract
This paper presents how to apply the stochastic collocation technique to assets that can not move below a boundary. It shows that the polynomial collocation towards a lognormal distribution does not work well. Then, the potentials issues of the related collocated local volatility model (CLV) are explored. Finally, a simple analytical expression for the Dupire local volatility derived from the option prices modelled by stochastic collocation is given.
This paper describes another extension of the Local Variance Gamma model originally proposed by P. Carr in 2008, and then further elaborated on by Carr and Nadtochiy, 2017 (CN2017), and Carr and Itkin, 2018 (CI2018). As compared with the latest version of the model developed in CI2018 and called the ELVG (the Expanded Local Variance Gamma model), here we provide two innovations. First, in all previous papers the model was constructed based on a Gamma time-changed {it arithmetic} Brownian motion: with no drift in CI2017, and with drift in CI2018, and the local variance to be a function of the spot level only. In contrast, here we develop a {it geometric} version of this model with drift. Second, in CN2017 the model was calibrated to option smiles assuming the local variance is a piecewise constant function of strike, while in CI2018 the local variance is a piecewise linear} function of strike. In this paper we consider 3 piecewise linear models: the local variance as a function of strike, the local variance as function of log-strike, and the local volatility as a function of strike (so, the local variance is a piecewise quadratic function of strike). We show that for all these new constructions it is still possible to derive an ordinary differential equation for the option price, which plays a role of Dupires equation for the standard local volatility model, and, moreover, it can be solved in closed form. Finally, similar to CI2018, we show that given multiple smiles the whole local variance/volatility surface can be recovered which does not require solving any optimization problem. Instead, it can be done term-by-term by solving a system of non-linear algebraic equations for each maturity which is fast.
In the classical model of stock prices which is assumed to be Geometric Brownian motion, the drift and the volatility of the prices are held constant. However, in reality, the volatility does vary. In quantitative finance, the Heston model has been successfully used where the volatility is expressed as a stochastic differential equation. In addition, we consider a regime switching model where the stock volatility dynamics depends on an underlying process which is possibly a non-Markov pure jump process. Under this model assumption, we find the locally risk minimizing pricing of European type vanilla options. The price function is shown to satisfy a Heston type PDE.
We study the shapes of the implied volatility when the underlying distribution has an atom at zero and analyse the impact of a mass at zero on at-the-money implied volatility and the overall level of the smile. We further show that the behaviour at small strikes is uniquely determined by the mass of the atom up to high asymptotic order, under mild assumptions on the remaining distribution on the positive real line. We investigate the structural difference with the no-mass-at-zero case, showing how one can--theoretically--distinguish between mass at the origin and a heavy-left-tailed distribution. We numerically test our model-free results in stochastic models with absorption at the boundary, such as the CEV process, and in jump-to-default models. Note that while Lees moment formula tells that implied variance is at most asymptotically linear in log-strike, other celebrated results for exact smile asymptotics such as Benaim and Friz (09) or Gulisashvili (10) do not apply in this setting--essentially due to the breakdown of Put-Call duality.
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 Heston stochastic volatility model is a standard model for valuing financial derivatives, since it can be calibrated using semi-analytical formulas and captures the most basic structure of the market for financial derivatives with simple structure in time-direction. However, extending the model to the case of time-dependent parameters, which would allow for a parametrization of the market at multiple timepoints, proves more challenging. We present a simple and numerically efficient approach to the calibration of the Heston stochastic volatility model with piecewise constant parameters. We show that semi-analytical formulas can also be derived in this more complex case and combine them with recent advances in computational techniques for the Heston model. Our numerical scheme is based on the calculation of the characteristic function using Gauss-Kronrod quadrature with an additional control variate that stabilizes the numerical integrals. We use our method to calibrate the Heston model with piecewise constant parameters to the foreign exchange (FX) options market. Finally, we demonstrate improvements of the Heston model with piecewise constant parameters upon the standard Heston model in selected cases.