R program of Dirichlet process