old MCMC stuff now working
Model selection with Laplace approximation round 1