Learn R Programming

Compositional (version 5.5)

Simulation of compositional data from the folded normal distribution: Simulation of compositional data from the folded model normal distribution

Description

Simulation of compositional data from the folded model normal distribution.

Usage

rfolded(n, mu, su, a)

Arguments

n

The sample size.

mu

The mean vector.

su

The covariance matrix.

a

The value of \(\alpha\).

Value

A matrix with compositional data.

Details

A sample from the folded model is generated.

References

Tsagris M. and Stewart C. (2020). A folded model for compositional data analysis. Australian and New Zealand Journal of Statistics, 62(2): 249-277. https://arxiv.org/pdf/1802.07330.pdf

See Also

alfa, alpha.mle, a.est

Examples

Run this code
# NOT RUN {
s <-  c(0.1490676523, -0.4580818209,  0.0020395316, -0.0047446076, -0.4580818209,
1.5227259250,  0.0002596411,  0.0074836251,  0.0020395316,  0.0002596411,
0.0365384838, -0.0471448849, -0.0047446076,  0.0074836251, -0.0471448849,
0.0611442781)
s <- matrix(s, ncol = 4)
m <- c(1.715, 0.914, 0.115, 0.167)
x <- rfolded(100, m, s, 0.5)
a.est(x)
# }

Run the code above in your browser using DataLab