In our experiment we used nd=30, n′=48, n=48, T=2, θ=0.05, and α=0.05 as the parameter values, unless otherwise stated. Furthermore, peptideIndex was set to represent a random subset of nd peptides, and phosphorylated was set to be a random binary vector of the same length as peptideIndex, i.e., nd.