Quantifying the Uncertainty in Model Parameters Using Gaussian Process-Based Markov Chain Monte Carlo: An Application to Cardiac Electrophysiological Models

06/02/2020
by   Jwala Dhamala, et al.
12

Estimation of patient-specific model parameters is important for personalized modeling, although sparse and noisy clinical data can introduce significant uncertainty in the estimated parameter values. This importance source of uncertainty, if left unquantified, will lead to unknown variability in model outputs that hinder their reliable adoptions. Probabilistic estimation model parameters, however, remains an unresolved challenge because standard Markov Chain Monte Carlo sampling requires repeated model simulations that are computationally infeasible. A common solution is to replace the simulation model with a computationally-efficient surrogate for a faster sampling. However, by sampling from an approximation of the exact posterior probability density function (pdf) of the parameters, the efficiency is gained at the expense of sampling accuracy. In this paper, we address this issue by integrating surrogate modeling into Metropolis Hasting (MH) sampling of the exact posterior pdfs to improve its acceptance rate. It is done by first quickly constructing a Gaussian process (GP) surrogate of the exact posterior pdfs using deterministic optimization. This efficient surrogate is then used to modify commonly-used proposal distributions in MH sampling such that only proposals accepted by the surrogate will be tested by the exact posterior pdf for acceptance/rejection, reducing unnecessary model simulations at unlikely candidates. Synthetic and real-data experiments using the presented method show a significant gain in computational efficiency without compromising the accuracy. In addition, insights into the non-identifiability and heterogeneity of tissue properties can be gained from the obtained posterior distributions.

READ FULL TEXT

page 8

page 10

research
10/13/2021

Fast Posterior Estimation of Cardiac Electrophysiological Model Parameters via Bayesian Active Learning

Probabilistic estimation of cardiac electrophysiological model parameter...
research
11/17/2019

Iterative Construction of Gaussian Process Surrogate Models for Bayesian Inference

A new algorithm is developed to tackle the issue of sampling non-Gaussia...
research
09/09/2019

Machine learning accelerates parameter optimization and uncertainty assessment of a land surface model

The performance of land surface models (LSMs) strongly depends on their ...
research
10/06/2022

Fast and robust parameter estimation with uncertainty quantification for the cardiac function

Parameter estimation and uncertainty quantification are crucial in compu...
research
11/23/2020

Probabilistic modeling of discrete structural response with application to composite plate penetration models

Discrete response of structures is often a key probabilistic quantity of...
research
07/10/2018

Surrogate-Based Inverse Modeling of the Hydrological System: An Adaptive Approach Considering Surrogate Structural Error

Inverse modeling is vital for an improved hydrological prediction. Howev...
research
07/10/2018

Surrogate-Based Bayesian Inverse Modeling of the Hydrological System: An Adaptive Approach Considering Surrogate Structural Error

Inverse modeling is vital for an improved hydrological prediction. Howev...

Please sign up or login with your details

Forgot password? Click here to reset