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.

Partial differential equation about the prevalence of a chronic disease in the presence of duration dependency

On Keiding's Equation and its relation to differential equations about prevalence and incidence in chronic disease epidemiology

On characteristics of an ordinary differential equation and a related inverse problem in epidemiology

New ways of estimating excess mortality of chronic diseases: Insights from the illness-death model

Change rates and prevalence of a dichotomous variable: simulations and applications

A new stochastic differential equation modelling incidence and prevalence with an application to systemic lupus erythematosus in England and Wales, 1995

Stochastic differential equation for modelling health related quality of life

Simulation of Populations in a Time-, Age- and Duration Dependent Illness-Death Model

Estimation of excess mortality in a chronic condition from current status data with disease duration: simulation study about need for long-term care

Illness-death model with renewal

An identifiability problem in a state model for partly undetected chronic diseases

Application of the chemical master equation and its analytical solution to the illness-death model

State model for partly undetected non-communicable diseases (NCDs)

Modelling physiologically structured populations: renewal equations and partial differential equations

An Epidemic Model with Time Delay Determined by the Disease Duration

Time-Fractional Differential Equations with Variable Order using RDTM and ADM: Application to Infectious-Disease Model

Stochastic Partial Differential Equation SEIRS Epidemic Models: Well-posedness and Longtime Behavior

Maximum likelihood estimation for aggregate current status data: Simulation study using the illness-death model for chronic diseases with duration dependency

Modelling count data with partial differential equation models in biology

Delay differential equations for the spatially-resolved simulation of epidemics with specific application to COVID-19

A nonstandard numerical scheme for a novel SECIR integro-differential equation-based model allowing nonexponentially distributed stay times