Presented at 7th Valencia International Meeting on Bayesian Statistics, June 27, 2002, Tenerife, Spain

## Use of Probability Gradients in Hybrid MCMC and a New Convergence Test

Kenneth M. Hanson
Los Alamos National Laboratory

### Abstract

Based on a presentation given at the recent Valencia 7 meeting, this talk brings together several themes that are important in simulation codes, namely, Markov Chain Monte Carlo, hybrid MCMC, and adjoint differentiation of computer codes.

The hybrid Markov Chain Monte Carlo technique affords a robust means for sampling multidimensional probability density functions with high efficiency, provided one can calculate the gradient of phi = minus-log-probability in a time comparable to calculating the probability itself. The latter condition is met using the technique of adjoint differentiation. Consider a calculation of the probability in terms of a sequence of operations, which can be called the forward calculation. Then the implementation of the adjoint differentiation technique amounts to carrying out the chain rule for differentiation on that sequence, but in the reverse direction. Adjoint differentiation (AD) is very useful when the forward calculation involves a sequence of modules, which individually are simple, but collectively perform a complex calculation. Such is often the case for complex modeling situations, in which case it is recommended to implement AD using the computer code. See "Operation of the Bayes Inference Engine," K. M. Hanson et al., in Maximum Entropy and Bayesian Methods, pp. 309-318 (Kluwer Academic, 1999).

The gradients of phi may also be used to form a sensitive statistic to test the convergence of the MCMC sequence. The statistic is the ratio of sample estimates for the variance of the distribution calculated two different ways, one based on using integration by parts and the gradient of phi and the other based on the normal second-moment calculation. The efficiency of the hybrid MCMC technique and the usefulness of the convergence test will be demonstrated for simple multivariate normal distributions.

adjoint differentiation, Markov Chain Monte Carlo, hybrid MCMC, Hamiltonian method, Hamiltonian dynamics, leapfrog method, Metropolis MCMC, convergence test

Get presentation (pdf, 318 KB)