http://andrewgelman.com/2012/12/06/yes-checking-calibration-of-probability-forecasts-is-part-of-bayesian-statistics/

E(y.new | p.hat) = p.hat for any p.hat


[Binary logistic regression with rstanarm](https://avehtari.github.io/modelselection/diabetes.html)