1. Implement parallel MCMC?
2. Include posterior variance estimates for regression change point.
3. Produce function to give modal partition (if return.mcmc=TRUE), maybe?