Generates random variables from a Dirichlet sampling model.

rvdirichlet(n = 1, alpha)

Arguments

n

integer: number of vectors to generate

alpha

the parameter vector; may be random

Value

A random vector (rv object) of length n.

Details

The Dirichlet distribution is a generalization of the Beta distribution. (If alpha is of length two, rvdirichlet draws from the Beta model.)

References

Kerman, J. and Gelman, A. (2007). Manipulating and Summarizing Posterior Simulations Using Random Variable Objects. Statistics and Computing 17:3, 235-244.

See also vignette("rv").

Author

Jouni Kerman jouni@kerman.com

Examples

if (FALSE) {

  a <- rvdirichlet(1, alpha=c(6, 3, 1)) # 
  sum(a) # one with probability 1   
  }