Adding pure dirichlet estimator