rmixdiri: Simulation of compositional data from mixtures of Dirichlet...

In Compositional: Compositional Data Analysis

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)

Compositional documentation built on Oct. 23, 2023, 5:09 p.m.