R/dirichlet.simul.R

Defines functions dirichlet.simul

Documented in dirichlet.simul

## File Name: dirichlet.simul.R
## File Version: 0.131


#-- simulate from a Dirichlet distribution
dirichlet.simul <- function( alpha )
{
    N <- nrow(alpha)
    K <- ncol(alpha)
    ygamma <- 0*alpha
    for (ii in 1L:K){
        ygamma[,ii] <- stats::rgamma( n=N, shape=alpha[,ii] )
    }
    x <- ygamma / rowSums(ygamma)
    return(x)
}
alexanderrobitzsch/sirt documentation built on March 18, 2024, 1:29 p.m.