However, as this is highly computationally demanding, we estimate it using sampling.
for every M = 1, ..., 20 and x i , we sample Y i 1000 times (with replacement), and estimate P (y ∈ Y i ) as the covered probability mass P D i {y : y ∈ Y i }.