Analytic Marginalization over Binary Variables in Physics Data

TL;DR

This paper introduces a method to efficiently marginalize over binary correction variables in physics data analysis by leveraging the Ising model, enabling accurate likelihood calculations without exponential complexity.

Contribution

It establishes an exact correspondence between binary marginalization and the Ising model, providing new tools for fast, accurate likelihood evaluation in complex data analyses.

Findings

Exact binary correction marginalization reduces to an Ising model form.

Proposed approximation schemes are computationally efficient.

Application to supernova data shows negligible impact of host-galaxy classification uncertainty.

Abstract

In many data analyses, each measurement may come with a simple yes/no correction; for example, belonging to one of two populations or being contaminated or not. Ignoring such binary effects may bias the results, while accounting for them explicitly quickly becomes infeasible as each of the $N$ data points introduces an additional parameter, resulting in an exponentially growing number of possible configurations ($2^N$). We show that, under generic conditions, an exact treatment of these binary corrections leads to a mathematical form identical to the well-known Ising model from statistical physics. This connection opens up a powerful set of tools developed for the Ising model, enabling fast and accurate likelihood calculations. We present efficient approximation schemes with minimal computational cost and demonstrate their effectiveness in applications, including Type Ia supernova calibration, where we show that the uncertainty in host-galaxy mass classification has negligible impact on the inferred value of the Hubble constant.