When a graph can be decomposed into clusters of well connected subgraphs, it is possible to speed up random walks taking advantage of the topology of the graph. In this work, a new random walk scheme is introduced and a condition is given when the new random walk performs better than the Metropolis algorithm.

We illustrate for 4D SU(2) and U(1) lattice gauge theory that sampling with a biased Metropolis scheme is essentially equivalent to using the heat bath algorithm. Only, the biased Metropolis method can also be applied when an efficient heat bath algorithm does not exist. For the examples discussed the biased Metropolis algorithm is also better suited for parallelization than the heat bath algor...

The waste-recycling Monte Carlo (WR) algorithm introduced by physicists is a modification of the (multi-proposal) Metropolis-Hastings algorithm, which makes use of all the proposals in the empirical mean, whereas the standard (multi-proposal) Metropolis-Hastings algorithm only uses the accepted proposals. In this paper, we extend the WR algorithm into a general control variate technique and exh...

Monte Carlo (MC) sampling methods are widely applied in Bayesian inference, system simulation and optimization problems. The Markov Chain Monte Carlo (MCMC) algorithms are a well-known class of MC methods which generate a Markov chain with the desired invariant distribution. In this document, we focus on the Metropolis-Hastings (MH) sampler, which can be considered as the atom of the MCMC techn...

Necessary and su cient conditions for geometric con vergence in the relative supremum norm of the Metropolis Hastings simulation algorithm with a general generating function are estab lished An explicit expression for the convergence rate is given Introduction This paper discusses the convergence rate for the Metropolis Hastings simulation algorithm proposed in Hastings The Metropolis Hastings ...

The adaptive Metropolis (AM) algorithm of Haario, Saksman and Tamminen [Bernoulli 7 (2001) 223-242] uses the estimated covariance of the target distribution in the proposal distribution. This paper introduces a new robust adaptive Metropolis algorithm estimating the shape of the target distribution and simultaneously coercing the acceptance rate. The adaptation rule is computationally simple ad...

Estimation of small failure probabilities is one of the most important and challenging problems in reliability engineering. In cases of practical interest, the failure probability is given by a high-dimensional integral. Since multivariate integration suffers from the curse of dimensionality, the usual numerical methods are inapplicable. Over the past decade, the civil engineering research comm...

Autoregressive fractional integrated moving average modeling strategy was used to model the daily average temperature (DAT) series of Sokoto metropolis for the period of 01/01/2003 to 03/04/2007. The time plot suggests that there is persistence dependence in the series. The order of fractional integration was found to be 0.6238841. The correct model for the daily average temperature data (DAT) ...