@classmethod
def probs(cls, x, y, rope=0, *, nsamples=50000, **kwargs):
    """
        Compute and return probabilities

        Args:
            x (np.array): a vector of scores for the first model
            y (np.array): a vector of scores for the second model
            rope (float): the width of the region of practical equivalence (default: 0)
            nsamples (int): the number of samples (default: 50000)

        Returns:
            `(p_left, p_rope, p_right)` if `rope > 0`;
            otherwise `(p_left, p_right)`.
        """
    return cls(x, y, rope, nsamples=nsamples, **kwargs).probs(rope > 0)
