An MCMC Method to Sample from Lattice Distributions

Abstract

We introduce a Markov Chain Monte Carlo (MCMC) algorithm to generate samples from probability distributions supported on a d-dimensional lattice = BZd, where B is a full-rank matrix. Specifically, we consider lattice distributions P in which the probability at a lattice point is proportional to a given probability density function, f, evaluated at that point. To generate samples from P, it suffices to draw samples from a pull-back measure PZd defined on the integer lattice. The probability of an integer lattice point under PZd is proportional to the density function π = |(B)|f B. The algorithm we present in this paper for sampling from PZd is based on the Metropolis-Hastings framework. In particular, we use π as the proposal distribution and calculate the Metropolis-Hastings acceptance ratio for a well-chosen target distribution. We can use any method, denoted by ALG, that ideally draws samples from the probability density π, to generate a proposed state. The target distribution is a piecewise sigmoidal distribution, chosen such that the coordinate-wise rounding of a sample drawn from the target distribution gives a sample from PZd. When ALG is ideal, we show that our algorithm is uniformly ergodic if -(π) satisfies a gradient Lipschitz condition.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…