Efficient sampling from the Bingham distribution

09/30/2020
by   Rong Ge, et al.
0

We give a algorithm for exact sampling from the Bingham distribution p(x)∝(x^⊤ A x) on the sphere 𝒮^d-1 with expected runtime of poly(d, λ_max(A)-λ_min(A)). The algorithm is based on rejection sampling, where the proposal distribution is a polynomial approximation of the pdf, and can be sampled from by explicitly evaluating integrals of polynomials over the sphere. Our algorithm gives exact samples, assuming exact computation of an inverse function of a polynomial. This is in contrast with Markov Chain Monte Carlo algorithms, which are not known to enjoy rapid mixing on this problem, and only give approximate samples. As a direct application, we use this to sample from the posterior distribution of a rank-1 matrix inference problem in polynomial time.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
05/16/2020

Polynomial-time approximation algorithms for the antiferromagnetic Ising model on line graphs

We present a polynomial-time Markov chain Monte Carlo algorithm for esti...
research
03/08/2017

Polynomial Time Algorithms for Dual Volume Sampling

We study dual volume sampling, a method for selecting k columns from an ...
research
11/26/2019

Approximating the Permanent by Sampling from Adaptive Partitions

Computing the permanent of a non-negative matrix is a core problem with ...
research
12/10/2021

Sampling from Discrete Energy-Based Models with Quality/Efficiency Trade-offs

Energy-Based Models (EBMs) allow for extremely flexible specifications o...
research
04/22/2023

Posterior Sampling from the Spiked Models via Diffusion Processes

Sampling from the posterior is a key technical problem in Bayesian stati...
research
05/30/2022

Notes on the runtime of A* sampling

The challenge of simulating random variables is a central problem in Sta...
research
10/20/2018

A Polynomial Time MCMC Method for Sampling from Continuous DPPs

We study the Gibbs sampling algorithm for continuous determinantal point...

Please sign up or login with your details

Forgot password? Click here to reset