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.

Risk-neutral densities: advanced methods of estimating nonnormal options underlying asset prices and returns

Technical report : Risk-neutral density recovery via spectral analysis

Parsimonious Generative Machine Learning for Non-Gaussian Tail Modeling and Risk-Neutral Distribution Extraction

Risk-Neutral Generative Networks

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

A Non-Gaussian Approach to Risk Measures

Can the Random Walk Model Be Beaten in Out-Of-Sample Density Forecasts? Evidence from Intraday Foreign Exchange Rates

Risk measures with non-Gaussian fluctuations

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

International evidence of the forecasting ability of option‐implied distributions

Accounting for risk of non linear portfolios: a novel Fourier approach

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

Empirical distribution of daily stock returns of selected developing and emerging markets with application to financial risk management

Recovering implied minimum distance risk-neutral probability measure using GMD

Recovering Risk-Neutral Probability Density Functions Using Gaussian Mixture Distributions

Model risk for European-style stock index options

Risk-Neutral Market Simulation

Pricing European Options in Realistic Markets

Financial Stochastic Models Diffusion: From Risk-Neutral to Real-World Measure

Elements for a theory of financial risks

Novel Analytic Representations for Caps, Floors, Collars, and Exchange Options on Continuous Flows, Arbitrage‐Free Relations, and Optimal Investments