Utility of a novel error-stepping method to improve gradient-based parameter identification by increasing the smoothness of the local objective surface: A case-study of pulmonary mechanics

Accurate model parameter identification relies on accurate forward model simulations to guide convergence. However, some forward simulation methodologies lack the precision required to properly define the local objective surface and can cause failed parameter identification. The role of objective surface smoothness in identification of a pulmonary mechanics model was assessed using forward simulation from a novel error-stepping method and a proprietary Runge-Kutta method. The objective surfaces were compared via the identified parameter discrepancy generated in a Monte Carlo simulation and the local smoothness of the objective surfaces they generate. The error-stepping method generated significantly smoother error surfaces in each of the cases tested (p<0.0001) and more accurate model parameter estimates than the Runge-Kutta method in three of the four cases tested (p<0.0001) despite a 75% reduction in computational cost. Of note, parameter discrepancy in most cases was limited to a particular oblique plane, indicating a non-intuitive multi-parameter trade-off was occurring. The error-stepping method consistently improved or equalled the outcomes of the Runge-Kutta time-integration method for forward simulations of the pulmonary mechanics model. This study indicates that accurate parameter identification relies on accurate definition of the local objective function, and that parameter trade-off can occur on oblique planes resulting prematurely halted parameter convergence.

[1]  Kenneth Levenberg A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[2]  Thomas Desaive,et al.  A graphical method for practical and informative identifiability analyses of physiological models: A case study of insulin kinetics and sensitivity , 2011, Biomedical engineering online.

[3]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[4]  Christopher E. Hann,et al.  A minimal model of lung mechanics and model-based markers for optimizing ventilator treatment in ARDS patients , 2009, Comput. Methods Programs Biomed..

[5]  Yeong Shiong Chiew,et al.  A Time-Continuous Model of Respiratory Mechanics of ARDS Patients , 2013 .

[6]  K. Hickling,et al.  The pressure-volume curve is greatly modified by recruitment. A mathematical model of ARDS lungs. , 1998, American journal of respiratory and critical care medicine.

[7]  Thomas Desaive,et al.  Physiological modeling, tight glycemic control, and the ICU clinician: what are models and how can they affect practice? , 2011, Annals of intensive care.

[8]  Ursula Klingmüller,et al.  Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood , 2009, Bioinform..

[9]  S Schumann,et al.  AUTOPILOT-BT: a system for knowledge and model based mechanical ventilation. , 2008, Technology and health care : official journal of the European Society for Engineering and Medicine.

[10]  J. Geoffrey Chase,et al.  Structural Identifiability and Practical Applicability of an Alveolar Recruitment Model for ARDS Patients , 2012, IEEE Transactions on Biomedical Engineering.

[11]  Maria Pia Saccomani,et al.  DAISY: A new software tool to test global identifiability of biological and physiological systems , 2007, Comput. Methods Programs Biomed..

[12]  J. Geoffrey Chase,et al.  Characterisation of the iterative integral parameter identification method , 2011, Medical & Biological Engineering & Computing.