A Bayesian discovery procedure

Michele Guindani, Peter Müller, Song Zhang

Research output: Contribution to journalArticle

29 Scopus citations

Abstract

We discuss a Bayesian discovery procedure for multiple-comparison problems. We show that, under a coherent decision theoretic framework, a loss function combining true positive and false positive counts leads to a decision rule that is based on a threshold of the posterior probability of the alternative. Under a semiparametric model for the data, we show that the Bayes rule can be approximated by the optimal discovery procedure, which was recently introduced by Storey. Improving the approximation leads us to a Bayesian discovery procedure, which exploits the multiple shrinkage in clusters that are implied by the assumed non-parametric model. We compare the Bayesian discovery procedure and the optimal discovery procedure estimates in a simple simulation study and in an assessment of differential gene expression based on microarray data from tumour samples. We extend the setting of the optimal discovery procedure by discussing modifications of the loss function that lead to different single-thresholding statistics. Finally, we provide an application of the previous arguments to dependent (spatial) data.

Original languageEnglish (US)
Pages (from-to)905-925
Number of pages21
JournalJournal of the Royal Statistical Society. Series B: Statistical Methodology
Volume71
Issue number5
DOIs
Publication statusPublished - Nov 2009

    Fingerprint

Keywords

  • Bayes optimal rule
  • False discovery rate
  • Loss function
  • Multiple comparison

ASJC Scopus subject areas

  • Statistics and Probability
  • Statistics, Probability and Uncertainty

Cite this