Sparse Bayesian mass mapping with uncertainties: local credible intervals
Abstract
ABSTRACT Until recently, mass-mapping techniques for weak gravitational lensing convergence reconstruction have lacked a principled statistical framework upon which to quantify reconstruction uncertainties, without making strong assumptions of Gaussianity. In previous work, we presented a sparse hierarchical Bayesian formalism for convergence reconstruction that addresses this shortcoming. Here, we draw on the concept of local credible intervals (cf. Bayesian error bars) as an extension of the uncertainty quantification techniques previously detailed. These uncertainty quantification techniques are benchmarked against those recovered via Px-MALA – a state-of-the-art proximal Markov chain Monte Carlo (MCMC) algorithm. We find that, typically, our recovered uncertainties are everywhere conservative (never underestimate the uncertainty, yet the approximation error is bounded above), of similar magnitude and highly correlated with those recovered via Px-MALA. Moreover, we demonstrate an increase in computational efficiency of $$\mathcal {O}(10^6)$$ when using our sparse Bayesian approach over MCMC techniques. This computational saving is critical for the application of Bayesian uncertainty quantification to large-scale stage IV surveys such as LSST and Euclid.
- Authors:
-
- Mullard Space Science Laboratory, University College London, Holmbury Hill Rd, Dorking RH5 6NT, UK
- Maxwell Institute for Mathematical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, UK
- Publication Date:
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1581019
- Grant/Contract Number:
- AC02-05CH11231; AC02-76SF00515
- Resource Type:
- Published Article
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Additional Journal Information:
- Journal Name: Monthly Notices of the Royal Astronomical Society Journal Volume: 492 Journal Issue: 1; Journal ID: ISSN 0035-8711
- Publisher:
- Oxford University Press
- Country of Publication:
- United Kingdom
- Language:
- English
Citation Formats
Price, M. A., Cai, X., McEwen, J. D., Pereyra, M., Kitching, T. D., and LSST Dark Energy Science Collaboration. Sparse Bayesian mass mapping with uncertainties: local credible intervals. United Kingdom: N. p., 2019.
Web. doi:10.1093/mnras/stz3453.
Price, M. A., Cai, X., McEwen, J. D., Pereyra, M., Kitching, T. D., & LSST Dark Energy Science Collaboration. Sparse Bayesian mass mapping with uncertainties: local credible intervals. United Kingdom. doi:10.1093/mnras/stz3453.
Price, M. A., Cai, X., McEwen, J. D., Pereyra, M., Kitching, T. D., and LSST Dark Energy Science Collaboration. Tue .
"Sparse Bayesian mass mapping with uncertainties: local credible intervals". United Kingdom. doi:10.1093/mnras/stz3453.
@article{osti_1581019,
title = {Sparse Bayesian mass mapping with uncertainties: local credible intervals},
author = {Price, M. A. and Cai, X. and McEwen, J. D. and Pereyra, M. and Kitching, T. D. and LSST Dark Energy Science Collaboration},
abstractNote = {ABSTRACT Until recently, mass-mapping techniques for weak gravitational lensing convergence reconstruction have lacked a principled statistical framework upon which to quantify reconstruction uncertainties, without making strong assumptions of Gaussianity. In previous work, we presented a sparse hierarchical Bayesian formalism for convergence reconstruction that addresses this shortcoming. Here, we draw on the concept of local credible intervals (cf. Bayesian error bars) as an extension of the uncertainty quantification techniques previously detailed. These uncertainty quantification techniques are benchmarked against those recovered via Px-MALA – a state-of-the-art proximal Markov chain Monte Carlo (MCMC) algorithm. We find that, typically, our recovered uncertainties are everywhere conservative (never underestimate the uncertainty, yet the approximation error is bounded above), of similar magnitude and highly correlated with those recovered via Px-MALA. Moreover, we demonstrate an increase in computational efficiency of $\mathcal {O}(10^6)$ when using our sparse Bayesian approach over MCMC techniques. This computational saving is critical for the application of Bayesian uncertainty quantification to large-scale stage IV surveys such as LSST and Euclid.},
doi = {10.1093/mnras/stz3453},
journal = {Monthly Notices of the Royal Astronomical Society},
number = 1,
volume = 492,
place = {United Kingdom},
year = {2019},
month = {12}
}
DOI: 10.1093/mnras/stz3453
Works referenced in this record:
FAst STatistics for weak Lensing (FASTLens): fast method for weak lensing statistics and map making
journal, May 2009
- Pires, S.; Starck, J. -L.; Amara, A.
- Monthly Notices of the Royal Astronomical Society, Vol. 395, Issue 3
Maximum-a-posteriori estimation with unknown regularisation parameters
conference, August 2015
- Pereyra, Marcelo; Bioucas-Dias, Jose M.; Figueiredo, Mario A. T.
- 2015 23rd European Signal Processing Conference (EUSIPCO)
PRISM: Sparse recovery of the primordial power spectrum
journal, June 2014
- Paykari, P.; Lanusse, F.; Starck, J. -L.
- Astronomy & Astrophysics, Vol. 566
Testing the cosmic shear spatially-flat universe approximation with generalized lensing and shear spectra
journal, July 2018
- Taylor, Peter L.; Kitching, Thomas D.; McEwen, Jason D.
- Physical Review D, Vol. 98, Issue 2
Sparse Reconstruction of the Merging A520 Cluster System
journal, September 2017
- Peel, Austin; Lanusse, François; Starck, Jean-Luc
- The Astrophysical Journal, Vol. 847, Issue 1
Weak lensing: Dark Matter, Dark Energy and Dark Gravity
journal, October 2009
- Heavens, Alan
- Nuclear Physics B - Proceedings Supplements, Vol. 194
An Introduction To Compressive Sampling
journal, March 2008
- Candes, E. J.; Wakin, M. B.
- IEEE Signal Processing Magazine, Vol. 25, Issue 2
Three-Dimensional Reconstruction of the Density Field: an svd Approach to Weak-Lensing Tomography
journal, January 2011
- VanderPlas, J. T.; Connolly, A. J.; Jain, B.
- The Astrophysical Journal, Vol. 727, Issue 2
The integrated bispectrum and beyond
journal, February 2017
- Munshi, Dipak; Coles, Peter
- Journal of Cosmology and Astroparticle Physics, Vol. 2017, Issue 02
HEALPix: A Framework for High‐Resolution Discretization and Fast Analysis of Data Distributed on the Sphere
journal, April 2005
- Gorski, K. M.; Hivon, E.; Banday, A. J.
- The Astrophysical Journal, Vol. 622, Issue 2
Compressed sensing
journal, April 2006
- Donoho, D. L.
- IEEE Transactions on Information Theory, Vol. 52, Issue 4
Maximum-a-Posteriori Estimation with Bayesian Confidence Regions
journal, January 2017
- Pereyra, Marcelo
- SIAM Journal on Imaging Sciences, Vol. 10, Issue 1
Proximal Markov chain Monte Carlo algorithms
journal, May 2015
- Pereyra, Marcelo
- Statistics and Computing, Vol. 26, Issue 4
Uncertainty quantification for radio interferometric imaging: II. MAP estimation
journal, August 2018
- Cai, Xiaohao; Pereyra, Marcelo; McEwen, Jason D.
- Monthly Notices of the Royal Astronomical Society, Vol. 480, Issue 3
Efficient optimal reconstruction of linear fields and band-powers from cosmological data
journal, October 2019
- Horowitz, B.; Seljak, U.; Aslanyan, G.
- Journal of Cosmology and Astroparticle Physics, Vol. 2019, Issue 10
Dark Matter Halos in the Standard Cosmological Model: Results from the Bolshoi Simulation
journal, October 2011
- Klypin, Anatoly A.; Trujillo-Gomez, Sebastian; Primack, Joel
- The Astrophysical Journal, Vol. 740, Issue 2
Hierarchical cosmic shear power spectrum inference
journal, December 2015
- Alsing, Justin; Heavens, Alan; Jaffe, Andrew H.
- Monthly Notices of the Royal Astronomical Society, Vol. 455, Issue 4
Characterizing the nonlinear growth of large-scale structure in the Universe
journal, July 2000
- Coles, Peter; Chiang, Lung-Yih
- Nature, Vol. 406, Issue 6794
Weak gravitational lensing
journal, January 2001
- Bartelmann, Matthias; Schneider, Peter
- Physics Reports, Vol. 340, Issue 4-5
Mapping the dark matter with weak gravitational lensing
journal, February 1993
- Kaiser, Nick; Squires, Gordon
- The Astrophysical Journal, Vol. 404
Dark Energy Survey Year 1 results: curved-sky weak lensing mass map
journal, January 2018
- Chang, C.; Pujol, A.; Mawdsley, B.
- Monthly Notices of the Royal Astronomical Society, Vol. 475, Issue 3
Exponential Convergence of Langevin Distributions and Their Discrete Approximations
journal, December 1996
- Roberts, Gareth O.; Tweedie, Richard L.
- Bernoulli, Vol. 2, Issue 4
Improving weak lensing mass map reconstructions using Gaussian and sparsity priors: application to DES SV
journal, May 2018
- Jeffrey, N.; Abdalla, F. B.; Lahav, O.
- Monthly Notices of the Royal Astronomical Society, Vol. 479, Issue 3
Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging: SARA for radio-interferometric imaging
journal, October 2012
- Carrillo, R. E.; McEwen, J. D.; Wiaux, Y.
- Monthly Notices of the Royal Astronomical Society, Vol. 426, Issue 2
Scalable Bayesian Uncertainty Quantification in Imaging Inverse Problems via Convex Optimization
journal, January 2019
- Repetti, Audrey; Pereyra, Marcelo; Wiaux, Yves
- SIAM Journal on Imaging Sciences, Vol. 12, Issue 1
Cosmology with weak lensing surveys
journal, June 2008
- Munshi, D.; Valageas, P.; Vanwaerbeke, L.
- Physics Reports, Vol. 462, Issue 3
Uncertainty quantification for radio interferometric imaging – I. Proximal MCMC methods
journal, July 2018
- Cai, Xiaohao; Pereyra, Marcelo; McEwen, Jason D.
- Monthly Notices of the Royal Astronomical Society, Vol. 480, Issue 3
High resolution weak lensing mass mapping combining shear and flexion
journal, June 2016
- Lanusse, F.; Starck, J. -L.; Leonard, A.
- Astronomy & Astrophysics, Vol. 591