Detecting episodes of star formation using Bayesian model selection
arXiv:2012.02285 · doi:10.1093/mnras/stab138
Abstract
Bayesian model comparison frameworks can be used when fitting models to data in order to infer the appropriate model complexity in a data-driven manner. We aim to use them to detect the correct number of major episodes of star formation from the analysis of the spectral energy distributions (SEDs) of galaxies, modeled after 3D-HST galaxies at z ~ 1. Starting from the published stellar population properties of these galaxies, we use kernel density estimates to build multivariate input parameter distributions to obtain realistic simulations. We create simulated sets of spectra of varying degrees of complexity (identified by the number of parameters), and derive SED fitting results and evidences for pairs of nested models, including the correct model as well as more simplistic ones, using the BAGPIPES codebase with nested sampling algorithm MultiNest. We then ask the question: is it true - as expected in Bayesian model comparison frameworks - that the correct model has larger evidence?} Our results indicate that the ratio of evidences (the Bayes factor) is able to identify the correct underlying model in the vast majority of cases. The quality of the results improves primarily as a function of the total S/N in the SED. We also compare the Bayes factors obtained using the evidence to those obtained via the Savage-Dickey Density Ratio (SDDR), an analytic approximation which can be calculated using samples from regular Markov Chain Monte Carlo methods. We show that the SDDR ratio can satisfactorily replace a full evidence calculation provided that the sampling density is sufficient.
18 pages, 10 figures
References in corpus (43)
- Scikit-learn: Machine Learning in Python
- The Sloan Digital Sky Survey: Technical Summary
- The Dust Content and Opacity of Actively Star-Forming Galaxies
- MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
- Theory of Star Formation
- dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
- CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey
- The Cosmic Evolution Survey (COSMOS) -- Overview
- Simulating Galaxy Formation with the IllustrisTNG Model
- The propagation of uncertainties in stellar population synthesis modeling I: The relevance of uncertain aspects of stellar evolution and the IMF to the derived physical properties of galaxies
- Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
- Physical Models of Galaxy Formation in a Cosmological Framework
- A simple model to interpret the ultraviolet, optical and infrared emission from galaxies
- Bayes in the sky: Bayesian inference and model selection in cosmology
- Introducing the Illustris Project: the evolution of galaxy populations across cosmic time
- An ultra-deep near-infrared spectrum of a compact quiescent galaxy at z=2.2
- A Comprehensive Analysis of Uncertainties Affecting the Stellar Mass - Halo Mass Relation for 0<z<4
- 3D-HST WFC3-selected Photometric Catalogs in the Five CANDELS/3D-HST Fields: Photometry, Photometric Redshifts and Stellar Masses
- Measuring and Modelling the Redshift Evolution of Clustering: the Hubble Deep Field North
- Modeling the Panchromatic Spectral Energy Distributions of Galaxies
- UltraVISTA: a new ultra-deep near-infrared survey in COSMOS
- Inferring the star-formation histories of massive quiescent galaxies with BAGPIPES: Evidence for multiple quenching mechanisms
- Deriving Physical Properties from Broadband Photometry with Prospector: Description of the Model and a Demonstration of its Accuracy Using 129 Galaxies in the Local Universe
- Modelling and interpreting spectral energy distributions of galaxies with BEAGLE
- A Nested Sampling Algorithm for Cosmological Model Selection
- Fitting the integrated Spectral Energy Distributions of Galaxies
- Applications of Bayesian model selection to cosmological parameters
- On the mass function of star clusters
- An Older, More Quiescent Universe from Panchromatic SED Fitting of the 3D-HST Survey
- A Critical Assessment of Stellar Mass Measurement Methods
- Non-parametric Star Formation History Reconstruction with Gaussian Processes I: Counting Major Episodes of Star Formation
- Breaking the Curve with CANDELS: A Bayesian Approach to Reveal the Non-Universality of the Dust-Attenuation Law at High Redshift
- Efficient Bayesian inference for multimodal problems in cosmology
- Improving the Estimation of Star formation Rates and Stellar Population Ages of High-redshift Galaxies from Broadband Photometry
- Reconstruction of Galaxy Star Formation Histories through SED Fitting: The Dense Basis Approach
- SED fitting with Markov Chain Monte Carlo: Methodology and Application to z=3.1 Lyman Alpha Emitting Galaxies
- On the importance of using appropriate spectral models to derive physical properties of galaxies at 0.7<z<2.8
- BayeSED: A General Approach to Fitting the Spectral Energy Distribution of Galaxies
- Decoding spectral energy distributions of dust-obscured starburst-AGN
- A Comprehensive Bayesian Discrimination of the Simple Stellar Population Model, Star Formation History and Dust Attenuation Law in the Spectral Energy Distribution Modeling of Galaxies
- Simultaneous Estimation of Photometric Redshifts and SED Parameters: Improved Techniques and a Realistic Error Budget
- A hierarchical Bayesian approach for reconstructing the Initial Mass Function of Single Stellar Populations
- Hierarchical Bayesian inference of the Initial Mass Function in Composite Stellar Populations