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

Abstract

Using a fully Bayesian approach, Gaussian Process regression is extended to include marginalisation over the kernel choice and kernel hyperparameters. In addition, Bayesian model comparison via the evidence enables direct kernel comparison. The calculation of the joint posterior was implemented with a transdimensional sampler which simultaneously samples over the discrete kernel choice and their hyperparameters by embedding these in a higher-dimensional space, from which samples are taken using nested sampling. Kernel recovery and mean function inference were explored on synthetic data from exoplanet transit light curve simulations. Subsequently, the method was extended to marginalisation over mean functions and noise models and applied to the inference of the present-day Hubble parameter, $H_0$, from real measurements of the Hubble parameter as a function of redshift, derived from the cosmologically model-independent cosmic chronometer and $\Lambda$CDM-dependent baryon acoustic oscillation observations. The inferred $H_0$ values from the cosmic chronometers, baryon acoustic oscillations and combined datasets are $H_0=66\pm 6\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$, $H_0=67\pm 10\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$ and $H_0=69\pm 6\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$, respectively. The kernel posterior of the cosmic chronometers dataset prefers a non-stationary linear kernel. Finally, the datasets are shown to be not in tension with $\ln R=12.17\pm 0.02$.

arXiv comment17 pages, 11 figures

References in corpus (29)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• In the Realm of the Hubble tension $-$ a Review of Solutions
• Fast and scalable Gaussian process modeling with applications to astronomical time series
• Cosmic Distances Calibrated to 1% Precision with Gaia EDR3 Parallaxes and Hubble Space Telescope Photometry of 75 Milky Way Cepheids Confirm Tension with LambdaCDM
• PolyChord: next-generation nested sampling
• Hubble Parameter and Baryon Acoustic Oscillation Measurement Constraints on the Hubble Constant, the Deviation from the Spatially-Flat $\Lambda$cdm Model, The Deceleration-Acceleration Transition Redshift, and Spatial Curvature
• Dynamic nested sampling: an improved algorithm for parameter estimation and evidence calculation
• A Gaussian process framework for modelling instrumental systematics: application to transmission spectroscopy
• Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets
• $H_0$ from cosmic chronometers and Type Ia supernovae, with Gaussian Processes and the novel Weighted Polynomial Regression method
• PyTransit: Fast and Easy Exoplanet Transit Modelling in Python
• On the evolution of the Hubble constant with the SNe Ia Pantheon Sample and Baryon Acoustic Oscillations: a feasibility study for GRB-cosmology in 2030
• Nested sampling for physical scientists
• Non-parametric spatial curvature inference using late-universe cosmological probes
• Quantifying tensions in cosmological parameters: Interpreting the DES evidence ratio
• Evidence for a Lower Value for $H_0$ from Cosmic Chronometers Data?
• Detecting extrasolar planets from stellar radial velocities using Bayesian evidence
• Nested Sampling Methods
• Bayesian model selection without evidences: application to the dark energy equation-of-state
• Quantifying dimensionality: Bayesian cosmological model complexities
• Constraints on the Dark Side of the Universe and Observational Hubble Parameter Data
• Testing the effect of $H_0$ on $f{\sigma }_8$ tension using a Gaussian Process method
• nestcheck: diagnostic tests for nested sampling calculations
• Bayesian evidence for the tensor-to-scalar ratio $r$ and neutrino masses $m_{\nu }$: Effects of uniform vs logarithmic priors
• Towards a model-independent reconstruction approach for late-time Hubble data
• Kernel Selection for Gaussian Process in Cosmology: with Approximate Bayesian Computation Rejection and Nested Sampling
• Maximum-Entropy Priors with Derived Parameters in a Specified Distribution
• Nested sampling with plateaus
• Bayesian evidence-driven likelihood selection for sky-averaged 21-cm signal extraction

Cited by in corpus (1)

• Resonances in reflective Hamiltonian Monte Carlo