This MCMC method takes a data numeric vector (Y) and assigns the elements of Y
to a (potentially infinite) number of normal distributions. The individual normal distributions from a mixture of normals can be inferred.
Following the method described in Escobar (1994)