Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations
Biological system's dynamics are increasingly studied with nonlinear ordinary differential equations, whose parameters are estimated from input/output experimental data. Structural identifiability analysis addresses the theoretical question whether the inverse problem of recovering the unknown...
Main Authors: | , |
---|---|
Format: | Article |
Language: | English |
Published: |
AIMS Press
2019-07-01
|
Series: | Mathematical Biosciences and Engineering |
Subjects: | |
Online Access: | https://www.aimspress.com/article/10.3934/mbe.2019322?viewType=HTML |
_version_ | 1811234571330519040 |
---|---|
author | Maria Pia Saccomani Karl Thomaseth |
author_facet | Maria Pia Saccomani Karl Thomaseth |
author_sort | Maria Pia Saccomani |
collection | DOAJ |
description | Biological system's dynamics are increasingly studied with nonlinear ordinary differential equations, whose parameters are estimated from input/output experimental data. Structural identifiability analysis addresses the theoretical question whether the inverse problem of recovering the unknown parameters from noise-free data is uniquely solvable (global), or if there is a finite (local), or an infinite number (non identifiable) of parameter values that generate identical input/output trajectories. In contrast, practical identifiability analysis aims to assess whether the experimental data provide information on the parameter estimates in terms of precision and accuracy. A main difference between the two identifiability approaches is that the former is mostly carried out analytically and provides exact results at a cost of increased computational complexity, while the latter is usually numerically tested by calculating statistical confidence regions and relies on decision thresholds. Here we focus on local identifiability, a critical issue in biological modeling. This is the case when a model has multiple parameter solutions which equivalently describe the input/output data, but predict different behaviours of the unmeasured variables, often those of major interest. We present theoretical background and applications to locally identifiable ODE models described by rational functions. We show how structural identifiability analysis completes the practical identifiability results. In particular we propose an algorithmic approach, implemented with our software DAISY, to calculate all numerical parameter solutions and to predict the corresponding behaviour of the unmeasured variables, which otherwise would remain hidden. A case study of a locally identifiable HIV model shows that one should be aware of the presence of multiple parameter solutions to comprehensively describe the biological system and avoid biological misinterpretation of the results. |
first_indexed | 2024-04-12T11:39:43Z |
format | Article |
id | doaj.art-5e0bb63bee3647d3903b3f5e1d87ad87 |
institution | Directory Open Access Journal |
issn | 1551-0018 |
language | English |
last_indexed | 2024-04-12T11:39:43Z |
publishDate | 2019-07-01 |
publisher | AIMS Press |
record_format | Article |
series | Mathematical Biosciences and Engineering |
spelling | doaj.art-5e0bb63bee3647d3903b3f5e1d87ad872022-12-22T03:34:43ZengAIMS PressMathematical Biosciences and Engineering1551-00182019-07-011666438645310.3934/mbe.2019322Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretationsMaria Pia Saccomani0Karl Thomaseth11. Department of Information Engineering, University of Padova, Padova, 35131 Italy2. National Research Council (IEIIT-CNR) c/o Department of Information Engineering, University of Padova, Padova, 35131 ItalyBiological system's dynamics are increasingly studied with nonlinear ordinary differential equations, whose parameters are estimated from input/output experimental data. Structural identifiability analysis addresses the theoretical question whether the inverse problem of recovering the unknown parameters from noise-free data is uniquely solvable (global), or if there is a finite (local), or an infinite number (non identifiable) of parameter values that generate identical input/output trajectories. In contrast, practical identifiability analysis aims to assess whether the experimental data provide information on the parameter estimates in terms of precision and accuracy. A main difference between the two identifiability approaches is that the former is mostly carried out analytically and provides exact results at a cost of increased computational complexity, while the latter is usually numerically tested by calculating statistical confidence regions and relies on decision thresholds. Here we focus on local identifiability, a critical issue in biological modeling. This is the case when a model has multiple parameter solutions which equivalently describe the input/output data, but predict different behaviours of the unmeasured variables, often those of major interest. We present theoretical background and applications to locally identifiable ODE models described by rational functions. We show how structural identifiability analysis completes the practical identifiability results. In particular we propose an algorithmic approach, implemented with our software DAISY, to calculate all numerical parameter solutions and to predict the corresponding behaviour of the unmeasured variables, which otherwise would remain hidden. A case study of a locally identifiable HIV model shows that one should be aware of the presence of multiple parameter solutions to comprehensively describe the biological system and avoid biological misinterpretation of the results.https://www.aimspress.com/article/10.3934/mbe.2019322?viewType=HTMLbiological systemslocal identifiabilitymultiple parameter solutionsdifferential equationparameter estimationhiv model |
spellingShingle | Maria Pia Saccomani Karl Thomaseth Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations Mathematical Biosciences and Engineering biological systems local identifiability multiple parameter solutions differential equation parameter estimation hiv model |
title | Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations |
title_full | Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations |
title_fullStr | Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations |
title_full_unstemmed | Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations |
title_short | Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations |
title_sort | calculating all multiple parameter solutions of ode models to avoid biological misinterpretations |
topic | biological systems local identifiability multiple parameter solutions differential equation parameter estimation hiv model |
url | https://www.aimspress.com/article/10.3934/mbe.2019322?viewType=HTML |
work_keys_str_mv | AT mariapiasaccomani calculatingallmultipleparametersolutionsofodemodelstoavoidbiologicalmisinterpretations AT karlthomaseth calculatingallmultipleparametersolutionsofodemodelstoavoidbiologicalmisinterpretations |