English

Weighted mining of massive collections of $p$-values by convex optimization

Methodology 2017-06-20 v2 Applications

Abstract

Researchers in data-rich disciplines---think of computational genomics and observational cosmology---often wish to mine large bodies of pp-values looking for significant effects, while controlling the false discovery rate or family-wise error rate. Increasingly, researchers also wish to prioritize certain hypotheses, for example those thought to have larger effect sizes, by upweighting, and to impose constraints on the underlying mining, such as monotonicity along a certain sequence. We introduce Princessp, a principled method for performing weighted multiple testing by constrained convex optimization. Our method elegantly allows one to prioritize certain hypotheses through upweighting and to discount others through downweighting, while constraining the underlying weights involved in the mining process. When the pp-values derive from monotone likelihood ratio families like the Gaussian means model, the new method allows exact solution of an important optimal weighting problem previously thought to be nonconvex and computationally infeasible. Our method scales to massive dataset sizes. We illustrate the applications of Princessp on a series of standard genomics datasets and offer comparisons with several previous `standard' methods. Princessp offers both ease of operation and the ability to scale to extremely large problem sizes. The method is available as open-source software from http://github.com/dobriban/pvalue_weighting_matlab .

Keywords

Cite

@article{arxiv.1603.05334,
  title  = {Weighted mining of massive collections of $p$-values by convex optimization},
  author = {Edgar Dobriban},
  journal= {arXiv preprint arXiv:1603.05334},
  year   = {2017}
}

Comments

This is an entirely rewritten version of the paper. The title of the paper, the name of the method, and the introduction have been changed, with the goal of making the paper more accessible and appealing to practitioners. New sections on monotone likelihood ratio families and two-sided tests have been added, which expand the scope of the method