def METHOD_NAME(f, n_iters=100, seed=None):
    """
    Transforms a function that computes a sample statistic to a
    function that estimates the corresponding population statistic and
    uncertainty via bootstrap.
    In this version, the positional arguments to the function are
    assumed to be equal-length arrays, with the ith index of the jth
    array representing quantity j measured at observation i; the data
    are assumed correlated across quantities at a fixed observation,
    and uncorrelated across observations of a fixed quantity.
    Parameters
    ----------
    f : callable
        Function of one or more arrays returning a scalar. The
        positional arguments are assumed to be 1-d arrays of
        consistent length to be resampled at each iteration. Keyword
        arguments are passed through but are not resampled.
    n_iter : int
        Number of bootstrap iterations
    seed : int, optional
        Random seed
    Returns
    -------
    callable
        Function with the same input signature as `f`, returning a
        pair of scalars: (estimate, error)
    """
    return _bootstrap(samples_correlated, f, n_iters, seed)