Estimates of derivatives of (log) densities and related objects

Abstract

We estimate the density and its derivatives using a local polynomial approximation to the logarithm of an unknown density function 𝑓 . The estimator is guaranteed to be nonnegative and achieves the same optimal rate of convergence in the interior as on the boundary of the support of 𝑓 . The estimator is therefore well–suited to applications in which nonnegative density estimates are required, such as in semiparametric maximum likelihood estimation. In addition, we show that our estimator compares favorably with other kernel–based methods, both in terms of asymptotic performance and computational ease. Simulation results confirm that our method can perform similarly or better in finite samples compared to these alternative methods when they are used with optimal inputs, i.e. an Epanechnikov kernel and optimally chosen bandwidth sequence. We provide code in several languages.

• 📚 Paper download
• Code (Julia)
• Code (Python, R, Stata, Matlab)
• 💡 Suggestions for the authors