changed order for speed of fitting
does posterior predictive checking as log density and difference on y axis