To avoid concerns about MCMC convergence we can change variables, integrate out two parameters, and grid search over the remaining ones. See: - **Notebook entry on [Analytic calculations](http://carlboettiger.info/2013/06/04/analytic-marginalizing-for-posteriors.html)** which change variables and integrate - **Notebook entry on [Grid-based search](http://carlboettiger.info/2013/06/05/semi-analytic-posteriors.html)** The latter still a work in progress, being developed in analysis script: [semi-analytic-posteriors](https://github.com/cboettig/nonparametric-bayes/blob/master/inst/examples/semi-analytic-posteriors.md)