Also, have a look at http://docs.pymc.io/notebooks/dp_mix.html if you have not done so. I think cell [26] and [37] in the notebook gives much better visualization of how good the model fitting and posterior prediction given the fit than just plotting the mean.