Abstract:Partial differential equation (PDE) models are often used to study biological phenomena involving movement-birth-death processes, including ecological population dynamics and the invasion of populations of biological cells. Count data, by definition, is non-negative, and count data relating to biological populations is often bounded above by some carrying capacity that arises through biological competition for space or nutrients. Parameter estimation, parameter identifiability, and making model predictions usually involves working with a measurement error model that explicitly relating experimental measurements with the solution of a mathematical model. In many biological applications, a typical approach is to assume the data are normally distributed about the solution of the mathematical model. Despite the widespread use of the standard additive Gaussian measurement error model, the assumptions inherent in this approach are rarely explicitly considered or compared with other options. Here, we interpret scratch assay data, involving migration, proliferation and delays in a population of cancer cells using a reaction-diffusion PDE model. We consider relating experimental measurements to the PDE solution using a standard additive Gaussian measurement error model alongside a comparison to a more biologically realistic binomial measurement error model. While estimates of model parameters are relatively insensitive to the choice of measurement error model, model predictions for data realisations are very sensitive. The standard additive Gaussian measurement error model leads to biologically inconsistent predictions, such as negative counts and counts that exceed the carrying capacity across a relatively large spatial region within the experiment. Furthermore, the standard additive Gaussian measurement error model requires estimating an additional parameter compared to the binomial measurement error model. In contrast, the binomial measurement error model leads to biologically plausible predictions and is simpler to implement. We provide open source Julia software on GitHub to replicate all calculations in this work, and we explain how to generalise our approach to deal with coupled PDE models with several dependent variables through a multinomial measurement error model, as well as pointing out other potential generalisations by linking our work with established practices in the field of generalised linear models.

A general approach for analysis of physiologically structured population models: the R package ‘PSPManalysis’

Steady-state analysis of structured population models

Physiologically based demographic model/GIS analyses of thirteen invasive species in Africa: why the biology matters

Long-term analysis of phenotypically structured models

Dynamic Analysis of a Predator and Prey Model with Some Computational Simulations

A framework of population dynamics from first principles

Ecological non-linear state space model selection via adaptive particle Markov chain Monte Carlo (AdPMCMC)

A simple model of population dynamics with beneficial and harmful interaction networks for empirical applications.

A tutorial on spatiotemporal partially observed Markov process models via the R package spatPomp

A general DDE framework to describe insect populations: Why delays are so important?

metaRange: A framework to build mechanistic range models

Multidimensional Stochastic Process Model and its Applications to Analysis of Longitudinal Data with Genetic Information

Modelling count data with partial differential equation models in biology

Applying allometric scaling to predator-prey systems

Advances in global sensitivity analyses of demographic-based species distribution models to address uncertainties in dynamic landscapes

ecode: An R package to investigate community dynamics in ordinary differential equation systems

EEMtoolbox: A user-friendly R package for flexible ensemble ecosystem modeling

Exploring the emergence and evolution of population patterns of leisure-time physical activity through agent-based modelling

A general modelling framework for open wildlife populations based on the polya tree prior

Mathematical models for age-structured population dynamics: An overview

Integrating eco‐evolutionary dynamics into matrix population models for structured populations: Discrete and continuous frameworks