def preserve_subset_statistic()

in statistickway_threshold.py [0:0]


def preserve_subset_statistic(queries):
    @jit
    def compute_statistic(D):
        np.concatenate(
            [(1 - np.prod(1 - D[:, q], 2)).sum(0) for q in np.array_split(queries, 10)]
        )
        return (
            np.concatenate(
                [
                    (1 - np.prod(1 - D[:, q], 2)).sum(0)
                    for q in np.array_split(queries, 10)
                ]
            )
            / D.shape[0]
        )

    return compute_statistic