A flexible method for estimating luminosity functions via Kernel Density Estimation -- II. Generalization and Python implementation

Abstract

We propose a generalization of our previous KDE (kernel density estimation) method for estimating luminosity functions (LFs). This new upgrade further extend the application scope of our KDE method, making it a very flexible approach which is suitable to deal with most of bivariate LF calculation problems. From the mathematical point of view, usually the LF calculation can be abstracted as a density estimation problem in the bounded domain of $\{Z_1<z<Z_2,L>f_{\lim }(z)\}$. We use the transformation-reflection KDE method ($\hat{\varphi }$) to solve the problem, and introduce an approximate method (${\hat{\varphi }}_1$) based on one-dimensional KDE to deal with the small sample size case. In practical applications, the different versions of LF estimators can be flexibly chosen according to the Kolmogorov-Smirnov test criterion. Based on 200 simulated samples, we find that for both cases of dividing or not dividing redshift bins, especially for the latter, our method performs significantly better than the traditional binning method ${\hat{\varphi }}_{\mathrm{bin}}$. Moreover, with the increase of sample size $n$, our LF estimator converges to the true LF remarkably faster than ${\hat{\varphi }}_{\mathrm{bin}}$. To implement our method, we have developed a public, open-source Python Toolkit, called \texttt{kdeLF}. With the support of \texttt{kdeLF}, our KDE method is expected to be a competitive alternative to existing nonparametric estimators, due to its high accuracy and excellent stability. \texttt{kdeLF} is available at \url{http://github.com/yuanzunli/kdeLF} with extensive documentation available at \url{http://kdelf.readthedocs.org/en/latest~}.

23 pages, 12 figures, accepted by ApJS

References in corpus (8)

• Kernel density estimation via diffusion
• Radio sources in the 6dFGS: Local luminosity functions at 1.4 GHz for star-forming galaxies and radio-loud AGN
• A lack of evolution in the very bright-end of the galaxy luminosity function from z = 8-10
• A Flexible Method of Estimating Luminosity Functions
• MOJAVE: Monitoring of Jets in AGN with VLBA Experiments. IV. The Parent Luminosity Function of Radio-Loud Blazars
• Gamma-ray Luminosity and Photon Index Evolution of FSRQ Blazars and Contribution to the Gamma-ray Background
• A flexible method for estimating luminosity functions via Kernel Density Estimation
• What determines the shape of the local (z<0.1) infrared galaxy luminosity function?