def METHOD_NAME(Y, X, m, num_resamples, conf_level, y_resamples):
    s = np.zeros(num_resamples)
    N = len(Y)
    r = np.random.randint(N, size=(num_resamples, y_resamples))
    for i in range(num_resamples):
        r_i = r[i, :]
        s[i] = sobol_first(Y[r_i], X[r_i], m)
    return norm.ppf(0.5 + conf_level / 2) * s.std(ddof=1)