docs: def estimate_tau_int(chains, **kwargs):
      """
      Estimate the integrated auto-correlation time for all parameters in a chain.
      """
      return autocorr.integrated_time(np.mean(chains, axis=0), **kwargs)