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.

Generalized sensitivity functions for size-structured population models

Generalized Finite Difference Method for a Class of Multidimensional Space-Fractional Diffusion Equations

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Derivative-Based Global Sensitivity Analysis for Models with High-Dimensional Inputs and Functional Outputs

Parameter Estimation for Dynamical Systems with Discrete Events and Logical Operations

On the structural sensitivity of some diffusion-reaction models of population dynamics

Deep Generalized Green's Functions

Modelling count data with partial differential equation models in biology

A discontinuous Galerkin method for systems of stochastic differential equations with applications to population biology, finance, and physics

Uncertainty quantification for generalized Langevin dynamics

On the existence of non-trivial steady-state size-distributions for a class of flocculation equations

Fast pick-freeze estimation of Sobol' sensitivity maps using basis expansions

Universally Sloppy Parameter Sensitivities in Systems Biology

Large population limit of interacting population dynamics via generalized gradient structures

Mixed-Type Discontinuous Galerkin Approach for Solving the Generalized FitzHugh–Nagumo Reaction–Diffusion Model

Parameter estimation and model selection for stochastic differential equations for biological growth

Analytical Extensions of the Generalized Moment Method for the Smoluchowski Coagulation Equation

Estimation of Parameter Distributions for Reaction-Diffusion Equations with Competition using Aggregate Spatiotemporal Data

Anomalous dimensions of the Smoluchowski coagulation equation

A New Approximation Approach for Analytically Solving the Population Balance Equation Due to Thermophoretic Coagulation

Computationally Scalable Bayesian SPDE Modeling for Censored Spatial Responses