Option prices can be used to extract the implied risk-neutral density functions of the future underlying asset prices and returns. These market expectations provide valuable information that can be helpful to policy makers and investors. We tested the accuracy and stability of four nonstructural models in estimating the "true" risk-neutral density functions from option prices: the density functional based on the confluent hypergeometric function, the mixture of lognormal distributions, the smoothed implied volatility smile and Edgeworth expansions. The "true" risk-neutral density is unknown, so it was generated using the Carr–Geman–Madan–Yor (CGMY) Gamma-Ornstein–Uhlenbeck (Gamma-OU) model, a structural model able to generate flexible "true" risk-neutral densities. We observed that the density functional based on the confluent hypergeometric function and mixture of lognormal distributions outperformed the smoothed implied volatility smile and the Edgeworth expansion models in capturing the "true" risk-neutral density. The smoothed implied volatility smile had the best performance in terms of stability. This work aims to be more exhaustive than previous studies, by testing a wider variety of nonstructural models. In addition, we tested the effectiveness of these models in more demanding conditions that are closer to reality. In fact, the use of the CGMY Gamma-OU model to build the "true" risk-neutral density allows densities with a higher probability of extreme events and more leptokurtosis, and also generates sudden, discontinuous moves in prices, which has advantages over the diffusion models used in previous studies.