We give the first analysis of a systematic scan version of the Metropolis algorithm. Our examples include generating random elements of a Coxeter group with probability determined by the length function. The analysis is based on interpreting Metropolis walks in terms of the multiplication in the Iwahori-Hecke algebra.

This paper presents a geometric approach to estimating subspaces as elements of complex Grassmann-manifold, with each subspace represented by its unique, complex projection matrix. Variation between the subspaces is modeled by rotating their projection matrices via the action of unitary matrices [elements of the unitary group U( )]. Subspace estimation or tracking then corresponds to inferences...

We give the first analysis of a systematic scan version of the Metropolis algorithm. Our examples include generating random elements of a Coxeter group with probability determined by the length function. The analysis is based on interpreting Metropolis walks in terms of the multiplication in the Iwahori-Hecke algebra.

The Metropolis Light Transport algorithm is a variant of the classic Metropolis method used in statistical physics. A variance analysis of the Metropolis Light Transport algorithm is presented that bounds its variance in terms of the number of paths used and the intrinsic correlation between samples. It is shown that the variance of a pixel is where is the number of samples for the entire image...

(First draft March 2005; revised November 2005) Abstract This paper extends some adaptive schemes that have been developed for the Random Walk Metropolis algorithm to more general versions of the Metropolis-Hastings (MH) algorithm, particularly to the Metropolis Adjusted Langevin algorithm of Roberts and Tweedie (1996). Our simulations show that the adaptation drastically improves the performan...

We introduce an adaptive output-sensitive Metropolis-Hastings algorithm for probabilistic models expressed as programs, Adaptive Lightweight Metropolis-Hastings (AdLMH). The algorithm extends Lightweight Metropolis-Hastings (LMH) by adjusting the probabilities of proposing random variables for modification to improve convergence of the program output. We show that AdLMH converges to the correct...

This paper describes modeling discrete event systems in Metropolis using one of the key concepts employed by Metropolis, orthogonalization of design aspects, which in this particular case, is the orthogonalization between capability and cost. To support the orthogonalization, quantity annotation mechanism is introduced. This paper formally analyzes simulation strategies for quantity annotation ...

We show that sampling with a biased Metropolis scheme is essentially equivalent to using the heatbath algorithm. However, the biased Metropolis method can also be applied when an efficient heatbath algorithm does not exist. This is first illustrated with an example from high energy physics (lattice gauge theory simulations). We then illustrate the Rugged Metropolis method, which is based on a s...

We consider Markov chain Monte Carlo algorithms which combine Gibbs updates with Metropolis-Hastings updates, resulting in a conditional Metropolis-Hastings sampler. We develop conditions under which this sampler will be geometrically or uniformly ergodic. We apply our results to an algorithm for drawing Bayesian inferences about the entire sample path of a diffusion process, based only upon di...

This paper gives geometric tools: comparison, Nash and Sobolev inequalities for pieces of the relevent Markov operators, that give useful bounds on rates of convergence for the Metropolis algorithm. As an example, we treat the random placement of N hard discs in the unit square, the original application of the Metropolis algorithm.