R/rgendiri.R

In Compositional: Compositional Data Analysis

Defines functions rgendiri

Documented in rgendiri

rgendiri <- function(n, a, b) {
  ##  n is the sample size
  ##  a is the parameters vector
  D <- length(a)
  y <- rangen::colRgamma(n, D, shape = a, rate = b)
  y / Rfast::rowsums(y)  ## Dirichlet simulated values
}

Try the Compositional package in your browser