Bayesian evidence estimation from posterior samples with normalizing flows

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Srinivasan, Rahul, Crisostomi, Marco, Trotta, Roberto, Barausse, Enrico, Breschi, Matteo
Formato: Preprint
Publicado: 2024
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866915048654372864
author Srinivasan, Rahul
Crisostomi, Marco
Trotta, Roberto
Barausse, Enrico
Breschi, Matteo
author_facet Srinivasan, Rahul
Crisostomi, Marco
Trotta, Roberto
Barausse, Enrico
Breschi, Matteo
contents We propose a novel method ($floZ$), based on normalizing flows, to estimate the Bayesian evidence (and its numerical uncertainty) from a pre-existing 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 for estimating the evidence: nested sampling (which computes the evidence as its main target) and a $k$-nearest-neighbors technique that produces evidence estimates from posterior samples. Provided representative samples from the target posterior are available, our method is more robust to posterior distributions with sharp features, especially in higher dimensions. For a simple multivariate Gaussian, we demonstrate its accuracy for up to 200 dimensions with $10^5$ posterior samples. $floZ$ has wide applicability, e.g., to estimate evidence from variational inference, Markov Chain Monte Carlo samples, or any other method that delivers samples and their likelihood from the unnormalized posterior density. As a physical application, we use $floZ$ to compute the Bayes factor for the presence of the first overtone in the ringdown signal of the gravitational wave data of GW150914, finding good agreement with nested sampling.
format Preprint
id arxiv_https___arxiv_org_abs_2404_12294
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Bayesian evidence estimation from posterior samples with normalizing flows
Srinivasan, Rahul
Crisostomi, Marco
Trotta, Roberto
Barausse, Enrico
Breschi, Matteo
Machine Learning
Cosmology and Nongalactic Astrophysics
General Relativity and Quantum Cosmology
We propose a novel method ($floZ$), based on normalizing flows, to estimate the Bayesian evidence (and its numerical uncertainty) from a pre-existing 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 for estimating the evidence: nested sampling (which computes the evidence as its main target) and a $k$-nearest-neighbors technique that produces evidence estimates from posterior samples. Provided representative samples from the target posterior are available, our method is more robust to posterior distributions with sharp features, especially in higher dimensions. For a simple multivariate Gaussian, we demonstrate its accuracy for up to 200 dimensions with $10^5$ posterior samples. $floZ$ has wide applicability, e.g., to estimate evidence from variational inference, Markov Chain Monte Carlo samples, or any other method that delivers samples and their likelihood from the unnormalized posterior density. As a physical application, we use $floZ$ to compute the Bayes factor for the presence of the first overtone in the ringdown signal of the gravitational wave data of GW150914, finding good agreement with nested sampling.
title Bayesian evidence estimation from posterior samples with normalizing flows
topic Machine Learning
Cosmology and Nongalactic Astrophysics
General Relativity and Quantum Cosmology
url https://arxiv.org/abs/2404.12294