R/DirRnd.R

DirRnd <- function(alpha)
## Dirichlet random numbers generator
{
    y <- rgamma(length(alpha), alpha)
    return(y/sum(y))
}
siamakz/BNPseq documentation built on May 15, 2019, 4:28 p.m.