initial Dirichlet Process code