Data-based identifiability analysis of non-linear dynamical models

• Published in: Bioinformatics 2007-10, Vol.23 (19), p.2612-2618
• Main Authors: Hengl, S., Kreutz, C., Timmer, J., Maiwald, T.
• Format: Article
• Language: English
• Subjects: Algorithms; Artificial Intelligence; Biological and medical sciences; Computer Simulation; Databases, Factual; Fundamental and applied biological sciences. Psychology; General aspects; Information Storage and Retrieval - methods; Mathematics in biology. Statistical analysis. Models. Metrology. Data processing in biology (general aspects); Models, Biological; Nonlinear Dynamics; Pattern Recognition, Automated - methods; Proteome - metabolism; Signal Transduction - physiology
• ISSN: 1367-4803; 1460-2059; 1367-4811
• DOI: 10.1093/bioinformatics/btm382

Motivation: Mathematical modelling of biological systems is becoming a standard approach to investigate complex dynamic, non-linear interaction mechanisms in cellular processes. However, models may comprise non-identifiable parameters which cannot be unambiguously determined. Non-identifiability manifests itself in functionally related parameters, which are difficult to detect. Results: We present the method of mean optimal transformations, a non-parametric bootstrap-based algorithm for identifiability testing, capable of identifying linear and non-linear relations of arbitrarily many parameters, regardless of model size or complexity. This is performed with use of optimal transformations, estimated using the alternating conditional expectation algorithm (ACE). An initial guess or prior knowledge concerning the underlying relation of the parameters is not required. Independent, and hence identifiable parameters are determined as well. The quality of data at disposal is included in our approach, i.e. the non-linear model is fitted to data and estimated parameter values are investigated with respect to functional relations. We exemplify our approach on a realistic dynamical model and demonstrate that the variability of estimated parameter values decreases from 81 to 1% after detection and fixation of structural non-identifiabilities. Availability: Our algorithm is written in Matlab and R. It is available from the authors on request. An implementation of ACE, written in Matlab as well as in C, is available online at www.stefanhengl.de Contact: [e-mail] Supplementary information: Supplementary data are available at Bioinformatics online.