Page 191 - Artificial Intelligence for Computational Modeling of the Heart
P. 191

Chapter 5 Machine learning methods for robust parameter estimation  163




                     a reinforcement learning method to train an agent that estimates
                     model parameters using a learned strategy. Both methods reach
                     either similar or superior performance compared to highly en-
                     gineered inverse problem techniques, while being more efficient
                     and precise.


                     5.2 A regression approach to model

                          parameter estimation
                        A first data driven approach consists in learning a regression
                     model to estimate parameters of a physiological model directly
                     from clinical measurements, as previously described in [220].
                     Without loss of generality, the method is presented in the context
                     of cardiac electrophysiology (EP), but it could be applied to mod-
                     els describing other aspects of cardiac function or other organs.
                     Furthermore, we will also illustrate how uncertainty estimates can
                     be calculated by leveraging the large database used to train the
                     regression model.


                     5.2.1 Data-driven estimation of myocardial electrical
                           diffusivity
                        The monodomain model introduced in section 1.2.2 is used
                     without loss of generality as any other EP model could be also
                     employed. The EP model y = f(x) takes as input the diffusivity
                     coefficients x = (c Myo ,c LV ,c RV ), among other parameters that are
                     fixed in this example, and returns the ECG features y = (QRS ,α),
                                                                           d
                     where QRSd is the QRS complex duration measured on 12-lead
                     ECG traces, and α is the electrical axis of the QRS complex. The
                     goal is to learn a function g that approximates the inverse problem
                     x = g(z) ≈ f −1 (z), namely that takes as input measured 12-lead
                     ECG parameters (z), and returns an estimate of the model param-
                     eters x.
                        In this setting, though, ground truth data, namely pairs of clin-
                     ical measurements (z) and diffusivity values (x), is not available.
                     We therefore leverage the computational model to create a large
                     database of pairs (y,x), assuming the model is accurate enough
                     to provide realistic ECG features. The database is created by uni-
                     formly sampling the space of parameters x, feeding the model
                     with the sampled parameters and collecting the respective out-
                     puts y. Once the database is ready, one could directly learn a
                     regression model x = g(y). Unfortunately this task can be made
                     difficult by the lack of observability of the parameters. Two sets
                     of parameters x 1 and x 2 could give very similar observed ECG
   186   187   188   189   190   191   192   193   194   195   196