Abstract
We present a polynomial-time perfect sampler for the Q-Ising with a vertex-independent noise. The Q-Ising, one of the generalized models of the Ising, arose in the context of Bayesian image restoration in statistical mechanics. We study the distribution of Q-Ising on a two-dimensional square lattice over n vertices, that is, we deal with a discrete state space {1,…,Q}n for a positive integer Q. Employing the Q-Ising (having a parameter β) as a prior distribution, and assuming a Gaussian noise (having another parameter α), a posterior is obtained from the Bayes’ formula. Furthermore, we generalize it: the distribution of noise is not necessarily a Gaussian, but any vertex-independent noise. We first present a Gibbs sampler from our posterior, and also present a perfect sampler by defining a coupling via a monotone update function. Then, we show O(nlog n) mixing time of the Gibbs sampler for the generalized model under a condition that β is sufficiently small (whatever the distribution of noise is). In case of a Gaussian, we obtain another more natural condition for rapid mixing that α is sufficiently larger than β. Thereby, we show that the expected running time of our sampler is O(nlog n).
Similar content being viewed by others
References
Bubley R, Dyer M (1997) Path coupling: a technique for proving rapid mixing in Markov chains. In: Proceedings of the 38th annual symposium on foundations of computer science (FOCS 1997), pp 223–231
Geman S, Geman D (1984) Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images. IEEE Trans Pattern Anal Mach Intell 6:721–741
Gibbs AL (2000) Bounding the convergence time of the Gibbs sampler in Bayesian image restoration. Biometrika 87(4):749–766
Gibbs AL (2004) Convergence in the Wasserstein metric for Markov chain Monte Carlo algorithms with applications to image restoration. Stochastic Models 20:473–492
Inoue J, Carlucci DM (2001) Image restoration using the Q-Ising spin glass. Phys Rev E 64:036121
Levin DA, Peres Y, Wilmer EL (2008) Markov chains and mixing times. American Mathematical Society, Providence
Propp J, Wilson D (1996) Exact sampling with coupled Markov chains and applications to statistical mechanics. Random Struct Algorithms 9:223–252
Tanaka K, Inoue J, Titterington DM (2003) Probabilistic image processing by means of the Bethe approximation for the Q-Ising model. J Phys A, Math Gen 36:11023–11035
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Yamamoto, M., Kijima, S. & Matsui, Y. A polynomial-time perfect sampler for the Q-Ising with a vertex-independent noise. J Comb Optim 22, 392–408 (2011). https://doi.org/10.1007/s10878-010-9309-7
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10878-010-9309-7