Generate random draws from a Dirichlet distribution
rdir(n, alpha)
n
The number of observations.
alpha
A vector of scale parameters, such that E[p_j] = α_j/∑_iα_i.
vector
An n column matrix containing the observations.
Matt Taddy mataddy@gmail.com
rdir(3,rep(1,6))