auto_keq fix and MCMC fixes

auto_keq now works with index

There are also many MCMC fixes. numpy.percentile has been removed and instead arviz.hdi is used. This handles asymmetric distributions much better. Generation of the KDE error model doesn't remove outliers now. Due to how the sampling is done the outliers are not actually outliers.