bilby.core.sampler.dynesty.dynesty_stats_plot

bilby.core.sampler.dynesty.dynesty_stats_plot(sampler)[source]

Plot diagnostic statistics from a dynesty run

The plotted quantities per iteration are:

  • nc: the number of likelihood calls

  • scale: the number of accepted MCMC steps if using bound="live" or bound="live-multi", otherwise, the scale applied to the MCMC steps

  • lifetime: the number of iterations a point stays in the live set

There is also a histogram of the lifetime compared with the theoretical distribution. To avoid edge effects, we discard the first 6 * nlive

Parameters:
sampler: dynesty.sampler.Sampler

The sampler object containing the run history.

Returns:
fig: matplotlib.pyplot.figure.Figure

Figure handle for the new plot

axs: matplotlib.pyplot.axes.Axes

Axes handles for the new plot