Density of the folded normal distribution {Compositional}

Density of the folded model normal distribution

Description

Density of the folded model normal distribution.

Usage

dfolded(x, a, p, mu, su, logged = TRUE)

Arguments

x

A vector or a matrix with compositional data. No zeros are allowed.

a

The value of \alpha.

p

The probability inside the simplex of the folded model.

mu

The mean vector.

su

The covariance matrix.

logged

A boolean variable specifying whether the logarithm of the density values to be returned. It is set to TRUE by default.

Details

Density values of the folded model.

Value

The density value(s).

Author(s)

Michail Tsagris.

R implementation and documentation: Michail Tsagris [e-mail].

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

rfolded, a.est, folded.contour

Examples

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)
mod <- a.est(x)
den <- dfolded(x, mod$best, mod$p, mod$mu, mod$su)