ENH: Dirichlet distribution

Implements a dirichlet distribution in the stats module.