Bayesian Analysis

Bayesian analysis of comparative microarray experiments by model averaging

Marco F. Ramoni, Paola Sebastiani, and Hui Xie

Full-text: Open access


A major challenge to the statistical analysis of microarray data is the small number of samples — limited by both cost and sample availability — compared to the large number of genes, now soaring into the tens of thousands per experiment. This situation is made even more difficult by the complex nature of the empirical distributions of gene expression measurements and the necessity to limit the number of false detections due to multiple comparisons. This paper introduces a novel Bayesian method for the analysis of comparative experiments performed with oligonucleotide microarrays. Our method models gene expression data by log-normal and gamma distributions with hierarchical prior distributions on the parameters of interest, and uses model averaging to compute the posterior probability of differential expression. An initial approximate Bayesian analysis is used to identify genes that have a large probability of differential expression, and this list of candidate genes is further refined by using stochastic computations. We assess the performance of this method using real data sets and show that it has an almost negligible false positive rate in small sample experiments that leads to a better detection performance.

Article information

Bayesian Anal., Volume 1, Number 4 (2006), 707-732.

First available in Project Euclid: 22 June 2012

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Microarrays gene expression gamma distribution log-normal distribution model averaging true and false positive rates false discovery rate


Sebastiani, Paola; Xie, Hui; Ramoni, Marco F. Bayesian analysis of comparative microarray experiments by model averaging. Bayesian Anal. 1 (2006), no. 4, 707--732. doi:10.1214/06-BA123.

Export citation