Abstract:Compositional data, which is data consisting of fractions or probabilities, is common in many fields including ecology, economics, physical science and political science. If these data would otherwise be normally distributed, their spread can be conveniently represented by a multivariate normal distribution truncated to the non-negative space under a unit simplex. Here this distribution is called the simplex-truncated multivariate normal distribution. For calculations on truncated distributions, it is often useful to obtain rapid estimates of their integral, mean and covariance; these quantities characterising the truncated distribution will generally possess different values to the corresponding non-truncated distribution. In this paper, three different approaches that can estimate the integral, mean and covariance of any simplex-truncated multivariate normal distribution are described and compared. These three approaches are (1) naive rejection sampling, (2) a method described by Gessner et al. that unifies subset simulation and the Holmes-Diaconis-Ross algorithm with an analytical version of elliptical slice sampling, and (3) a semi-analytical method that expresses the integral, mean and covariance in terms of integrals of hyperrectangularly-truncated multivariate normal distributions, the latter of which are readily computed in modern mathematical and statistical packages. Strong agreement is demonstrated between all three approaches, but the most computationally efficient approach depends strongly both on implementation details and the dimension of the simplex-truncated multivariate normal distribution. For computations in low-dimensional distributions, the semi-analytical method is fast and thus should be considered. As the dimension increases, the Gessner et al. method becomes the only practically efficient approach of the methods tested here.

Approximate null distribution of the largest root in multivariate analysis

Accuracy of the Tracy--Widom limits for the extreme eigenvalues in white Wishart matrices

Computing the Tracy--Widom Distribution for Arbitrary

Computing the Tracy-Widom Distribution for Arbitrary $β>0$

An asymptotic refinement of the Gauss-Lucas Theorem for random polynomials with i.i.d. roots

A negative binomial approximation to the distribution of the sum of maxima of indicator random variables

Expected Euler characteristic method for the largest eigenvalue: (Skew-)orthogonal polynomial approach

A Tracy-Widom Empirical Estimator For Valid P-values With High-Dimensional Datasets

Optimal estimation of the null distribution in large-scale inference

Almost Worst Case Distributions in Multiple Priors Models

Large deviations of extreme eigenvalues of generalized sample covariance matrices

Asymptotics for non-degenerate multivariate $U$-statistics with estimated nuisance parameters under the null and local alternative hypotheses

Distributions of differences of Riemann zeta zeros

Joint distribution of primes in multiple short intervals

Large deviations for the largest singular value of sparse non-Hermitian matrices

Maximum interpoint distance of high-dimensional random vectors

On the Null Distribution of Bayes Factors in Linear Regression

Fluctuation of the Largest Eigenvalue of a Kernel Matrix with application in Graphon-based Random Graphs

Integral, mean and covariance of the simplex-truncated multivariate normal distribution

Optimal tail estimates in $β$-ensembles and applications to last passage percolation

A Fast Algorithm for Computing Macaulay Null Spaces of Bivariate Polynomial Systems