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.

Probability Weighting Functions Obtained from Hong Kong Index Option Market

Pricing and Hedging Foreign Equity Options under Hawkes Jump–diffusion Processes

A jump-diffusion model for option pricing under fuzzy environments

Prediction of Volatility Based on Realized-Garch-kernel-type Models: Evidence from China and the U.S.

The double exponential jump diffusion model for pricing European options under fuzzy environments

On Weighted Lower and Upper Possibilistic Means and Variances of Fuzzy Numbers and Its Application in Decision

Pricing European Currency Options in a Fractional Brownian Environment

Empirical Pricing Kernels Obtained from the Uk Index Options Market

Volatility Risk Premia:Evidence from the Hong Kong Warrant Market

The Variable Precision Method for Elicitation of Probability Weighting Functions.

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

Probability Weighting in Commodity Futures Markets

On the relative performance of some parametric and nonparametric estimators of option prices

Recovering implied minimum distance risk-neutral probability measure using GMD

Recovering Risk-Neutral Probability Density Functions Using Gaussian Mixture Distributions

Asset pricing for the lottery-like security under probability weighting: Based on generalized Wang transform

GWT: A Novel Probability Weighting Function for Stock Price Distortions

Option Pricing under Double Heston Jump-Diffusion Model with Approximative Fractional Stochastic Volatility

Nonparametric Predictive Inference for Asian options

Subjective Probability Adjusted Coherent Financial Risk Measurement Based on General Expected Utility Theory

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