TITLE:

IMPLEMENTING COMPONENTWISE HASTINGS ALGORITHMS (No. 73)

DATE:

Friday, March 4th, 2005

TIME:

3:30 PM

LOCATION:

GMCS 214

SPEAKER:

Richard Levine, Department of Mathematics and Statistics, San Diego State University

ABSTRACT:

Markov chain Monte Carlo (MCMC) routines have revolutionized the
application of Monte Carlo methods in statistical application and
statistical computing methodology. The Hastings sampler, encompassing both the Gibbs and Metropolis samplers as special cases, is the most commonly applied MCMC algorithm. The performance of the Hastings sampler relies heavily on the choice of sweep strategy, that is, the method by which the components or blocks of the random variable X of interest are visited and updated, and the choice of proposal distribution, that is the distribution from which candidate variates be drawn for the accept-reject rule each iteration of the algorithm. We focus on the random sweep strategy, where the components are updated in a random order, and random proposal
distributions, where the proposal distribution is characterized by
a randomly generated parameter. We develop an adaptive Hastings
sampler which learns from and adapts to random variates generated during the algorithm towards choosing the optimal random sweep strategy and proposal distribution for the problem at hand. As part of the development, we prove convergence of the random variates to the
distribution of interest and discuss practical implementations of the
methods. We illustrate the results presented by applying the adaptive
componentwise Hastings samplers developed to sample multivariate Gaussian target distributions and Bayesian frailty models.

HOST:

Jose Castillo

DOWNLOAD: