MCMC方法在贝叶斯统计中运用很多，MIT发布的EMCEE是实现的比较好的。介绍页面在下面。源代码中examples里的代码可以帮助理解各种功能，特别是line.py 列出了最小二乘法，最大似然法和MCMC方法进行线性拟合的测试结果。
此方法最重要的问题是需要会按照自己的需要改写似然函数。参考文献：1008.4686
http://dan.iel.fm/emcee/current/
MCMC方法的科普介绍可以参考
http://www.cnblogs.com/ywl925/archive/2013/06/05/3118875.html