Density values of a mixture of Dirichlet distributions {Compositional}

Density values of a mixture of Dirichlet distributions

Description

Density values of a mixture of Dirichlet distributions.

Usage

dmixdiri(x, a, prob, logged = TRUE)

Arguments

x

A vector or a matrix with compositional data. Zeros are not allowed.

a

A matrix where each row contains the parameters of each Dirichlet component.

prob

A vector with the mixing probabilities.

logged

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

Details

The density of the mixture of Dirichlet distribution for a vector or a matrix of compositional data is returned.

Value

A vector with the density values.

Author(s)

Michail Tsagris.

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

References

Ye X., Yu Y. K. and Altschul S. F. (2011). On the inference of Dirichlet mixture priors for protein sequence comparison. Journal of Computational Biology, 18(8), 941-954.

See Also

rmixdiri, mixdiri.contour

Examples

a <- matrix( c(12, 30, 45, 32, 50, 16), byrow = TRUE,ncol = 3)
prob <- c(0.5, 0.5)
x <- rmixdiri(100, a, prob)$x
f <- dmixdiri(x, a, prob)