simultaneously throughout a range of sparsity classes,
provided the level q for the false discovery rate is set to q < 0.5. Generalizations of this methodology to non-orthogonal
designs differ mainly in the generation of the p-values for testing Hj0
: βj = 0,