Bayesian optimal discovery procedure for simultaneous significance testing.
BACKGROUND
In high throughput screening, such as differential gene expression screening, drug sensitivity screening, and genome-wide RNAi screening, tens of thousands of tests need to be conducted simultaneously. However, the number of replicate measurements per test is extremely small, rarely exceeding 3. Several current approaches demonstrate that test statistics with shrinking variance estimates have more power over the traditional t statistic.


RESULTS
We propose a Bayesian hierarchical model to incorporate the shrinkage concept by introducing a mixture structure on variance components. The estimates from the Bayesian model are utilized in the optimal discovery procedure (ODP) proposed by Storey in 2007, which was shown to have optimal performance in multiple significance tests. We compared the performance of the Bayesian ODP with several competing test statistics.


CONCLUSIONS
We have conducted simulation studies with 2 to 6 replicates per gene. We have also included test results from two real datasets. The Bayesian ODP outperforms the other methods in our study, including the original ODP. The advantage of the Bayesian ODP becomes more significant when there are few replicates per test. The improvement over the original ODP is based on the fact that Bayesian model borrows strength across genes in estimating unknown parameters. The proposed approach is efficient in computation due to the conjugate structure of the Bayesian model. The R code (see Additional file 1) to calculate the Bayesian ODP is provided.
Cao
Jing
J
Department of Statistical Science, Southern Methodist University, Dallas, Texas, USA. jcao@smu.edu
Xie
Xian-Jin
XJ
Zhang
Song
S
Whitehurst
Angelique
A
White
Michael A
MA
D000465:Algorithms; D001499:Bayes Theorem; D019295:Computational Biology / Q000379:methods*; D003198:Computer Simulation; D016208:Databases, Factual; D020869:Gene Expression Profiling / Q000379:methods; D006801:Humans; D008957:Models, Genetic; D015233:Models, Statistical; D008962:Models, Theoretical; D010363:Pattern Recognition, Automated / Q000379:methods; D011381:Programming Languages; D034622:RNA Interference; D012372:ROC Curve; D015203:Reproducibility of Results