Journal
STATISTICAL SCIENCE
Volume 25, Issue 3, Pages 388-407Publisher
INST MATHEMATICAL STATISTICS-IMS
DOI: 10.1214/10-STS339
Keywords
EM algorithm; empirical Bayes; Laplace approximation; LEMMA; LIMMA; linear mixed models; local false discovery rate; microarray analysis; mixture model; two-groups model
Categories
Funding
- NSF [DMS 0805865]
- NIH [R01-GM083606-01]
Ask authors/readers for more resources
A two-groups mixed-effects model for the comparison of (normalized) microarray data from two treatment groups is considered. Most competing parametric methods that have appeared in the literature are obtained as special cases or by minor modification of the proposed model. Approximate maximum likelihood fitting is accomplished via a fast and scalable algorithm, which we call LEMMA (Laplace approximated EM Microarray Analysis). The posterior odds of treatment x gene interactions, derived from the model, involve shrinkage estimates of both the interactions and of the gene specific error variances. Genes are classified as being associated with treatment based on the posterior odds and the local false discovery rate (f.d.r.) with a fixed cutoff. Our model-based approach also allows one to declare the non-null status of a gene by controlling the false discovery rate (FDR). It is shown in a detailed simulation study that the approach outperforms well-known competitors. We also apply the proposed methodology to two previously analyzed microarray examples. Extensions of the proposed method to paired treatments and multiple treatments are also discussed.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available