Optimal multiple testing under a Gaussian prior on the effect sizes

Abstract

We develop a new method for large-scale frequentist multiple testing with Bayesian prior information. We find optimal p -value weights that maximize the average power of the weighted Bonferroni method. Due to the nonconvexity of the optimization problem, previous methods that account for uncertain prior information are suitable for only a small number of tests. For a Gaussian prior on the effect sizes, we give an efficient algorithm that is guaranteed to find the optimal weights nearly exactly. Our method can discover new loci in genome-wide association studies and compares favourably to competitors. An open-source implementation is available.