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~}.

arXiv comment23 pages, 12 figures, accepted by ApJS

References in corpus (34)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Astropy: A Community Python Package for Astronomy
• emcee: The MCMC Hammer
• Kernel density estimation via diffusion
• The Sloan Digital Sky Survey Quasar Catalog V. Seventh Data Release
• The SDSS Quasar Survey: Quasar Luminosity Function from Data Release Three
• High-Redshift Quasars Found in Sloan Digital Sky Survey Commissioning Data IV: Luminosity Function from the Fall Equatorial Stripe Sampl
• Radio sources in the 6dFGS: Local luminosity functions at 1.4 GHz for star-forming galaxies and radio-loud AGN
• The Luminosity Function of Fermi-detected Flat-Spectrum Radio Quasars
• Evolution of the AGN UV luminosity function from redshift 7.5
• Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). V. Quasar Luminosity Function and Contribution to Cosmic Reionization at z = 6
• A lack of evolution in the very bright-end of the galaxy luminosity function from z = 8-10
• Exploring Reionization-Era Quasars III: Discovery of 16 Quasars at $6.4\lesssim z\lesssim 6.9$ with DESI Legacy Imaging Surveys and UKIRT Hemisphere Survey and Quasar Luminosity Function at $z\sim 6.7$
• First Results from the Lyman Alpha Galaxies in the Epoch of Reionization (LAGER) Survey: Cosmological Reionization at z ~ 7
• Evolution of the galaxy stellar mass functions and UV luminosity functions at $z=6-9$ in the Hubble Frontier Fields
• An improved method of constructing binned luminosity functions
• AGNfitter: A Bayesian MCMC approach to fitting spectral energy distributions of AGN
• The local radio-galaxy population at 20 GHz
• zeus: A Python implementation of Ensemble Slice Sampling for efficient Bayesian parameter inference
• The Ly$\alpha$ Luminosity Function and Cosmic Reionization at $z\sim$ 7.0: a Tale of Two LAGER Fields
• On the redshift cut-off for flat-spectrum radio sources
• The rest-frame UV luminosity function at $z\simeq 4$: a significant contribution of AGN to the bright-end of the galaxy population
• A Flexible Method of Estimating Luminosity Functions
• The HerMES sub-millimetre local and low-redshift luminosity functions
• Constructing a bivariate distribution function with given marginals and correlation: application to the galaxy luminosity function
• MOJAVE: Monitoring of Jets in AGN with VLBA Experiments. IV. The Parent Luminosity Function of Radio-Loud Blazars
• The Infrared Medium-deep Survey. VIII. Quasar Luminosity Function at $z\sim 5$
• Gamma-ray Luminosity and Photon Index Evolution of FSRQ Blazars and Contribution to the Gamma-ray Background
• A Statistical Method for Estimating Luminosity Functions using Truncated Data
• A mixture evolution scenario of the AGN radio luminosity function. II. Do low- and high-power radio-loud AGNs evolve differently?
• A graphical analysis of the systematic error of classical binned methods in constructing luminosity functions
• A flexible method for estimating luminosity functions via Kernel Density Estimation
• A catalogue of faint local radio AGN and the properties of their host galaxies
• What determines the shape of the local (z<0.1) infrared galaxy luminosity function?

Cited by in corpus (1)

• On the cosmic evolution of AGN obscuration and the X-ray luminosity function: XMM-Newton and Chandra spectral analysis of the 31.3 deg${}^2$ Stripe 82X