An empirical Bayes mixture method for effect size and false discovery rate estimation

Abstract

Many statistical problems involve data from thousands of parallel cases. Each case has some associated effect size, and most cases will have no effect. It is often important to estimate the effect size and the local or tail-area false discovery rate for each case. Most current methods do this separately, and most are designed for normal data. This paper uses an empirical Bayes mixture model approach to estimate both quantities together for exponential family data. The proposed method yields simple, interpretable models that can still be used nonparametrically. It can also estimate an empirical null and incorporate it fully into the model. The method outperforms existing effect size and false discovery rate estimation procedures in normal data simulations; it nearly acheives the Bayes error for effect size estimation. The method is implemented in an R package (mixfdr), freely available from CRAN.

Published in at http://dx.doi.org/10.1214/09-AOAS276 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

References in corpus (7)

• Needles and straw in haystacks: Empirical Bayes estimates of possibly sparse sequences
• Microarrays, Empirical Bayes and the Two-Groups Model
• Estimating the proportion of false null hypotheses among a large number of independently tested hypotheses
• In-season prediction of batting averages: A field test of empirical Bayes and Bayes methodologies
• Estimation and confidence sets for sparse normal mixtures
• On optimality of Bayesian testimation in the normal means problem
• Frequentist optimality of Bayesian wavelet shrinkage rules for Gaussian and non-Gaussian noise

Cited by in corpus (10)

• Replicability analysis for genome-wide association studies
• A nonparametric empirical Bayes framework for large-scale multiple testing
• Optimal shrinkage estimation of mean parameters in family of distributions with quadratic variance
• Optimal control of false discovery criteria in the two-group model
• Empirical Bayes estimation of posterior probabilities of enrichment
• Bayesian latent hierarchical model for transcriptomic meta-analysis to detect biomarkers with clustered meta-patterns of differential expression signals
• Testing the disjunction hypothesis using Voronoi diagrams with applications to genetics
• On nonparametric estimation of a mixing density via the predictive recursion algorithm
• Feature-specific inference for penalized regression using local false discovery rates
• Estimating the Mean and Variance of a High-dimensional Normal Distribution Using a Mixture Prior