A Statistical Method for Estimating Luminosity Functions using Truncated Data

Abstract

The observational limitations of astronomical surveys lead to significant statistical inference challenges. One such challenge is the estimation of luminosity functions given redshift $z$ and absolute magnitude $M$ measurements from an irregularly truncated sample of objects. This is a bivariate density estimation problem; we develop here a statistically rigorous method which (1) does not assume a strict parametric form for the bivariate density; (2) does not assume independence between redshift and absolute magnitude (and hence allows evolution of the luminosity function with redshift); (3) does not require dividing the data into arbitrary bins; and (4) naturally incorporates a varying selection function. We accomplish this by decomposing the bivariate density into nonparametric and parametric portions. There is a simple way of estimating the integrated mean squared error of the estimator; smoothing parameters are selected to minimize this quantity. Results are presented from the analysis of a sample of quasars.

30 pages, 9 figures, Accepted for publication in ApJ

References in corpus (1)

• The SDSS Quasar Survey: Quasar Luminosity Function from Data Release Three

Cited by in corpus (16)

• The Demographics of Broad-Line Quasars in the Mass-Luminosity Plane. I. Testing FWHM-Based Virial Black Hole Masses
• Shedding Light on the Galaxy Luminosity Function
• A Flexible Method of Estimating Luminosity Functions
• Determining Quasar Black Hole Mass Functions from their Broad Emission Lines: Application to the Bright Quasar Survey
• 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
• 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
• Constructing a multivariate distribution function with a vine copula: toward multivariate luminosity and mass functions
• The bivariate K-band-submillimetre luminosity functions of the local HRS galaxy sample
• Effect of spectral index distribution on estimating the AGN radio luminosity function
• Completeness III: identifying characteristic systematics and evolution in galaxy redshift surveys
• Cross-disciplinary research in astronomy
• GPU-Accelerated Hierarchical Bayesian Inference with Application to Modeling Cosmic Populations: CUDAHM
• Statistics in astronomy
• A flexible method for estimating luminosity functions via Kernel Density Estimation -- II. Generalization and Python implementation