dmixdiri: Density values of a mixture of Dirichlet distributions

In Compositional: Compositional Data Analysis

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)

Compositional documentation built on Oct. 9, 2024, 5:10 p.m.