Integral Bayesian symbolic regression for optimal discovery of governing equations from scarce and noisy data

TL;DR

This paper introduces an integral Bayesian symbolic regression approach that effectively learns governing differential equations directly from raw, noisy, and scarce time-series data, outperforming existing methods in synthetic and real-world biological applications.

Contribution

It presents a novel integral Bayesian symbolic regression method that bypasses derivative estimation, enabling robust discovery of governing equations from limited and noisy data sets.

Findings

Accurately recovers ground-truth models in synthetic benchmarks.

Makes quasi-optimal predictions across all noise regimes.

Discovers novel microbial growth equations outperforming classical models.

Abstract

Understanding how systems evolve over time often requires discovering the differential equations that govern their behavior. Automatically learning these equations from experimental data is challenging when the data are noisy or limited, and existing approaches struggle, in particular, with the estimation of unobserved derivatives. Here, we introduce an integral Bayesian symbolic regression method that learns governing equations directly from raw time-series data, without requiring manual assumptions or error-prone derivative estimation. By sampling the space of symbolic differential equations and evaluating them via numerical integration, our method robustly identifies governing equations even from noisy or scarce data. We show that this approach accurately recovers ground-truth models in synthetic benchmarks, and that it makes quasi-optimal predictions of system dynamics for all noise regimes. Applying this method to bacterial growth experiments across multiple species and substrates, we discover novel growth equations that outperform classical models in accurately capturing all phases of microbial proliferation, including lag, exponential, and saturation. Unlike standard approaches, our method reveals subtle shifts in growth dynamics, such as double ramp-ups or non-canonical transitions, offering a deeper, data-driven understanding of microbial physiology.