Search

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., Comment: 18 pages, 10 figures

Cosmic voids and their corresponding redshift-projected mass densities, known as troughs, play an important role in our attempt to model the large-scale structure of the Universe. Understanding these structures enables us to compare the standard model with alternative cosmologies, constrain the dark energy equation of state, and distinguish between different gravitational theories. In this paper, we extend the subspace-constrained mean shift algorithm, a recently introduced method to estimate density ridges, and apply it to 2D weak lensing mass density maps from the Dark Energy Survey Y1 data release to identify curvilinear filamentary structures. We compare the obtained ridges with previous approaches to extract trough structure in the same data, and apply curvelets as an alternative wavelet-based method to constrain densities. We then invoke the Wasserstein distance between noisy and noiseless simulations to validate the denoising capabilities of our method. Our results demonstrate the viability of ridge estimation as a precursor for denoising weak lensing observables to recover the large-scale structure, paving the way for a more versatile and effective search for troughs., Comment: 12 pages, 5 figures, accepted for publication in MNRAS