Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Skip to main content

A Homotopy Method for Parameter Estimation of Nonlinear Differential Equations with Multiple Optima

  • Published:
Journal of Scientific Computing Aims and scope Submit manuscript

Abstract

A numerical method for estimating multiple parameter values of nonlinear systems arising from biology is presented. The uncertain parameters are modeled as random variables. Then the solutions are expressed as convergent series of orthogonal polynomial expansions in terms of the input random parameters. Homotopy continuation method is employed to solve the resulting polynomial system, and more importantly, to compute the multiple optimal parameter values. Several numerical examples, from a single equation to problems with relatively complicated forms of governing equations, are used to demonstrate the robustness and effectiveness of this numerical method.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4

Similar content being viewed by others

References

  1. Bates, D., Hauenstein, J., Sommese, A., Wampler, C.: Numerically Solving Polynomial Systems with Bertini, vol. 25. SIAM, Philadelphia (2013)

    MATH  Google Scholar 

  2. Blanchard, E.D., Sandu, A., Sandu, C.: A polynomial chaos-based Kalman filter approach for parameter estimation of mechanical systems. J. Dyn. Syst. Meas. Control 132(6), 061404 (2010)

    Article  Google Scholar 

  3. Coelho, F.C., Codeco, C.T., Gomes, M.G.: A Bayesian framework for parameter estimation in dynamical models. PLoS ONE 6(5), e19616 (2011)

    Article  Google Scholar 

  4. Dobson, I.: Computing a closest bifurcation instability in multidimensional parameter space. J. Nonlinear Sci. 3(1), 307–327 (1993)

    Article  MathSciNet  MATH  Google Scholar 

  5. Friedman, A., Hao, W.: Mathematical modeling of liver fibrosis. Math. Biosci. Eng. 14(1), 143–164 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  6. Gallas, J.: Structure of the parameter space of the Hénon map. Phys. Rev. Lett. 70(18), 2714 (1993)

    Article  Google Scholar 

  7. Gardner, T.S., Cantor, C.R., Collins, J.J.: Construction of a genetic toggle switch in Escherichia coli. Nature 403(6767), 339–342 (2000)

    Article  Google Scholar 

  8. Gomes, W., Beck, A., da Silva Jr, C.: Modeling random corrosion processes via polynomial chaos expansion. J. Braz. Soc. Mech. Sci. Eng. 34(SPE2), 561–568 (2012)

    Google Scholar 

  9. Hao, W., Crouser, E.D., Friedman, A.: Mathematical model of sarcoidosis. Proc. Natl. Acad. Sci. U.S.A. 111(45), 16065–16070 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  10. Hao, W., Friedman, A.: Mathematical model on Alzheimer’s disease. BMC Syst. Biol. 10(1), 108 (2016)

    Article  Google Scholar 

  11. Hao, W., Hauenstein, J., Hu, B., Liu, Y., Sommese, A., Zhang, Y.-T.: Bifurcation for a free boundary problem modeling the growth of a tumor with a necrotic core. Nonlinear Anal. Real World Appl. 13(2), 694–709 (2012)

    Article  MathSciNet  MATH  Google Scholar 

  12. Hao, W., Hauenstein, J., Hu, B., Sommese, A.: A three-dimensional steady-state tumor system. Appl. Math. Comput. 218(6), 2661–2669 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  13. Hao, W., Hauenstein, J., Hu, B., Sommese, A.: A bootstrapping approach for computing multiple solutions of differential equations. J. Comput. Appl. Math. 258, 181–190 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  14. Hao, W., Hauenstein, J., Shu, C.-W., Sommese, A., Xu, Z., Zhang, Y.-T.: A homotopy method based on weno schemes for solving steady state problems of hyperbolic conservation laws. J. Comput. Phys. 250, 332–346 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  15. Hao, W., Marsh, C., Friedman, A.: A mathematical model of idiopathic pulmonary fibrosis. PLoS ONE 10(9), e0135097 (2015)

    Article  Google Scholar 

  16. Hao, W., Nepomechie, R., Sommese, A.: Completeness of solutions of bethe’s equations. Phys. Rev. E 88(5), 052113 (2013)

    Article  Google Scholar 

  17. Hao, W., Nepomechie, R., Sommese, A.: Singular solutions, repeated roots and completeness for higher-spin chains. J. Stat. Mech. Theory Exp. 2014(3), P03024 (2014)

    Article  MathSciNet  Google Scholar 

  18. Hao, W., Schlesinger, L.S., Friedman, A.: Modeling granulomas in response to infection in the lung. PLoS ONE 11(3), e0148738 (2016)

    Article  Google Scholar 

  19. Jaakkola, T.S., Jordan, M.I.: Bayesian parameter estimation via variational methods. Stat. Comput. 10(1), 25–37 (2000)

    Article  Google Scholar 

  20. Kramer, S.C., Sorenson, H.W.: Bayesian parameter estimation. IEEE Trans. Autom. Control 33(2), 217–222 (1988)

    Article  MATH  Google Scholar 

  21. Lal, A., Lash, A.E., Altschul, S.F., Velculescu, V., Zhang, L., McLendon, R.E., Marra, M.A., Prange, C., Morin, P.J., Polyak, K., Papadopoulos, N., Vogelstein, B., Kinzler, K.W., Strausberg, R.L., Riggins, G.J.: A public database for gene expression in human cancers. Cancer Res. 59(21), 5403–5407 (1999)

    Google Scholar 

  22. Li, J., Xiu, D.: A generalized polynomial chaos based ensemble Kalman filter with high accuracy. J. Comput. phys. 228(15), 5454–5469 (2009)

    Article  MathSciNet  MATH  Google Scholar 

  23. Lillacci, G., Khammash, M.: Parameter estimation and model selection in computational biology. PLoS Comput. Biol. 6(3), e1000696 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  24. Pence, B.L., Fathy, H.K., Stein, J.L.: A maximum likelihood approach to recursive polynomial chaos parameter estimation. In: American Control Conference (ACC), pp. 2144–2151. IEEE (2010)

  25. Sommese, A., Verschelde, J., Wampler, C.: Numerical algebraic geometry. In: The Mathematics of Numerical Analysis, vol. 32 of Lectures in Applied Mathematics. Citeseer (1996)

  26. Sommese, A., Wampler, C.: The Numerical Solution of Systems of Polynomials Arising in Engineering and Science, vol. 99. World Scientific, Singapore (2005)

    Book  MATH  Google Scholar 

  27. Sturrock, M., Hao, W., Schwartzbaum, J., Rempala, G.A.: A mathematical model of pre-diagnostic glioma growth. J. Theor. Biol. 380, 299–308 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  28. Tarantola, A.: Inverse Problem Theory and Methods for Model Parameter Estimation. SIAM, Philadelphia (2005)

    Book  MATH  Google Scholar 

  29. Vanlier, J., Tiemann, C., Hilbers, P., van Riel, N.: Parameter uncertainty in biochemical models described by ordinary differential equations. Math. Biosci. 246(2), 305–314 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  30. Voit, E., Chou, I.: Parameter estimation in canonical biological systems models. Int. J. Syst. Synth. Biol. 1, 1–19 (2010)

    Google Scholar 

  31. Wu, H., Wang, F., Chang, M.: Dynamic sensitivity analysis of biological systems. BMC Bioinf. 9(12), S17 (2008)

    Article  Google Scholar 

  32. Xiu, D., Karniadakis, G.E.: Modeling uncertainty in steady state diffusion problems via generalized polynomial chaos. Comput. Methods Appl. Mech. Eng. 191(43), 4927–4948 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  33. Xiu, D., Karniadakis, G.E.: The Wiener–Askey polynomial chaos for stochastic differential equations. SIAM J. Sci. Comput. 24(2), 619–644 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  34. Xiu, D., Karniadakis, G.E.: Modeling uncertainty in flow simulations via generalized polynomial chaos. J. Comput. Phys. 187(1), 137–167 (2003)

    Article  MathSciNet  MATH  Google Scholar 

  35. Zhan, C., Yeung, L.: Parameter estimation in systems biology models using spline approximation. BMC Syst. Biol. 5(1), 14 (2011)

    Article  Google Scholar 

  36. Zhusubaliyev, Z., Mosekilde, E.: Bifurcations and Chaos in Piecewise-Smooth Dynamical Systems: Applications to Power Converters, Relay and Pulse-Width Modulated Control Systems, and Human Decision-Making Behavior, vol. 44. World Scientific, Singapore (2003)

    Google Scholar 

Download references

Acknowledgements

The authors have been supported by the American Heart Association under Grant 17SDG33660722.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Wenrui Hao.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hao, W. A Homotopy Method for Parameter Estimation of Nonlinear Differential Equations with Multiple Optima. J Sci Comput 74, 1314–1324 (2018). https://doi.org/10.1007/s10915-017-0518-4

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10915-017-0518-4

Keywords