Markov chain Monte Carlo simulation using the DREAM software package: Theory, concepts, and MATLAB implementation. Bayesian inference has found widespread application and use in science and engineering to reconcile Earth system models with data, including prediction in space (interpolation), prediction in time (forecasting), assimilation of observations and deterministic/stochastic model output, and inference of the model parameters. Bayes theorem states that the posterior probability, p ( H | Y ¿ ) of a hypothesis, H is proportional to the product of the prior probability, p(H) of this hypothesis and the likelihood, L ( H | Y ¿ ) of the same hypothesis given the new observations, Y ¿ , or p ( H | Y ¿ ) ¿ p ( H ) L ( H | Y ¿ ) . In science and engineering, H often constitutes some numerical model, ¿(x) which summarizes, in algebraic and differential equations, state variables and fluxes, all knowledge of the system of interest, and the unknown parameter values, x are subject to inference using the data Y ¿ . Unfortunately, for complex system models the posterior distribution is often high dimensional and analytically intractable, and sampling methods are required to approximate the target. In this paper I review the basic theory of Markov chain Monte Carlo (MCMC) simulation and introduce a MATLAB toolbox of the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm developed by Vrugt et¿al. (2008a, 2009a) and used for Bayesian inference in fields ranging from physics, chemistry and engineering, to ecology, hydrology, and geophysics. This MATLAB toolbox provides scientists and engineers with an arsenal of options and utilities to solve posterior sampling problems involving (among others) bimodality, high-dimensionality, summary statistics, bounded parameter spaces, dynamic simulation models, formal/informal likelihood functions (GLUE), diagnostic model evaluation, data assimilation, Bayesian model averaging, distributed computation, and informative/noninformative prior distributions. The DREAM toolbox supports parallel computing and includes tools for convergence analysis of the sampled chain trajectories and post-processing of the results. Seven different case studies illustrate the main capabilities and functionalities of the MATLAB toolbox. A MATLAB toolbox of the DREAM algorithm is presented.This toolbox provides users with a great arsenal of options for Bayesian inference.The toolbox supports/implements parallel computing.Seven different case studies illustrate the main capabilities and functionalities of the toolbox.
Keywords for this software
References in zbMATH (referenced in 9 articles )
Showing results 1 to 9 of 9.
- Drovandi, Christopher; Nott, David J.; Pagendam, Daniel E.: A semiautomatic method for history matching using sequential Monte Carlo (2021)
- Padmanabha, Govinda Anantha; Zabaras, Nicholas: Solving inverse problems using conditional invertible neural networks (2021)
- Saibaba, Arvind K.; Prasad, Pranjal; de Sturler, Eric; Miller, Eric; Kilmer, Misha E.: Randomized approaches to accelerate MCMC algorithms for Bayesian inverse problems (2021)
- Jokhadze, Valeriane; Schmidt, Wolfgang M.: Measuring model risk in financial risk management and pricing (2020)
- Li, Jing; Tartakovsky, Alexandre M.: Gaussian process regression and conditional polynomial chaos for parameter estimation (2020)
- Zhang, Xinlei; Xiao, Heng; Gomez, Thomas; Coutier-Delgosha, Olivier: Evaluation of ensemble methods for quantifying uncertainties in steady-state CFD applications with small ensemble sizes (2020)
- Mara, Thierry A.; Fahs, Marwan; Shao, Qian; Younes, Anis: Random sampling from joint probability distributions defined in a Bayesian framework (2019)
- Wang, Chao; Chen, Qian; Gerlach, Richard: Bayesian realized-GARCH models for financial tail risk forecasting incorporating the two-sided Weibull distribution (2019)
- Karabatsos, George; Leisen, Fabrizio: An approximate likelihood perspective on ABC methods (2018)