Convergence rates for a hierarchical Gibbs sampler

From MaRDI portal
(Redirected from Publication:502894)



Abstract: We establish some results for the rate of convergence in total variation of a Gibbs sampler to its equilibrium distribution. This sampler is motivated by a hierarchical Bayesian inference construction for a gamma random variable. Our results apply to a wide range of parameter values in the case that the hierarchical depth is 3 or 4, and are more restrictive for depth greater than 4. Our method involves showing a relationship between the total variation of two ordered copies of our chain and the maximum of the ratios of their respective co-ordinates. We construct auxiliary stochastic processes to show that this ratio does converge to 1 at a geometric rate.


Specific results are given for the rate of convergence of a Markov chain resulting from the application of Gibbs sampling to a hierarchical gamma model. In particular, observations are assumed to come from conditional iid sampling from a gamma distribution with fixed shape parameter and random scale parameter. In turn, the inverse scale parameter is assumed to be sampled from a higher level gamma distribution, and so on 3 or 4 times. The results are rather technical and the proofs rely on repeated applications of coupling and of properties of Markov processes. Explicit formulas for the convergence rates of Gibbs samplers are not common; here, the convergence considered is in the total variation metric.











This page was built for publication: Convergence rates for a hierarchical Gibbs sampler

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q502894)