Fast sampling from β-ensembles

by   Guillaume Gautier, et al.

We study sampling algorithms for β-ensembles with time complexity less than cubic in the cardinality of the ensemble. Following Dumitriu Edelman (2002), we see the ensemble as the eigenvalues of a random tridiagonal matrix, namely a random Jacobi matrix. First, we provide a unifying and elementary treatment of the tridiagonal models associated to the three classical Hermite, Laguerre and Jacobi ensembles. For this purpose, we use simple changes of variables between successive reparametrizations of the coefficients defining the tridiagonal matrix. Second, we derive an approximate sampler for the simulation of β-ensembles, and illustrate how fast it can be for polynomial potentials. This method combines a Gibbs sampler on Jacobi matrices and the diagonalization of these matrices. In practice, even for large ensembles, only a few Gibbs passes suffice for the marginal distribution of the eigenvalues to fit the expected theoretical distribution. When the conditionals in the Gibbs sampler can be simulated exactly, the same fast empirical convergence is observed for the fluctuations of the largest eigenvalue. Our experimental results support a conjecture by Krishnapur et al. (2016), that the Gibbs chain on Jacobi matrices of size N mixes in 𝒪(log(N)).


page 1

page 2

page 3

page 4


Stable distributions and domains of attraction for unitarily invariant Hermitian random matrix ensembles

We consider random matrix ensembles on the Hermitian matrices that are h...

Convergence rate of a collapsed Gibbs sampler for crossed random effects models

In this paper, we analyze the convergence rate of a collapsed Gibbs samp...

Spectral Ergodicity in Deep Learning Architectures via Surrogate Random Matrices

In this work a novel method to quantify spectral ergodicity for random m...

Extended L-ensembles: a new representation for Determinantal Point Processes

Determinantal point processes (DPPs) are a class of repulsive point proc...

Approximating the Spectral Gap of the Pólya-Gamma Gibbs Sampler

The self-adjoint, positive Markov operator defined by the Pólya-Gamma Gi...

Backfitting for large scale crossed random effects regressions

Regression models with crossed random effect error models can be very ex...

Finite Mixtures of ERGMs for Ensembles of Networks

Ensembles of networks arise in many scientific fields, but currently the...

Please sign up or login with your details

Forgot password? Click here to reset