Estimation of risk-neutral densities using positive convolution approximation

Abstract

This paper proposes a new nonparametric method for estimating the conditional risk-neutral density (RND) from a cross-section of option prices. The idea of the method is to -t option prices by -nding the optimal density in a special admissible set. The admissible set consists of functions, each of which may be represented as a convolution of a positive kernel with another density. The method is termed the positive convolution approximation (PCA). The important properties of PCA are that it (1) is completely agnostic about the data generating process, (2) controls against over-tting while allowing for small samples, (3) always produces arbitragefree estimators, and (4) is computationally simple. In a Monte-Carlo experiment, PCA is compared to several popular methods: mixtures of lognormals (with one, two, and three lognormals), Hermite polynomials, two regularization methods (for the RND and for implied volatilities), and sigma shape polynomials. PCA is found to be a promising alternative to the competitors. c 2003 Elsevier B.V. All rights reserved.