Compositional (version 5.5)

Density values of a mixture of Dirichlet distributions: 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.

Value

A vector with the density values.

Details

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

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

# NOT RUN {
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)
# }