DREAM(D): an adaptive Markov Chain Monte Carlo simulation algorithm to solve discrete, noncontinuous, and combinatorial posterior parameter estimation problems.

J.A. Vrugt, C.J.F. ter Braak

Research output: Contribution to journalArticleAcademicpeer-review

124 Citations (Scopus)

Abstract

Formal and informal Bayesian approaches have found widespread implementation and use in environmental modeling to summarize parameter and predictive uncertainty. Successful implementation of these methods relies heavily on the availability of efficient sampling methods that approximate, as closely and consistently as possible the (evolving) posterior target distribution. Much of this work has focused on continuous variables that can take on any value within their prior defined ranges. Here, we introduce theory and concepts of a discrete sampling method that resolves the parameter space at fixed points. This new code, entitled DREAM(D) uses the recently developed DREAM algorithm (Vrugt et al., 2008, 2009a, b) as its main building block but implements two novel proposal distributions to help solve discrete and combinatorial optimization problems. This novel MCMC sampler maintains detailed balance and ergodicity, and is especially designed to resolve the emerging class of optimal experimental design problems. Three different case studies involving a Sudoku puzzle, soil water retention curve, and rainfall – runoff model calibration problem are used to benchmark the performance of DREAM(D). The theory and concepts developed herein can be easily integrated into other (adaptive) MCMC algorithms.
Original languageEnglish
Pages (from-to)3701-3713
JournalHydrology and Earth System Sciences
Volume15
Issue number12
DOIs
Publication statusPublished - 2011

Keywords

  • rainfall-runoff models
  • metropolis algorithm
  • differential evolution
  • global optimization
  • uncertainty
  • spaces
  • mcmc

Fingerprint

Dive into the research topics of 'DREAM(D): an adaptive Markov Chain Monte Carlo simulation algorithm to solve discrete, noncontinuous, and combinatorial posterior parameter estimation problems.'. Together they form a unique fingerprint.

Cite this