Abstract:This paper applies classic linear inverse theory to the estimation of the implied risk neutral probability density function (PDF) from option prices. To overcome non-uniqueness and instability inherent in the option inverse problem, smoothness requirement for the shape of a PDF and an initial model are introduced by a penalty function. Positivity constraints are included as a hard bond on the PDF values. Then the option inverse problem becomes a non-negative least-squares problem which can be solved by the classic methods such as the non-negative least squares program of Lawson and Hanson (1974). The best solution is not the solution that gives best fit, but the solution that gives the optimal trade-off between the goodness of fit and smoothness of the estimated risk natural PDF. The proposed inversion technique is compared to the models of Black-Scholes (BS), a mixture of two lognormals (MLN), Jarrow and Rudd's Edgeworth expansion (JR), and jump diffusion (JD) for the estimation of the PDF from the option prices associated with the September 2007 NYMEX natural gas futures. It is found that the inversion technique not only gives best goodness of fit, but also the significantly better model resolution. BS, JD and MLN models basically cannot resolve the densities far away from the strikes where option prices are observed and can resolve long wavelength features of the densities inside the strikes where option prices are observed. On the other hand, the inversion model can resolve not only the significant details of the densities inside the strikes where option prices are observed, but also the long wavelength features of the densities away from the strikes where option prices are observed. The empirical study for the last three months of the September 2007 futures contract shows that the shapes of the estimated PDFs become more symmetric as the futures contract is closer to the expiration date. The dispersion of the estimated PDFs decreases with decreasing the time to expire, indicating the resolution of uncertainty with time.

Inversion of Option Prices for Implied Risk Neutral Probability Density Functions: General Theory and Its Applications to the Natural Gas Market

Inversion of Option Prices for Implied Risk-Neutral Probability Density Functions: General Theory and Its Applications to the Natural Gas Market

A jump-diffusion model for option pricing under fuzzy environments

Recovering Risk-Neutral Probability Density Functions Using Gaussian Mixture Distributions

A Jump-Diffusion Approach to Modelling Vulnerable Option Pricing

Risk-neutral density extraction from option prices: improved pricing with mixture density networks

Recovering implied minimum distance risk-neutral probability measure using GMD

Shape-Preserving Interpolation and Smoothing for Options Market Implied Volatility

On implied volatility for options—Some reasons to smile and more to correct

Probability Weighting Functions Obtained from Hong Kong Index Option Market

Deep Hedging: Learning Risk-Neutral Implied Volatility Dynamics

Risk-Neutral Generative Networks

High-resolution path-integral development of financial options

Stochastic Path-Dependent Volatility Models for Price-Storage Dynamics in Natural Gas Markets and Discrete-Time Swing Option Pricing

Jump Diffusion-Informed Neural Networks with Transfer Learning for Accurate American Option Pricing under Data Scarcity

Pricing European Options in Realistic Markets

On the implied volatility of Inverse options under stochastic volatility models

Solving the inverse problem of option pricing based on a black-scholes model using a regularized-gauss-newton method

Volatility Parametrizations with Random Coefficients: Analytic Flexibility for Implied Volatility Surfaces

Total variation regularization analysis for inverse volatility option pricing problem

Multivariate Modeling of Natural Gas Spot Trading Hubs Incorporating Futures Market Realized Volatility