# Quantitative Assessment of Biological Dynamics with Aggregate Data

## TL;DR

This paper introduces a new method for estimating model parameters using aggregate data, improving accuracy in biological modeling.

## Contribution

A novel Bayesian learning framework using modified Hamiltonian Monte Carlo and an elliptical slice sampler for parameter estimation.

## Key findings

- The framework outperforms least-squares fitting in parameter estimation for ODE models.
- It successfully uses real and synthetic data from microbial growth experiments.
- Results show robust data assimilation and improved model accuracy.

## Abstract

We develop and apply a learning framework for parameter estimation in initial value problems that are assessed only indirectly via aggregate data such as sample means and/or standard deviations. Our comprehensive framework follows Bayesian principles and consists of specialized Markov chain Monte Carlo computational schemes that rely on modified Hamiltonian Monte Carlo to align with constraints induced by summary statistics and a novel elliptical slice sampler adapted to the parameters of biological models. We benchmark our methods with synthetic data on microbial growth in batch culture and test them with real growth curve data from laboratory replication experiments on Prochlorococcus microbes. The results indicate that our learning framework can utilize experimental or historical data and lead to robust parameter estimation and data assimilation in ODE models that outperform least-squares fitting.