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.

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

Some Results on the Truncated Multivariate Skew-Normal Distribution

Independence properties of the truncated multivariate elliptical distributions

The alchemy of probability distributions: beyond Gram-Charlier expansions, and a skew-kurtotic-normal distribution from a rank transmutation map

Probability boxes on totally preordered spaces for multivariate modelling

A General Inferential Framework for Singly-Truncated Bivariate Normal Models with Applications in Economics

On Moments of Folded and Doubly Truncated Multivariate Extended Skew-Normal Distributions

Fast Simulation of Hyperplane-Truncated Multivariate Normal Distributions

Multivariate Geometric Skew-Normal Distribution

Some properties of the unified skew-normal distribution

On Bayesian predictive density estimation for skew-normal distributions

Moments Calculation for the Doubly Truncated Multivariate Normal Density

A Tractable Probability Distribution with Applications in Three-Dimensional Statistics

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

Statistical applications of the multivariate skew normal distribution

Statistical applications of the multivariate skew-normal distribution

Inference in Randomized Least Squares and PCA via Normality of Quadratic Forms

Mean mixtures of normal distributions: properties, inference and application

Generating Multivariate Ordinal Data via Entropy Principles

Scale mixtures of multivariate centered skew-normal distributions

The orthogonal skew model: computationally efficient multivariate skew-normal and skew-t distributions with applications to model-based clustering