On the Application of Bayesian Leave-One-Out Cross-Validation to Exoplanet Atmospheric Analysis

TL;DR

This paper introduces Bayesian leave-one-out cross-validation to evaluate exoplanet atmospheric models at the data point level, enhancing model criticism and reliability assessment for spectroscopic data analysis.

Contribution

It applies Bayesian LOO cross-validation to exoplanet spectra, providing a new, interpretable metric (elpd_LOO) for assessing model performance and data influence at the individual data point level.

Findings

Previous H$^-$ detections depend on a single data point.

elpd_LOO offers a robust out-of-sample predictive accuracy measure.

Method demonstrated on synthetic and real exoplanet spectra.

Abstract

Over the last decade, exoplanetary transmission spectra have yielded an unprecedented understanding about the physical and chemical nature of planets outside our solar system. Physical and chemical knowledge is mainly extracted via fitting competing models to spectroscopic data, based on some goodness-of-fit metric. However, current employed metrics shed little light on how exactly a given model is failing at the individual data point level and where it could be improved. As the quality of our data and complexity of our models increases, there is an urgent need to better understand which observations are driving our model interpretations. Here we present the application of Bayesian leave-one-out cross-validation to assess the performance of exoplanet atmospheric models and compute the expected log pointwise predictive density (elpd$_\text{LOO}$). elpd$_\text{LOO}$ estimates the out-of-sample predictive accuracy of an atmospheric model at data point resolution providing interpretable model criticism. We introduce and demonstrate this method on synthetic HST transmission spectra of a hot Jupiter. We apply elpd$_\text{LOO}$ to interpret current observations of HAT-P-41b and assess the reliability of recent inferences of H$^-$ in its atmosphere. We find that previous detections of H$^{-}$ are dependent solely on a single data point. This new metric for exoplanetary retrievals complements and expands our repertoire of tools to better understand the limits of our models and data. elpd$_\text{LOO}$ provides the means to interrogate models at the single data point level, a prerequisite for robustly interpreting the imminent wealth of spectroscopic information coming from JWST.