Abstract
We characterize the convergence of the Gibbs sampler which samples from the joint posterior distribution of parameters and missing data in hierarchical linear models with arbitrary symmetric error distributions. We show that the convergence can be uniform, geometric or subgeometric depending on the relative tail behavior of the error distributions, and on the parametrization chosen. Our theory is applied to characterize the convergence of the Gibbs sampler on latent Gaussian process models. We indicate how the theoretical framework we introduce will be useful in analyzing more complex models.
Citation
Omiros Papaspiliopoulos. Gareth Roberts. "Stability of the Gibbs sampler for Bayesian hierarchical models." Ann. Statist. 36 (1) 95 - 117, February 2008. https://doi.org/10.1214/009053607000000749
Information