Linear models and empirical bayes methods for assessing differential expression in microarray experiments
- PMID: 16646809
- DOI: 10.2202/1544-6115.1027
Linear models and empirical bayes methods for assessing differential expression in microarray experiments
Abstract
The problem of identifying differentially expressed genes in designed microarray experiments is considered. Lonnstedt and Speed (2002) derived an expression for the posterior odds of differential expression in a replicated two-color experiment using a simple hierarchical parametric model. The purpose of this paper is to develop the hierarchical model of Lonnstedt and Speed (2002) into a practical approach for general microarray experiments with arbitrary numbers of treatments and RNA samples. The model is reset in the context of general linear models with arbitrary coefficients and contrasts of interest. The approach applies equally well to both single channel and two color microarray experiments. Consistent, closed form estimators are derived for the hyperparameters in the model. The estimators proposed have robust behavior even for small numbers of arrays and allow for incomplete data arising from spot filtering or spot quality weights. The posterior odds statistic is reformulated in terms of a moderated t-statistic in which posterior residual standard deviations are used in place of ordinary standard deviations. The empirical Bayes approach is equivalent to shrinkage of the estimated sample variances towards a pooled estimate, resulting in far more stable inference when the number of arrays is small. The use of moderated t-statistics has the advantage over the posterior odds that the number of hyperparameters which need to estimated is reduced; in particular, knowledge of the non-null prior for the fold changes are not required. The moderated t-statistic is shown to follow a t-distribution with augmented degrees of freedom. The moderated t inferential approach extends to accommodate tests of composite null hypotheses through the use of moderated F-statistics. The performance of the methods is demonstrated in a simulation study. Results are presented for two publicly available data sets.
Similar articles
-
An empirical bayesian method for differential expression studies using one-channel microarray data.Stat Appl Genet Mol Biol. 2003;2:Article8. doi: 10.2202/1544-6115.1024. Epub 2003 Nov 29. Stat Appl Genet Mol Biol. 2003. PMID: 16646786
-
Quality optimised analysis of general paired microarray experiments.Stat Appl Genet Mol Biol. 2006;5:Article10. doi: 10.2202/1544-6115.1209. Epub 2006 Apr 21. Stat Appl Genet Mol Biol. 2006. PMID: 16646864
-
Assessing differential gene expression with small sample sizes in oligonucleotide arrays using a mean-variance model.Biometrics. 2007 Mar;63(1):41-9. doi: 10.1111/j.1541-0420.2006.00675.x. Biometrics. 2007. PMID: 17447928
-
Accurate ranking of differentially expressed genes by a distribution-free shrinkage approach.Stat Appl Genet Mol Biol. 2007;6:Article9. doi: 10.2202/1544-6115.1252. Epub 2007 Feb 23. Stat Appl Genet Mol Biol. 2007. PMID: 17402924
-
Semi-parametric differential expression analysis via partial mixture estimation.Stat Appl Genet Mol Biol. 2008;7(1):Article15. doi: 10.2202/1544-6115.1333. Epub 2008 Apr 28. Stat Appl Genet Mol Biol. 2008. PMID: 18454730
Cited by
-
MicroRNA‑101 inhibits renal tubular epithelial‑to‑mesenchymal transition by targeting TGF‑β1 type I receptor.Int J Mol Med. 2021 Jun;47(6):119. doi: 10.3892/ijmm.2021.4952. Epub 2021 May 6. Int J Mol Med. 2021. PMID: 33955520 Free PMC article.
-
A multi-layer inference approach to reconstruct condition-specific genes and their regulation.Bioinformatics. 2013 Jun 15;29(12):1541-52. doi: 10.1093/bioinformatics/btt186. Epub 2013 Apr 22. Bioinformatics. 2013. PMID: 23610368 Free PMC article.
-
The PTIP-Associated Histone Methyltransferase Complex Prevents Stress-Induced Maladaptive Cardiac Remodeling.PLoS One. 2015 May 22;10(5):e0127839. doi: 10.1371/journal.pone.0127839. eCollection 2015. PLoS One. 2015. PMID: 26001054 Free PMC article.
-
Poly-glutamine expanded huntingtin dramatically alters the genome wide binding of HSF1.J Huntingtons Dis. 2012;1(1):33-45. J Huntingtons Dis. 2012. PMID: 23293686 Free PMC article.
-
Loss of cardiolipin leads to perturbation of mitochondrial and cellular iron homeostasis.J Biol Chem. 2013 Jan 18;288(3):1696-705. doi: 10.1074/jbc.M112.428938. Epub 2012 Nov 28. J Biol Chem. 2013. PMID: 23192348 Free PMC article.
LinkOut - more resources
Full Text Sources
Other Literature Sources