No Arabic abstract
Following Boukai (2021) we present the Generalized Gamma (GG) distribution as a possible RND for modeling European options prices under Hestons (1993) stochastic volatility (SV) model. This distribution is seen as especially useful in situations in which the spots price follows a negatively skewed distribution and hence, Black-Scholes based (i.e. the log-normal distribution) modeling is largely inapt. We apply the GG distribution as RND to modeling current market option data on three large market-index ETFs, namely the SPY, IWM and QQQ as well as on the TLT (an ETF that tracks an index of long term US Treasury bonds). The current option chain of each of the three market-index ETFs shows of a pronounced skew of their volatility `smile which indicates a likely distortion in the Black-Scholes modeling of such option data. Reflective of entirely different market expectations, this distortion appears not to exist in the TLT option data. We provide a thorough modeling of the available option data we have on each ETF (with the October 15, 2021 expiration) based on the GG distribution and compared it to the option pricing and RND modeling obtained directly from a well-calibrated Hestons (1993) SV model (both theoretically and empirically, using Monte-Carlo simulations of the spots price). All three market-index ETFs exhibit negatively skewed distributions which are well-matched with those derived under the GG distribution as RND. The inadequacy of the Black-Scholes modeling in such instances which involve negatively skewed distribution is further illustrated by its impact on the hedging factor, delta, and the immediate implications to the retail trader. In contrast, for the TLT ETF, which exhibits no such distortion to the volatility `smile, the three pricing models (i.e. Hestons, Black-Scholes and Generalized Gamma) appear to yield similar results.
We consider Hestons (1993) stochastic volatility model for valuation of European options to which (semi) closed form solutions are available and are given in terms of characteristic functions. We prove that the class of scale-parameter distributions with mean being the forward spot price satisfies Hestons solution. Thus, we show that any member of this class could be used for the direct risk-neutral valuation of the option price under Hestons SV model. In fact, we also show that any RND with mean being the forward spot price that satisfies Hestons option valuation solution, must be a member of a scale-family of distributions in that mean. As particular examples, we show that one-paramet
This paper proposes the sample path generation method for the stochastic volatility version of CGMY process. We present the Monte-Carlo method for European and American option pricing with the sample path generation and calibrate model parameters to the American style S&P 100 index options market, using the least square regression method. Moreover, we discuss path-dependent options such as Asian and Barrier options.
We propose a general, very fast method to quickly approximate the solution of a parabolic Partial Differential Equation (PDEs) with explicit formulas. Our method also provides equaly fast approximations of the derivatives of the solution, which is a challenge for many other methods. Our approach is based on a computable series expansion in terms of a small parameter. As an example, we treat in detail the important case of the SABR PDE for $beta = 1$, namely $partial_{tau}u = sigma^2 big [ frac{1}{2} (partial^2_xu - partial_xu) + u rho partial_xpartial_sigma u + frac{1}{2} u^2 partial^2_sigma u , big ] + kappa (theta - sigma) partial_sigma$, by choosing $ u$ as small parameter. This yields $u = u_0 + u u_1 + u^2 u_2 + ldots$, with $u_j$ independent of $ u$. The terms $u_j$ are explicitly computable, which is also a challenge for many other, related methods. Truncating this expansion leads to computable approximations of $u$ that are in closed form, and hence can be evaluated very quickly. Most of the other related methods use the time $tau$ as a small parameter. The advantage of our method is that it leads to shorter and hence easier to determine and to generalize formulas. We obtain also an explicit expansion for the implied volatility in the SABR model in terms of $ u$, similar to Hagans formula, but including also the {em mean reverting term.} We provide several numerical tests that show the performance of our method. In particular, we compare our formula to the one due to Hagan. Our results also behave well when used for actual market data and show the mean reverting property of the volatility.
We derive a backward and forward nonlinear PDEs that govern the implied volatility of a contingent claim whenever the latter is well-defined. This would include at least any contingent claim written on a positive stock price whose payoff at a possibly random time is convex. We also discuss suitable initial and boundary conditions for those PDEs. Finally, we demonstrate how to solve them numerically by using an iterative finite-difference approach.
Using agent-based modelling, empirical evidence and physical ideas, such as the energy function and the fact that the phase space must have twice the dimension of the configuration space, we argue that the stochastic differential equations which describe the motion of financial prices with respect to real world probability measures should be of second order (and non-Markovian), instead of first order models `a la Bachelier--Samuelson. Our theoretical result in stochastic dynamical systems shows that one cannot correctly reduce second order models to first order models by simply forgetting about momenta. We propose some simple second order models, including a stochastic constrained n-oscillator, which can explain many market phenomena, such as boom-bust cycles, stochastic quasi-periodic behavior, and hot money going from one market sector to another.