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.

Parameter estimation for discretely-observed linear birth-and-death processes

Parameter Identification for Fractional Ornstein–Uhlenbeck Processes Based on Discrete Observation

Parameter estimation for fractional Ornstein–Uhlenbeck processes at discrete observation

Parameter estimation for a hidden linear birth and death process with immigration

Likelihood-based inference for discretely observed birth-death-shift processes, with applications to evolution of mobile genetic elements

Estimation of the parameters of a stochastic logistic growth model

Parameter Estimation for Generalized Diffusion Processes with Reflected Boundary

Consistent least squares estimation in population-size-dependent branching processes

Inferring Density-Dependent Population Dynamics Mechanisms through Rate Disambiguation for Logistic Birth-Death Processes

Parameter estimation and long-range dependence of the fractional binomial process

Parameter estimation in the stochastic SIR model via scaled geometric Brownian motion

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

Targeted estimation of state occupation probabilities for the non‐Markov illness‐death model

Estimation in a semiparametric modulated renewal process

Fitting birth-death processes to panel data with applications to bacterial DNA fingerprinting

Estimation of the lifetime distribution from fluctuations in Bellman-Harris processes

Robust Non-Parametric Mortality and Fertility Modelling and Forecasting: Gaussian Process Regression Approaches

Parameter inference and model selection in deterministic and stochastic dynamical models via approximate Bayesian computation: modeling a wildlife epidemic

Regularized Bidimensional Estimation of the Hazard Rate

Modelling count data with partial differential equation models in biology

Accelerating Bayesian inference for stochastic epidemic models using incidence data