TY - JOUR

T1 - Computational modelling with functional differential equations

T2 - Identification, selection, and sensitivity

AU - Baker, C. T.H.

AU - Bocharov, G. A.

AU - Paul, C. A.H.

AU - Rihan, F. A.

N1 - Funding Information:
* Corresponding author. E-mail addresses: cthbaker@ma.man.ac.uk, cthbaker@na-net.ornl.gov (C.T.H. Baker), bocharov@inm.ras.ru (G.A. Bocharov), chris@ma.man.ac.uk, cpaul@na-net.ornl.gov (C.A.H. Paul), f.rihan@salford.ac.uk, frihan@na-net.ornl.gov (F.A. Rihan). 1 Research Professor, Manchester University; Visiting Professor, University College Chester. 2 INM, Russian Academy of Sciences, Moscow; Honorary Research Fellow, Manchester; Leverhulme Visiting Professor, University College Chester. Supported by the Leverhulme Trust and the Russian Foundation for Basic Research. 3 Research Fellow. 4 University of Salford; Honorary Research Fellow, Manchester. Supported in part by NERC and the Egyptian Government.

PY - 2005/5

Y1 - 2005/5

N2 - Mathematical models based upon certain types of differential equations, functional differential equations, or systems of such equations, are often employed to represent the dynamics of natural, in particular biological, phenomena. We present some of the principles underlying the choice of a methodology (based on observational data) for the computational identification of, and discrimination between, quantitatively consistent models, using scientifically meaningful parameters. We propose that a computational approach is essential for obtaining meaningful models. For example, it permits the choice of realistic models incorporating a time-lag which is entirely natural from the scientific perspective. The time-lag is a feature that can permit a close reconciliation between models incorporating computed parameter values and observations. Exploiting the link between information theory, maximum likelihood, and weighted least squares, and with distributional assumptions on the data errors, we may construct an appropriate objective function to be minimized computationally. The minimizer is sought over a set of parameters (which may include the time-lag) that define the model. Each evaluation of the objective function requires the computational solution of the parametrized equations defining the model. To select a parametrized model, from amongst a family or hierarchy of possible best-fit models, we are able to employ certain indicators based on information-theoretic criteria. We can evaluate confidence intervals for the parameters, and a sensitivity analysis provides an expression for an information matrix, and feedback on the covariances of the parameters in relation to the best fit. This gives a firm basis for any simplification of the model (e.g., by omitting a parameter).

AB - Mathematical models based upon certain types of differential equations, functional differential equations, or systems of such equations, are often employed to represent the dynamics of natural, in particular biological, phenomena. We present some of the principles underlying the choice of a methodology (based on observational data) for the computational identification of, and discrimination between, quantitatively consistent models, using scientifically meaningful parameters. We propose that a computational approach is essential for obtaining meaningful models. For example, it permits the choice of realistic models incorporating a time-lag which is entirely natural from the scientific perspective. The time-lag is a feature that can permit a close reconciliation between models incorporating computed parameter values and observations. Exploiting the link between information theory, maximum likelihood, and weighted least squares, and with distributional assumptions on the data errors, we may construct an appropriate objective function to be minimized computationally. The minimizer is sought over a set of parameters (which may include the time-lag) that define the model. Each evaluation of the objective function requires the computational solution of the parametrized equations defining the model. To select a parametrized model, from amongst a family or hierarchy of possible best-fit models, we are able to employ certain indicators based on information-theoretic criteria. We can evaluate confidence intervals for the parameters, and a sensitivity analysis provides an expression for an information matrix, and feedback on the covariances of the parameters in relation to the best fit. This gives a firm basis for any simplification of the model (e.g., by omitting a parameter).

KW - Computation

KW - Data

KW - Differential equations

KW - Identifiability

KW - Information-theoretic criteria

KW - Modelling

KW - Objective function

KW - Parametric estimation

KW - Sensitivity

KW - Time-lag

KW - Well-posedness

UR - http://www.scopus.com/inward/record.url?scp=14844354032&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=14844354032&partnerID=8YFLogxK

U2 - 10.1016/j.apnum.2004.08.014

DO - 10.1016/j.apnum.2004.08.014

M3 - Article

AN - SCOPUS:14844354032

SN - 0168-9274

VL - 53

SP - 107

EP - 129

JO - Applied Numerical Mathematics

JF - Applied Numerical Mathematics

IS - 2-4

ER -