Technical report : Risk-neutral density recovery via spectral analysis

In this paper, we propose a new method for estimating the conditional risk-neutral density (RND) directly from a cross-section of put option bid-ask quotes. More precisely, we propose to view the RND recovery problem as an inverse problem. We first show that it is possible to define restricted put and call operators that admit a singular value decomposition (SVD), which we compute explicitly. We subsequently show that this new framework allows us to devise a simple and fast quadratic programming method to recover the smoothest RND whose corresponding put prices lie inside the bid-ask quotes. This method is termed the spectral recovery method (SRM). Interestingly, the SVD of the restricted put and call operators sheds some new light on the RND recovery problem. The SRM improves on other RND recovery methods in the sense that: - it is fast and simple to implement since it requires solution of a single quadratic program, while being fully nonparametric; - it takes the bid ask quotes as sole input and does not require any sort of calibration, smoothing or preprocessing of the data; - it is robust to the paucity of price quotes; - it returns the smoothest density giving rise to prices that lie inside the bid ask quotes. The estimated RND is therefore as well-behaved as can be; - it returns a closed form estimate of the RND on the interval [0,B] of the positive real line, where B is a positive constant that can be chosen arbitrarily. We thus obtain both the middle part of the RND together with its full left tail and part of its right tail. We confront this method to both real and simulated data and observe that it fares well in practice. The SRM is thus found to be a promising alternative to other RND recovery methods.