Abstract:The analysis of photometric large-scale structure data is often complicated by the need to account for many observational and astrophysical systematics. The elaborate models needed to describe them often introduce many ``nuisance parameters'', which can be a major inhibitor of an efficient parameter inference. In this paper we introduce an approximate method to analytically marginalise over a large number of nuisance parameters based on the Laplace approximation. We discuss the mathematics of the method, its relation to concepts such as volume effects and profile likelihood, and show that it can be further simplified for calibratable systematics by linearising the dependence of the theory on the associated parameters. We quantify the accuracy of this approach by comparing it with traditional sampling methods in the context of existing data from the Dark Energy Survey, as well as futuristic Stage-IV photometric data. The linearised version of the method is able to obtain parameter constraints that are virtually equivalent to those found by exploring the full parameter space for a large number of calibratable nuisance parameters, while reducing the computation time by a factor 3-10. Furthermore, the non-linearised approach is able to analytically marginalise over a large number of parameters, returning constraints that are virtually indistinguishable from the brute-force method in most cases, accurately reproducing both the marginalised uncertainty on cosmological parameters, and the impact of volume effects associated with this marginalisation. We provide simple recipes to diagnose when the approximations made by the method fail and one should thus resort to traditional methods. The gains in sampling efficiency associated with this method enable the joint analysis of multiple surveys, typically hindered by the large number of nuisance parameters needed to describe them.

On marginals and profiled posteriors for cosmological parameter estimation

Profile Likelihoods in Cosmology: When, Why and How illustrated with $Λ$CDM, Massive Neutrinos and Dark Energy

Marginal Post Processing of Bayesian Inference Products with Normalizing Flows and Kernel Density Estimators

PROSPECT: A profile likelihood code for frequentist cosmological parameter inference

Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix

Procoli: Profiles of cosmological likelihoods

Cosmology with 6 parameters in the Stage-IV era: efficient marginalisation over nuisance parameters

On the use of marginal posteriors in marginal likelihood estimation via importance-sampling

Improved marginal likelihood estimation via power posteriors and importance sampling

Extremely expensive likelihoods: A variational-Bayes solution for precision cosmology

A comparison of Bayesian sampling algorithms for high-dimensional particle physics and cosmology applications

Understanding posterior projection effects with normalizing flows

The future of cosmological likelihood-based inference: accelerated high-dimensional parameter estimation and model comparison

Informed Total-Error-Minimizing Priors: Interpretable cosmological parameter constraints despite complex nuisance effects

To Profile or To Marginalize -- A SMEFT Case Study

II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements

Cosmological Parameter Inference with Bayesian Statistics

Kernel-, mean- and noise-marginalised Gaussian processes for exoplanet transits and $H_0$ inference

SICRET: Supernova Ia Cosmology with truncated marginal neural Ratio EsTimation

Markov Chain Monte Carlo for Bayesian Parametric Galaxy Modeling in LSST

Bayesian model comparison in cosmology with Population Monte Carlo