Accelerated Bayesian SED Modeling using Amortized Neural Posterior
Estimation
- URL: http://arxiv.org/abs/2203.07391v1
- Date: Mon, 14 Mar 2022 18:00:03 GMT
- Title: Accelerated Bayesian SED Modeling using Amortized Neural Posterior
Estimation
- Authors: ChangHoon Hahn, Peter Melchior
- Abstract summary: We present an alternative scalable approach to rigorous Bayesian inference using Amortized Neural Posterior Estimation (ANPE)
ANPE is a simulation-based inference method that employs neural networks to estimate the posterior probability distribution.
We present, and publicly release, $rm SEDflow$, an ANPE method to produce posteriors of the recent Hahn et al. (2022) SED model from optical photometry.
- Score: 0.0
- License: http://creativecommons.org/licenses/by/4.0/
- Abstract: State-of-the-art spectral energy distribution (SED) analyses use a Bayesian
framework to infer the physical properties of galaxies from observed photometry
or spectra. They require sampling from a high-dimensional space of SED model
parameters and take $>10-100$ CPU hours per galaxy, which renders them
practically infeasible for analyzing the $billions$ of galaxies that will be
observed by upcoming galaxy surveys ($e.g.$ DESI, PFS, Rubin, Webb, and Roman).
In this work, we present an alternative scalable approach to rigorous Bayesian
inference using Amortized Neural Posterior Estimation (ANPE). ANPE is a
simulation-based inference method that employs neural networks to estimate the
posterior probability distribution over the full range of observations. Once
trained, it requires no additional model evaluations to estimate the posterior.
We present, and publicly release, ${\rm SED}{flow}$, an ANPE method to produce
posteriors of the recent Hahn et al. (2022) SED model from optical photometry.
${\rm SED}{flow}$ takes ${\sim}1$ $second~per~galaxy$ to obtain the posterior
distributions of 12 model parameters, all of which are in excellent agreement
with traditional Markov Chain Monte Carlo sampling results. We also apply ${\rm
SED}{flow}$ to 33,884 galaxies in the NASA-Sloan Atlas and publicly release
their posteriors: see https://changhoonhahn.github.io/SEDflow.
Related papers
- von Mises Quasi-Processes for Bayesian Circular Regression [57.88921637944379]
We explore a family of expressive and interpretable distributions over circle-valued random functions.
The resulting probability model has connections with continuous spin models in statistical physics.
For posterior inference, we introduce a new Stratonovich-like augmentation that lends itself to fast Markov Chain Monte Carlo sampling.
arXiv Detail & Related papers (2024-06-19T01:57:21Z) - Amortizing intractable inference in diffusion models for vision, language, and control [89.65631572949702]
This paper studies amortized sampling of the posterior over data, $mathbfxsim prm post(mathbfx)propto p(mathbfx)r(mathbfx)$, in a model that consists of a diffusion generative model prior $p(mathbfx)$ and a black-box constraint or function $r(mathbfx)$.
We prove the correctness of a data-free learning objective, relative trajectory balance, for training a diffusion model that samples from
arXiv Detail & Related papers (2024-05-31T16:18:46Z) - $floZ$: Improved Bayesian evidence estimation from posterior samples with normalizing flows [0.0]
We introduce $floZ$, an improved method for estimating the Bayesian evidence from a set of samples drawn from the unnormalized posterior distribution.
We validate it on distributions whose evidence is known analytically, up to 15 parameter space dimensions, and compare with two state-of-the-art techniques.
$floZ$ has wide applicability, e.g., to estimate the evidence from variational inference, Markov Chain Monte Carlo samples, or any other method that delivers samples from the unnormalized posterior density.
arXiv Detail & Related papers (2024-04-18T16:16:02Z) - Closed-form Filtering for Non-linear Systems [83.91296397912218]
We propose a new class of filters based on Gaussian PSD Models, which offer several advantages in terms of density approximation and computational efficiency.
We show that filtering can be efficiently performed in closed form when transitions and observations are Gaussian PSD Models.
Our proposed estimator enjoys strong theoretical guarantees, with estimation error that depends on the quality of the approximation and is adaptive to the regularity of the transition probabilities.
arXiv Detail & Related papers (2024-02-15T08:51:49Z) - Dynamical System Identification, Model Selection and Model Uncertainty Quantification by Bayesian Inference [0.8388591755871735]
This study presents a Bayesian maximum textitaposteriori (MAP) framework for dynamical system identification from time-series data.
arXiv Detail & Related papers (2024-01-30T12:16:52Z) - Neural Inference of Gaussian Processes for Time Series Data of Quasars [72.79083473275742]
We introduce a new model that enables it to describe quasar spectra completely.
We also introduce a new method of inference of Gaussian process parameters, which we call $textitNeural Inference$.
The combination of both the CDRW model and Neural Inference significantly outperforms the baseline DRW and MLE.
arXiv Detail & Related papers (2022-11-17T13:01:26Z) - Hierarchical Inference of the Lensing Convergence from Photometric
Catalogs with Bayesian Graph Neural Networks [0.0]
We introduce fluctuations on galaxy-galaxy lensing scales of $sim$1$''$ and extract random sightlines to train our BGNN.
For each test set of 1,000 sightlines, the BGNN infers the individual $kappa$ posteriors, which we combine in a hierarchical Bayesian model.
For a test field well sampled by the training set, the BGNN recovers the population mean of $kappa$ precisely and without bias.
arXiv Detail & Related papers (2022-11-15T00:29:20Z) - Inverting brain grey matter models with likelihood-free inference: a
tool for trustable cytoarchitecture measurements [62.997667081978825]
characterisation of the brain grey matter cytoarchitecture with quantitative sensitivity to soma density and volume remains an unsolved challenge in dMRI.
We propose a new forward model, specifically a new system of equations, requiring a few relatively sparse b-shells.
We then apply modern tools from Bayesian analysis known as likelihood-free inference (LFI) to invert our proposed model.
arXiv Detail & Related papers (2021-11-15T09:08:27Z) - Primordial non-Gaussianity from the Completed SDSS-IV extended Baryon
Oscillation Spectroscopic Survey I: Catalogue Preparation and Systematic
Mitigation [3.2855185490071444]
We investigate the large-scale clustering of the final spectroscopic sample of quasars from the recently completed extended Baryon Oscillation Spectroscopic Survey (eBOSS)
We develop a neural network-based approach to mitigate spurious fluctuations in the density field caused by spatial variations in the quality of the imaging data used to select targets for follow-up spectroscopy.
arXiv Detail & Related papers (2021-06-25T16:01:19Z) - The Generalized Lasso with Nonlinear Observations and Generative Priors [63.541900026673055]
We make the assumption of sub-Gaussian measurements, which is satisfied by a wide range of measurement models.
We show that our result can be extended to the uniform recovery guarantee under the assumption of a so-called local embedding property.
arXiv Detail & Related papers (2020-06-22T16:43:35Z)
This list is automatically generated from the titles and abstracts of the papers in this site.
This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.