Improving exoplanet mass characterisation with Bayesian model selection using the Learned Harmonic Mean Estimator
Abstract
Radial velocity (RV) analyses require modelling choices (such as eccentricity treatment, noise model, velocity trends, and number of planets) that can significantly affect derived planetary masses. Current practice often relies on information criteria to compare and select models, but these have known limitations: they lack the built-in Occam's razor of Bayesian model comparison, and they do not incorporate prior information. Computing the Bayesian evidence needed for Bayes factor model comparison has traditionally required dedicated algorithms such as nested sampling. The learned harmonic mean estimator (LHME) offers an alternative, estimating the Bayesian evidence directly from MCMC posterior samples, with less computational cost and with no modification to the fitting procedure. We present the first application of the LHME to RV model selection, fitting 18 model variants -- comparing circular and eccentric orbits, white noise and Gaussian Process noise models, and long-term velocity trends -- to six single-planet systems, and 72 variants to a seventh system for an N versus N+1 planet model comparison. We find that no single model is universally preferred, reinforcing the need for model comparison to select the most appropriate model for a system, thereby ensuring robust mass characterisation. The LHME, implemented in the open-source harmonic package, makes rigorous Bayesian model comparison accessible to existing MCMC-based RV workflows, and we encourage its wider use for other model comparisons in astrophysics.
Versions
Related
- ROSE (Red-giant Oscillations Spectra Estimator): A Modular Machine-Learning Framework for Automated Asteroseismic Characterisation. I. νmax and Δν from TESS
- Unsupervised selection and characterisation of little red dots in JWST surveys with manifold learning
- Estimator Precision and Selection for Up-the-Ramp Sampling with a High-Read-Noise InGaAs Camera
- Detecting multiple planets in Gaia data using Bayesian model selection
- Astrometric planet search around southern ultracool dwarfs. V. Discovery of the exoplanet candidate PALTA-1 b
BibTeX
@misc{dobson2026improving,
title = {{Improving exoplanet mass characterisation with Bayesian model selection using the Learned Harmonic Mean Estimator}},
author = {Ross S. Dobson and Vincent Van Eylen and Ioanna Manolopoulou and Kendall Sullivan},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2606.27252},
eprint = {2606.27252},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2606.27252},
}
RIS
TY - UNPB
TI - Improving exoplanet mass characterisation with Bayesian model selection using the Learned Harmonic Mean Estimator
AU - Ross S. Dobson
AU - Vincent Van Eylen
AU - Ioanna Manolopoulou
AU - Kendall Sullivan
PY - 2026
DA - 2026/09/28
PB - arXiv
DO - 10.48550/arxiv.2606.27252
UR - https://arxiv.org/abs/2606.27252
AB - Radial velocity (RV) analyses require modelling choices (such as eccentricity treatment, noise model, velocity trends, and number of planets) that can significantly affect derived planetary masses. Current practice often relies on information criteria to compare and select models, but these have known limitations: they lack the built-in Occam's razor of Bayesian model comparison, and they do not incorporate prior information. Computing the Bayesian evidence needed for Bayes factor model comparison has traditionally required dedicated algorithms such as nested sampling. The learned harmonic mean estimator (LHME) offers an alternative, estimating the Bayesian evidence directly from MCMC posterior samples, with less computational cost and with no modification to the fitting procedure. We present the first application of the LHME to RV model selection, fitting 18 model variants -- comparing circular and eccentric orbits, white noise and Gaussian Process noise models, and long-term velocity trends -- to six single-planet systems, and 72 variants to a seventh system for an N versus N+1 planet model comparison. We find that no single model is universally preferred, reinforcing the need for model comparison to select the most appropriate model for a system, thereby ensuring robust mass characterisation. The LHME, implemented in the open-source harmonic package, makes rigorous Bayesian model comparison accessible to existing MCMC-based RV workflows, and we encourage its wider use for other model comparisons in astrophysics.
ER -