DirRnd <- function(alpha) ## Dirichlet random numbers generator { y <- rgamma(length(alpha), alpha) return(y/sum(y)) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.