Gaussian process emulation for discontinuous response surfaces with applications for cardiac electrophysiology models

05/25/2018
by   Sanmitra Ghosh, et al.
0

Mathematical models of biological systems are beginning to be used for safety-critical applications, where large numbers of repeated model evaluations are required to perform uncertainty quantification and sensitivity analysis. Most of these models are nonlinear both in variables and parameters/inputs which has two consequences. First, analytic solutions are rarely available so repeated evaluation of these models by numerically solving differential equations incurs a significant computational burden. Second, many models undergo bifurcations in behaviour as parameters are varied. As a result, simulation outputs often contain discontinuities as we change parameter values and move through parameter/input space. Statistical emulators such as Gaussian processes are frequently used to reduce the computational cost of uncertainty quantification, but discontinuities render a standard Gaussian process emulation approach unsuitable as these emulators assume a smooth and continuous response to changes in parameter values. In this article, we propose a novel two-step method for building a Gaussian Process emulator for models with discontinuous response surfaces. We first use a Gaussian Process classifier to detect boundaries of discontinuities and then constrain the Gaussian Process emulation of the response surface within these boundaries. We introduce a novel `certainty metric' to guide active learning for a multi-class probabilistic classifier. We apply the new classifier to simulations of drug action on a cardiac electrophysiology model, to propagate our uncertainty in a drug's action through to predictions of changes to the cardiac action potential. The proposed two-step active learning method significantly reduces the computational cost of emulating models that undergo multiple bifurcations.

READ FULL TEXT

page 20

page 22

page 23

research
01/09/2018

Known Boundary Emulation of Complex Computer Models

Computer models are now widely used across a range of scientific discipl...
research
05/28/2018

Dirichlet-based Gaussian Processes for Large-scale Calibrated Classification

In this paper, we study the problem of deriving fast and accurate classi...
research
01/13/2020

Considering discrepancy when calibrating a mechanistic electrophysiology model

Uncertainty quantification (UQ) is a vital step in using mathematical mo...
research
02/15/2021

Nonintrusive Uncertainty Quantification for automotive crash problems with VPS/Pamcrash

Uncertainty Quantification (UQ) is a key discipline for computational mo...
research
12/15/2020

Active Learning for Deep Gaussian Process Surrogates

Deep Gaussian processes (DGPs) are increasingly popular as predictive mo...
research
05/23/2019

A Bulirsch-Stoer algorithm using Gaussian processes

In this paper, we treat the problem of evaluating the asymptotic error i...
research
05/03/2021

Uncertainty quantification and estimation in differential dynamic microscopy

Differential dynamic microscopy (DDM) is a form of video image analysis ...

Please sign up or login with your details

Forgot password? Click here to reset