Compositional (version 5.5)

Simulation of compositional data from mixtures of Dirichlet distributions: 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.

Value

A list including:

id

A numeric variable indicating the cluster of simulated vector.

x

A matrix containing the simulated compositional data.

Details

A sample from a Dirichlet mixture model is generated.

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

# 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)
# }