Determining and benchmarking risk neutral distributions implied from option prices

Risk neutral probability density functions (RNDs) play a central role in assessing models for stock market behavior. However, it remains challenging to distill a realistic estimate for the RND from empirical data. In this work we introduce a novel method to infer a RND estimate from observed option prices. Our method efficiently yields a realistic rational function approximation to the RND, it is flexible w.r.t. the shape of the underlying distribution and robust in the presence of noise. To show this, we first investigate how well a method can actually retrieve a known distribution from noisy option prices. Then we consider real market data and show how our method can be used to derive a single continuously differentiable RND estimate from empirical call and put option price data.