Kernel Selection for Gaussian Process in Cosmology: with Approximate Bayesian Computation Rejection and Nested Sampling

Abstract

Gaussian Process (GP) has gained much attention in cosmology due to its ability to reconstruct cosmological data in a model-independent manner. In this study, we compare two methods for GP kernel selection: Approximate Bayesian Computation (ABC) Rejection and nested sampling. We analyze three types of data: cosmic Chronometer data (CC), Type Ia Supernovae (SNIa), and Gamma Ray Burst (GRB), using five kernel functions. To evaluate the differences between kernel functions, we assess the strength of evidence using Bayes factors. Our results show that, for ABC Rejection, the Matérn kernel with $\nu$=5/2 (M52 kernel) outperformes the commonly used Radial Basis Function (RBF) kernel in approximating all three datasets. Bayes factors indicate that the M52 kernel typically supports the observed data better than the RBF kernel, but with no clear advantage over other alternatives. However, nested sampling gives different results, with the M52 kernel losing its advantage. Nevertheless, Bayes factors indicate no significant dependence of the data on each kernel.

arXiv comment15 pages, 4 tables, 6 figures, accepted for publication in ApJS. We draw a new conclusion which is different from the standard ABC-related methods(e.g. ABC-SMC) in kernel selection for Gaussian Process

References in corpus (35)

• Array Programming with NumPy
• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• The Complete Light-curve Sample of Spectroscopically Confirmed Type Ia Supernovae from Pan-STARRS1 and Cosmological Constraints from The Combined Pantheon Sample
• dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
• Constraints on the redshift dependence of the dark energy potential
• Bayes in the sky: Bayesian inference and model selection in cosmology
• A 6% measurement of the Hubble parameter at $z\sim 0.45$: direct evidence of the epoch of cosmic re-acceleration
• Cosmic Chronometers: Constraining the Equation of State of Dark Energy. I: H(z) Measurements
• Raising the bar: new constraints on the Hubble parameter with cosmic chronometers at z$\sim$2
• Improved constraints on the expansion rate of the Universe up to z~1.1 from the spectroscopic evolution of cosmic chronometers
• Adaptive approximate Bayesian computation
• Age-dating Luminous Red Galaxies observed with the Southern African Large Telescope
• PolyChord: nested sampling for cosmology
• Constraints on the equation of state of dark energy and the Hubble constant from stellar ages and the CMB
• PolyChord: next-generation nested sampling
• Reconstruction of dark energy and expansion dynamics using Gaussian processes
• CosmoSIS: modular cosmological parameter estimation
• Gaussian Process Cosmography
• A Possible Cepheid-Like Luminosity Estimator for the Long Gamma-Ray Bursts
• $H_0$ from cosmic chronometers and Type Ia supernovae, with Gaussian Processes and the novel Weighted Polynomial Regression method
• Cosmology with gamma-ray bursts: I. The Hubble diagram through the calibrated $E_{\mathrm{p},\mathrm{i}}$ - $E_{\mathrm{iso}}$ correlation
• Supernova Simulations and Strategies For the Dark Energy Survey
• Non-parametric spatial curvature inference using late-universe cosmological probes
• New Observational $H(z)$ Data from Full-Spectrum Fitting of Cosmic Chronometers in the LEGA-C Survey
• Likelihood-Free Cosmological Inference with Type Ia Supernovae: Approximate Bayesian Computation for a Complete Treatment of Uncertainty
• Approximate Bayesian Computation for Forward Modeling in Cosmology
• cosmoabc: Likelihood-free inference via Population Monte Carlo Approximate Bayesian Computation
• How to use GP: Effects of the mean function and hyperparameter selection on Gaussian Process regression
• astroABC: An Approximate Bayesian Computation Sequential Monte Carlo sampler for cosmological parameter estimation
• Non-parametric reconstruction of the cosmological \textit{jerk} parameter
• An improved method to test the Distance--Duality relation
• Sequential Monte Carlo with Adaptive Weights for Approximate Bayesian Computation
• Towards a model-independent reconstruction approach for late-time Hubble data
• Influence of the Bounds of the Hyperparameters on the Reconstruction of Hubble Constant with Gaussian Process
• Does $\Lambda$CDM really be in tension with the Hubble diagram data?

Cited by in corpus (17)

• The CosmoVerse White Paper: Addressing observational tensions in cosmology with systematics and fundamental physics
• Nonparametric late-time expansion history reconstruction and implications for the Hubble tension in light of recent DESI and type Ia supernovae data
• Testing the cosmic distance duality relation with Type Ia supernova and transverse BAO measurements
• Checking the second law at cosmic scales
• Reconstructing the growth index $\gamma$ with Gaussian Processes
• Multiple measurements on the cosmic curvature using Gaussian process regression without calibration and a cosmological model
• Is ${\omega }_0{\omega }_a$CDM a good model for the clumpy Universe?
• Cosmological parameter constraints using phenomenological symbolic expressions: On the significance of symbolic expression complexity and accuracy
• Probing the Cosmic Distance Duality Relation via Non-Parametric Reconstruction for High Redshifts
• Kernel-, mean- and noise-marginalised Gaussian processes for exoplanet transits and $H_0$ inference
• Testing the Cosmic Distance Duality Relation with Neural Kernel Gaussian Process Regression
• Estimating Cosmological Parameters and Reconstructing Hubble Constant with Artificial Neural Networks: A Test with covariance matrix and mock H(z)
• Model independent calibration for sound horizon combining observations of supernovae and transversal BAO measurements
• Probing the Hubble Tension with an Infinite-Future Condition on the Hubble Parameter
• Constraining the Hubble Constant with a Simulated Full Covariance Matrix Using Neural Networks
• Reconstruction of a dark energy model for the Dirac-Born-Infeld scalar field with the Hubble and DESI data via Gaussian process
• Optimizing Gaussian Process Kernels Using Nested Sampling and ABC Rejection for H(z) Reconstruction