Multivariate normal random values simulation on the simplex: Multivariate normal random values simulation on the simplex
Description
Multivariate normal random values simulation on the simplex.
Usage
rcompnorm(n, m, s, type = "alr")
Arguments
n
The sample size, a numerical value.
m
The mean vector in \(R^d\).
s
The covariance matrix in \(R^d\).
type
The alr (type = "alr") or the ilr (type = "ilr") is to be used for closing the Euclidean data onto the simplex.
Value
A matrix with the simulated data.
Details
The algorithm is straightforward, generate random values from a multivariate normal distribution in \(R^d\) and brings the
values to the simplex \(S^d\) using the inverse of a log-ratio transformation.
References
Aitchison J. (1986). The statistical analysis of compositional data. Chapman & Hall.