Exploiting the bootstrap method for quantifying parameter confidence intervals in dynamical systems.

A quantitative description of dynamical systems requires the estimation of uncertain kinetic parameters and an analysis of their precision. A method frequently used to describe the confidence intervals of estimated parameters is based on the Fisher-Information-Matrix. The application of this traditional method has two important shortcomings: (i) it gives only lower bounds for the variance of a parameter if the solution of the underlying model equations is non-linear in parameters. (ii) The resulting confidence interval is symmetric with respect to the estimated parameter. Here, we show that by applying the bootstrap method a better approximation of (possibly) asymmetric confidence intervals for parameters could be obtained. In contrast to previous applications devoted to non-parametric problems, a dynamical model describing a bio-chemical network is used to evaluate the method.

[1]  Steven P. Asprey,et al.  On the design of optimally informative dynamic experiments for model discrimination in multiresponse nonlinear situations , 2003 .

[2]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[3]  D Weuster-Botz,et al.  Automated sampling device for monitoring intracellular metabolite dynamics. , 1999, Analytical biochemistry.

[4]  D. Freedman,et al.  On the histogram as a density estimator:L2 theory , 1981 .

[5]  B. Efron,et al.  Bootstrap confidence intervals , 1996 .

[6]  Sandro Macchietto,et al.  Statistical tools for optimal dynamic model building , 2000 .

[7]  S. T. Buckland,et al.  An Introduction to the Bootstrap. , 1994 .

[8]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[9]  E. Gilles,et al.  The organization of metabolic reaction networks. III. Application for diauxic growth on glucose and lactose. , 2001, Metabolic engineering.

[10]  J. Glasner,et al.  Genome-wide expression profiling in Escherichia coli K-12. , 1999, Nucleic acids research.

[11]  Matthias Reuss,et al.  Optimal Experimental Design for Parameter Estimation in Unstructured Growth Models , 1994 .

[12]  Lennart Ljung,et al.  System Identification: Theory for the User , 1987 .

[13]  Daniel E. Zak,et al.  Importance of input perturbations and stochastic gene expression in the reverse engineering of genetic regulatory networks: insights from an identifiability analysis of an in silico network. , 2003, Genome research.

[14]  Ursula Klingmüller,et al.  Simulation Methods for Optimal Experimental Design in Systems Biology , 2003, Simul..