Correct way to calculate the model evidence

Hi everyone,

I'm using an EnsembleSampler to fit a series of models to the data. For each model, I'm interested in an estimate of the best fitting parameters, but I'd also like to calculate the odds ratios for all the models in order to understand which one is better describing the data.
To do so I need to recover the evidence, i.e. the integral of the un-normalized posterior. Is there a way to recover this information from the output of run_mcmc() when using an EnsembleSampler?
I know that PTSampler can do it but I saw that it was moved to another package. Since I'm working on a pre-existing code, I'd prefer not to adapt it to use the new ptemcee package if there is a reasonable way to get these data using the current implementation of the algorithm.

Thanks in advance for the help,

Cheers,
Enrico

The issue names EnsembleSampler, run_mcmc(), PTSampler, and ptemcee but does not mention repository files or tests. Start by checking whether run_mcmc() exposes evidence-related output; a contribution would need defined implementation and validation criteria, which the issue does not provide.