Simulation of compositional data from mixtures of Dirichlet distributions {Compositional}

Simulation of compositional data from mixtures of Dirichlet distributions

Description

Simulation of compositional data from mixtures of Dirichlet distributions.

Usage

rmixdiri(n, a, prob)

Arguments

n

The sample size.

a

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

prob

A vector with the mixing probabilities.

Details

A sample from a Dirichlet mixture model is generated.

Value

A list including:

id

A numeric variable indicating the cluster of simulated vector.

x

A matrix containing the simulated compositional data.

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

rmixcomp, 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)