Number of steps cut with fixed parameters

[tomasstolker]

Number of steps cut with fixed parameters

[tomasstolker]

Thanks for pointing out that issue @vandalt!

Right, this is something that I had implemented recently. The problem is that the ndim of the stored samples from run_mcmc is 3, while the ndim from run_multinest is 2. In the latter case the shape if (nsamples, nparams).

It would be great if you could create a PR which solves this issues and takes the two ndim cases into account.

By the way, when using FitModel, the run_multinest function has more features compared to run_mcmc and convergences more easily in case the posterior is multimodal. It also returns the evidence which can be used for model comparison.

[vandalt]

Right, I had not tried run_multinest so far so I did not think about that part, thank you.

I will create a PR that handles both cases. Ellipsis would probably be the best way to achieve this (i.e. samples[..., i]).

Thanks for the tip about run_multinest vs run_mcmc.