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)