compute_set_pvalues function

Compute BEAMR p-values for sets