Related papers: Fast sampling from $\beta$-ensembles

This work identifies a solvable (in the sense that spectral correlation functions can be expressed in terms of orthogonal polynomials), rotationally invariant random matrix ensemble with a logarithmic weakly confining potential. The…

Using the Coulomb gas method and standard methods of statistical physics, we compute analytically the joint cumulative probability distribution of the extreme eigenvalues of the Jacobi-MANOVA ensemble of random matrices, in the limit of…

There are several methods to treat ensembles of random matrices in symmetric spaces, circular matrices, chiral matrices and others. Orthogonal polynomials and the supersymmetry method are particular powerful techniques. Here, we present a…

The focus of this paper is on the probability, $E_\beta(0;J)$, that a set $J$ consisting of a finite union of intervals contains no eigenvalues for the finite $N$ Gaussian Orthogonal ($\beta=1$) and Gaussian Symplectic ($\beta=4$) Ensembles…

Here, using two real non-zero parameters $\lambda$ and $\mu$, we construct Gaussian pseudo-orthogonal ensembles of a large number $N$ of $n \times n$ ($n$ even and large) real pseudo-symmetric matrices under the metric $\eta$ using $…

Quantum dynamical simulations of statistical ensembles pose a significant computational challenge due to the fact that mixed states need to be represented. If the underlying dynamics is fully unitary, for example in ultrafast coherent…

We give a proof of the Universality Conjecture for orthogonal (beta=1) and symplectic (beta=4) random matrix ensembles of Laguerre-type in the bulk of the spectrum as well as at the hard and soft spectral edges. Our results are stated…

Let $\mathcal{P}_{\beta}^{(V)} (N_{\cal I})$ be the probability that a $N\times N$ $\beta$-ensemble of random matrices with confining potential $V(x)$ has $N_{\cal I}$ eigenvalues inside an interval ${\cal I}=[a,b]$ of the real line. We…

We study the global spectrum fluctuations for $\beta$-Hermite and $\beta$-Laguerre ensembles via the tridiagonal matrix models introduced in \cite{dumitriu02}, and prove that the fluctuations describe a Gaussian process on monomials. We…

This paper presents a Jacobi-type iteration for computing a given specified eigenpair of a symmetric matrix. For a certain class of diagonally dominant matrices, the procedure is shown to converge at a linear rate depending on how the…

The well-known Gumbel-Max Trick for sampling elements from a categorical distribution (or more generally a non-negative vector) and its variants have been widely used in areas such as machine learning and information retrieval. To sample a…

When performing Bayesian data analysis using a general linear mixed model, the resulting posterior density is almost always analytically intractable. However, if proper conditionally conjugate priors are used, there is a simple two-block…

We compute the gap probability that a circle of radius r around the origin contains exactly k complex eigenvalues. Four different ensembles of random matrices are considered: the Ginibre ensembles and their chiral complex counterparts, with…

We study the problem of approximating the eigenspectrum of a symmetric matrix $\mathbf A \in \mathbb{R}^{n \times n}$ with bounded entries (i.e., $\|\mathbf A\|_{\infty} \leq 1$). We present a simple sublinear time algorithm that…

The concept of probabilistic values, such as Beta Shapley values and weighted Banzhaf values, has gained recent attention in applications like feature attribution and data valuation. However, exact computation of these values is often…

The speed of convergence of the R-linear GMRES is bounded in terms of a polynomial approximation problem on a finite subset of the spectrum. This result resembles the classical GMRES convergence estimate except that the matrix involved is…

The MC$^3$ (Madigan and York, 1995) and Gibbs (George and McCulloch, 1997) samplers are the most widely implemented algorithms for Bayesian Model Averaging (BMA) in linear regression models. These samplers draw a variable at random in each…

We study the rate of convergence for the largest eigenvalue distributions in the Gaussian unitary and orthogonal ensembles to their Tracy-Widom limits. We show that one can achieve an $O(N^{-2/3})$ rate with particular choices of the…

The emergence of big data has led to so-called convergence complexity analysis, which is the study of how Markov chain Monte Carlo (MCMC) algorithms behave as the sample size, $n$, and/or the number of parameters, $p$, in the underlying…

We consider properties of determinants of some random symmetric matrices issued from multivariate statistics: Wishart/Laguerre ensemble (sample covariance matrices), Uniform Gram ensemble (sample correlation matrices) and Jacobi ensemble…