Ilya

HI Inti! You can try to initialize walkers near maximum likelihood values of the parameters with `emcee.utils.sampling_ball`. If your posterior density of the parameters has many modes then try parallel...

Really, Dan, could you give an example to emcee users somewhere? That issue with stochastic calculations inside emcee is a bit confusing... Thanks!

Hi! Check this [quastion](https://stackoverflow.com/q/35843519) on SO. This occurs in ``scipy.cluster.kmeans2``. Thus it is not ``nestle``-specific

Here the author explains why he has implemented sampyl and what is a problem with pymc3: http://matatat.org/ab-testing-with-sampyl.html

@Gabriel-p , there's nice and easy to use python implementation of the nested sampling [here](https://github.com/kbarbary/nestle) by @kbarbary. So if you are not sticked to MCMC and want to get evidence...

@madanh , what is the problem? [Check here](https://stackoverflow.com/questions/6022302/how-to-apply-unmerged-upstream-pull-requests-from-other-forks-into-my-fork) Edit: Oops, it is your's PR...