R/rdirichlet.R

In ricky-kotecha/rkHMM:

rdirichlet <- function(n, alpha) {
  # taken from gtools
  l <- length(alpha)
  x <- matrix(rgamma(l * n, alpha), ncol = l, byrow = TRUE)
  sm <- x %*% rep(1, l)
  res  <- x / as.vector(sm)
  res[is.nan(res)] <- 0
  res
}