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

TL;DR

This paper introduces an advanced, flexible kernel density estimation method for luminosity functions, implemented in an open-source Python toolkit, outperforming traditional binning methods especially with larger samples.

Contribution

It generalizes previous KDE approaches for luminosity functions, introduces a new approximate method for small samples, and provides a practical Python toolkit for improved estimation accuracy.

Findings

Our method outperforms traditional binning in simulations.

The estimator converges faster to the true LF with increasing sample size.

The Python toolkit facilitates easy application and 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_{\mathrm{lim}}(z) \}$. We use the transformation-reflection KDE method ($\hat{\phi}$) to solve the problem, and introduce an approximate method ($\hat{\phi}_{\mathrm{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{\phi}_{\mathrm{bin}}$. Moreover, with the increase of sample size $n$, our LF estimator converges to the true LF remarkably faster than $\hat{\phi}_{\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~}.