Gaussian process enhanced semi-automatic approximate Bayesian computation: parameter inference in a stochastic differential equation system for chemotaxis

Borowska, A. , Giurghita, D. and Husmeier, D. (2021) Gaussian process enhanced semi-automatic approximate Bayesian computation: parameter inference in a stochastic differential equation system for chemotaxis. Journal of Computational Physics, 429, 109999. (doi: 10.1016/j.jcp.2020.109999)

[img] Text
226089.pdf - Accepted Version
Available under License Creative Commons Attribution Non-commercial No Derivatives.

15MB

Abstract

Chemotaxis is a type of cell movement in response to a chemical stimulus which plays a key role in multiple biophysical processes, such as embryogenesis and wound healing, and which is crucial for understanding metastasis in cancer research. In the literature, chemotaxis has been modelled using biophysical models based on systems of nonlinear stochastic partial differential equations (NSPDEs), which are known to be challenging for statistical inference due to the intractability of the associated likelihood and the high computational costs of their numerical integration. Therefore, data analysis in this context has been limited to comparing predictions from NSPDE models to laboratory data using simple descriptive statistics. We present a statistically rigorous framework for parameter estimation in complex biophysical systems described by NSPDEs such as the one of chemotaxis. We adopt a likelihood-free approach based on approximate Bayesian computations with sequential Monte Carlo (ABC-SMC) which allows for circumventing the intractability of the likelihood. To find informative summary statistics, crucial for the performance of ABC, we propose to use a Gaussian process (GP) regression model. The interpolation provided by the GP regression turns out useful on its own merits: it relatively accurately estimates the parameters of the NSPDE model and allows for uncertainty quantification, at a very low computational cost. Our proposed methodology allows for a considerable part of computations to be completed before having observed any data, providing a practical toolbox to experimental scientists whose modes of operation frequently involve experiments and inference taking place at distinct points in time. In an application to externally provided synthetic data we demonstrate that the correction provided by ABC-SMC is essential for accurate estimation of some of the NSPDE model parameters and for more flexible uncertainty quantification.

Item Type:Articles
Status:Published
Refereed:Yes
Glasgow Author(s) Enlighten ID:Husmeier, Professor Dirk and Borowska, Dr Agnieszka and Giurghita, Miss Diana
Authors: Borowska, A., Giurghita, D., and Husmeier, D.
College/School:College of Science and Engineering > School of Mathematics and Statistics > Statistics
Journal Name:Journal of Computational Physics
Publisher:Elsevier
ISSN:0021-9991
ISSN (Online):1090-2716
Published Online:13 November 2020
Copyright Holders:Copyright © 2020 Elsevier Inc.
First Published:First published in Journal of Computational Physics 429: 109999
Publisher Policy:Reproduced in accordance with the publisher copyright policy

University Staff: Request a correction | Enlighten Editors: Update this record

Project CodeAward NoProject NamePrincipal InvestigatorFunder's NameFunder RefLead Dept
303232EPSRC Centre for Multiscale soft tissue mechanics with MIT and POLIMI (SofTMech-MP)Xiaoyu LuoEngineering and Physical Sciences Research Council (EPSRC)EP/S030875/1M&S - Mathematics
300982Exploiting Closed-Loop Aspects in Computationally and Data Intensive AnalyticsRoderick Murray-SmithEngineering and Physical Sciences Research Council (EPSRC)EP/R018634/1Computing Science
172141EPSRC Centre for Multiscale soft tissue mechanics with application to heart & cancerRaymond OgdenEngineering and Physical Sciences Research Council (EPSRC)EP/N014642/1M&S - Mathematics
305655Inference of cardio-mechanical parameters in real time: moving mathematical modelling into the clinicDirk HusmeierThe Royal Society of Edinburgh (ROYSOCED)62335M&S - Statistics