GpABC: a Julia package for approximate Bayesian computation with Gaussian process emulation

Bioinformatics. 2020 May 1;36(10):3286-3287. doi: 10.1093/bioinformatics/btaa078.

Abstract

Motivation: Approximate Bayesian computation (ABC) is an important framework within which to infer the structure and parameters of a systems biology model. It is especially suitable for biological systems with stochastic and nonlinear dynamics, for which the likelihood functions are intractable. However, the associated computational cost often limits ABC to models that are relatively quick to simulate in practice.

Results: We here present a Julia package, GpABC, that implements parameter inference and model selection for deterministic or stochastic models using (i) standard rejection ABC or sequential Monte Carlo ABC or (ii) ABC with Gaussian process emulation. The latter significantly reduces the computational cost.

Availability and implementation: https://github.com/tanhevg/GpABC.jl.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Bayes Theorem
  • Computer Simulation
  • Likelihood Functions
  • Monte Carlo Method
  • Normal Distribution
  • Systems Biology*