implemented mcmc for inference (scipy.optim BFGS too unstable)