HR: 1700h
AN: H14A-05 [Abstracts]
TI: Speed-up of Markov Chain Monte Carlo Simulation Using Self-Adaptive Different Evolution with Subspace Sampling
AU: * Vrugt, J A
EM: vrugt@lanl.gov
AF: Center for NonLinear Studies (CNLS), Los Alamos National Laboratory, Mail Stop T003,
Los Alamos, NM 87545,
AB:
Markov chain Monte Carlo (MCMC) methods are widely used in fields ranging from physics and chemistry, to
finance, economics and statistical inference for estimating the average properties of complex systems. The
convergence rate of MCMC schemes is often observed, however to be disturbingly low, limiting its practical use in
many applications. This is frequently caused by an inappropriate selection of the proposal distribution used to
generate trial moves. Here we show that significant improvements to the efficiency of MCMC algorithms can be
made by using a self-adaptive Differential Evolution search strategy within a population-based evolutionary
framework. This scheme differs fundamentally from existing MCMC algorithms, in that trial jumps are simply a
fixed multiple of the difference of randomly chosen members of the population using various genetic operators
that are adaptively updated during the search. In addition, the algorithm includes randomized subspace sampling
to further improve convergence and acceptance rate. Detailed balance and ergodicity of the algorithm are proved,
and hydrologic examples show that the proposed method significantly enhances the efficiency and applicability of
MCMC simulations to complex, multi-modal search problems.
DE: 1846 Model calibration (3333)
DE: 1873 Uncertainty assessment (3275)
DE: 1894 Instruments and techniques: modeling
SC: Hydrology [H]
MN: 2007 Fall Meeting