A linear model involving a mixture distribution is considered for the comparison of normalized microarray data from two treatment groups. Model fitting using an empirical Bayes approach has been shown to be both accurate and numerically stable. The posterior odds of treatment/gene interactions derived from the model involve shrinkage estimates of both the interactions and the gene-specific error variances, leading to powerful inference. We show that the same model can easily be fit under a fully Bayesian framework, allowing increased flexibility in terms of prior distributional assumptions and posterior inference.
Baldi P and Long AD (2001) A Bayesian framework for the analysis of microarray expression data: regularized t-test and statistical inferences of gene changes . Bioinformatics, 17, 509–19 .
2.
Bar H , Booth J , Schifano E and Wells M (2010) Laplace approximated EM microarray analysis: an empirical Bayes approach for comparative microarray experiments . Statistical Science. http://www.stat.cornell.edu/lemma/
3.
Bar H and Schifano E (2009) lemma: Laplace approximated EM microarray analysis. R package version 1.2-1. http://www.stat.cornell.edu/lemma
4.
Brot P , Richardson S and Radvanyi F (2002) Bayesian hierarchical model for identifying changes in gene expression from microarray experiments . Journal of Computational Biology, 9, 671–83 .
5.
Callow MJ , Dudoit S , Gong EL , Speed TP and Rubin EM (2000) Microarray expression profiling identifies genes with altered expression in HDL-deficient mice . Genome Research, 10, 2022–29 .
6.
Cui X , Hwang JTG , Qui J , Blades NJ and Churchill GA (2005) Improved statistical tests for differential gene expression by shrinking variance components . Biostatistics, 6, 59–75 .
7.
Efron B (2005) Local false discovery rates. http://www-stat.stanford.edu/brad/papers/False.pdf
8.
Efron B (2008) Microarrays, empirical bayes and the two groups model . Statistical Science, 23, 1–22 .
9.
Efron B , Turnbull BB and Narasimhan B (2008) locfdr: Computes local false discovery rates. R package version 1.1-6.
10.
Gelman A , Carlin JB , Stern HS and Rubin DB (2004) Bayesian data analysis, 2nd edn. North Carolina: Chapman and Hall .
11.
Kendziorski CM , Newton MA , Lan H and Gould MN (2003) On parametric empirical Bayes methods for comparing multiple groups using replicated gene expression profiles . Statistics in Medicine, 22, 3899–914 .
12.
Liu P (2006) Sample size calculation and empirical Bayes tests for microarray data. Ph.D. thesis, Cornell University, Ithaca, NY .
13.
Lonnstedt I and Britton T (2005) Hierarchical Bayes models for cDNA microarray gene expression . Biostatistics, 6, 279–91 .
14.
Lonnstedt I and Speed T (2002) Replicated microarray data . Statistica Sinica, 12, 31–46 .
15.
Newton MA , Kendziorski CM , Richmond CS , Blattner FR and Tsui KW (2001) On differential variability of expression ratios: improving statistical inference about gene expression changes from microarray data . Computational Biology, 8, 37–52 .
16.
Smyth GK (2004) Linear models for empirical bayes methods for assessing differential expression in microarray experiments . Statistical Applications in Genetics and Molecular Biology, 3, Article 2.
17.
Smyth GK (2005) Limma: linear models for microarray data. New York: Springer , pp. 397–420.
18.
Spiegelhalter D , Thomas A , Best N and Lunn D (2003) WinBUGS user manual Version 1.4. http://www.mrc-bsu.cam.ac.uk./bugs/winbugs/manual14.pdf
19.
Wright GW and Simon RM (2003) A random variance model for detection of differential gene expression in small microarray experiments . Bioinformatics, 19, 2448–55 .
