Building Unbiased Estimators from Non-Gaussian Likelihoods with Application to Shear Estimation

Abstract

We develop a general framework for generating estimators of a given quantity which are unbiased to a given order in the difference between the true value of the underlying quantity and the fiducial position in theory space around which we expand the likelihood. We apply this formalism to rederive the optimal quadratic estimator and show how the replacement of the second derivative matrix with the Fisher matrix is a generic way of creating an unbiased estimator (assuming choice of the fiducial model is independent of data). Next apply the approach to estimation of shear lensing, closely following the work of Bernstein and Armstrong (2014). Our first order estimator reduces to their estimator in the limit of zero shear, but it also naturally allows for the case of non-constant shear and the easy calculation of correlation functions or power spectra using standard methods. Both our first-order estimator and Bernstein and Armstrong's estimator exhibit a bias which is quadratic in true shear. Our third-order estimator is, at least in the realm of the toy problem of Bernstein and Armstrong, unbiased to 0.1% in relative shear errors $\Delta g/|g|$ for shears up to $|g|=0.2$.

arXiv comment8 pages, 3 figures, not too long; v2: matches version accepted by JCAP

References in corpus (7)

• Estimating the Power Spectrum of the Cosmic Microwave Background
• Radical Compression of Cosmic Microwave Background Data
• Bayesian Galaxy Shape Measurement for Weak Lensing Surveys -I. Methodology and a Fast Fitting Algorithm
• Galaxy-Galaxy Flexion: Weak Lensing to Second Order
• Bayesian Lensing Shear Measurement
• An Implementation of Bayesian Lensing Shear Measurement
• Cosmography and Power Spectrum Estimation: a Unified Approach

Cited by in corpus (4)

• An accurate and practical method for inference of weak gravitational lensing from galaxy images
• Cosmology with 6 parameters in the Stage-IV era: efficient marginalisation over nuisance parameters
• Large-scale structure perturbation theory without losing stream crossing
• Renormalization group computation of likelihood functions for cosmological data sets