powered by
Dirichlet random values simulation.
rdiri(n, a)
The sample size, a numerical value.
A numerical vector with the parameter values.
A matrix with the simulated data.
The algorithm is straightforward, for each vector, independent gamma values are generated and then divided by their total sum.
Ng Kai Wang, Guo-Liang Tian and Man-Lai Tang (2011). Dirichlet and related distributions: Theory, methods and applications. John Wiley & Sons.
Aitchison J. (1986). The statistical analysis of compositional data. Chapman & Hall.
diri.est, diri.nr, diri.contour
# NOT RUN { x <- rdiri( 100, c(5, 7, 1, 3, 10, 2, 4) ) diri.est(x) # }
Run the code above in your browser using DataLab