Density function and random number generation for the Dirichlet distribution

1 | ```
rdirichlet(n, alpha)
``` |

`n` |
number of random observations to draw |

`alpha` |
the Dirichlet distribution's parameters. Can be a vector (one set of parameters for all observations) or a matrix (a different set of parameters for each observation), see “Details” |

The Dirichlet distribution is a multidimensional generalization of the Beta distribution where each dimension is governed by an *alpha*-parameter.
Formally this is

*
D(α)=[Γ(∑α)/∏Γ(α)]∏ y^(α-1)
*

Usually, `alpha`

is a vector thus the same parameters will be used for all observations.
If `alpha`

is a matrix, a complete set of *alpha*-parameters must be supplied for each observation.

returns a matrix with random numbers according to the supplied alpha vector or matrix.

Chong Wu

1 2 | ```
X1 <- rdirichlet(100, c(5, 5, 10))
X1
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.