Penalized log-density estimation using Legendre polynomials

Abstract

In this article, we present a penalized log-density estimation method using Legendre polynomials with l(1) penalty to adjust estimate's smoothness. Re-expressing the logarithm of the density estimator via a linear combination of Legendre polynomials, we can estimate parameters by maximizing the penalized log-likelihood function. Besides, we proposed an implementation strategy that builds on the coordinate decent algorithm, together with the Bayesian information criterion (BIC). In particular, we derive a numerical solution to the maximum tuning parameter lambda(max) which leads to all zero coefficients and practically facilitates searching the optimal tuning parameter. Extensive simulation studies clearly show that our proposed estimator is computationally competitive with other existing nonparametric density estimators (e.g., kernel, kernel smooth and logspline estimators) benchmarked by the mean integrated squared errors (MISE) and the mean integrated absolute error (MIAE) under the experiment scenario of separated bimodal models in regard to the true density function. With an application to Old Faithful geyser data, our proposed method is found to effectively perform density estimation.