Added a function to simulate molecular rates under the dirichlet process.