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.

Independence properties of the truncated multivariate elliptical distributions

Some Theoretical and Computational Aspects of the Truncated Multivariate Skew-Normal/Independent Distributions

The multivariate tail-inflated normal distribution and its application in finance

The Asymptotic Distribution and Berry–Esseen Bound of a New Test for Independence in High Dimension with an Application to Stochastic Optimization

Multivariate double truncated expectation and covariance risk measures for elliptical distributions

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

Em Algorithm of the Truncated Multinormal Distribution with Linear Restriction on the Variables

Univariate likelihood projections and characterizations of the multivariate normal distribution

An analysis of multivariate measures of skewness and kurtosis of skew-elliptical distributions

A nonparametric test for elliptical distribution based on kernel embedding of probabilities

Some Results on the Truncated Multivariate Skew-Normal Distribution

Distribution-free tests of multivariate independence based on center-outward quadrant, Spearman, Kendall, and van der Waerden statistics

A Conditional Distribution Function-Based Measure for Independence and K-sample Tests in Multivariate Data

Multivariate doubly truncated moments for generalized skew-elliptical distributions with application to multivariate tail conditional risk measures

Testing High-dimensional Covariance Matrices under the Elliptical Distribution and Beyond

Dependence of multivariate extremes

RANK-BASED INDICES FOR TESTING INDEPENDENCE BETWEEN TWO HIGH-DIMENSIONAL VECTORS

Stochastic representations and probabilistic characteristics of multivariate skew-elliptical distributions

Multivariate tail covariance for generalized skew-elliptical distributions

A test for normality and independence based on characteristic function