Robust hyperparameter estimation protects against hypervariable genes and improves power to detect differential expression

Annals of Applied Statistics

Volumn 10, Issue 2, 2016, Pages 946-963

  Phipson, Belinda ^a   Lee, Stanley ^b   Majewski, Ian J ^b   Alexander, Warren S ^b   Smyth, Gordon K ^b  

^a MURDOCH CHILDREN'S RESEARCH INSTITUTE   (Australia)

^b UNIVERSITY OF MELBOURNE   (Australia)

Author keywords

Empirical Bayes; Gene expression; Microarrays; Outliers; Robustness

Indexed keywords

EID: 84979917177     PISSN: 19326157     EISSN: 19417330     Source Type: Journal    
DOI: 10.1214/16-AOAS920     Document Type: Article

References (42)

1. BALDI, P. and LONG, A. D. (2001). A Bayesian framework for the analysis of microarray expression data: Regularized t-test and statistical inferences of gene changes. Bioinformatics 17 509–519.
2. BARLOW, R. E., BARTHOLOMEW, D. J., BREMNER, J. M. and BRUNK, H. D. (1972). Statistical Inference Under Order Restrictions. Wiley, London. MR0326887
3. BENJAMINI, Y. and HOCHBERG, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B. Stat. Methodol. 57 289–300. MR1325392
4. BERGER, J. O. (1984). The robust Bayesian viewpoint. In Robustness of Bayesian Analyses (J. Kadane, ed.). Stud. Bayesian Econometrics 4 63–144. North-Holland, Amsterdam. With com-ments and with a reply by the author. MR0785367
5. BERGER, J. O. (1990). Robust Bayesian analysis: Sensitivity to the prior. J. Statist. Plann. Inference 25 303–328. MR1064429
6. BOLSTAD, B. M., IRIZARRY, R. A., ÅSTRand, M. and SPEED, T. P. (2003). A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19 185–193.
7. BRENT, R. P. (1973). Algorithms for Minimization Without Derivatives. Prentice-Hall, Englewood Cliffs, NJ.
8. CASELLA, G. (1985). An introduction to empirical Bayes data analysis. Amer. Statist. 39 83–87. MR0789118
9. CHEN, Y., LUN, A. T. L. and SMYTH, G. K. (2014). Differential expression analysis of complex RNA-seq experiments using edgeR. In Statistical Analysis of Next Generation Sequence Data (S. Datta and D. S. Nettleton, eds.) 51–74. Springer, New York.
10. CLEVELand, W. S. (1979). Robust locally weighted regression and smoothing scatterplots. J. Amer. Statist. Assoc. 74 829–836. MR0556476
11. EFRON, B. and MORRIS, C. (1972). Limiting the risk of Bayes and empirical Bayes estimators. II. The empirical Bayes case. J. Amer. Statist. Assoc. 67 130–139. MR0323015
12. EFRON, B. and MORRIS, C. (1973). Stein’s estimation rule and its competitors—An empirical Bayes approach. J. Amer. Statist. Assoc. 68 117–130. MR0388597
13. GAVER, D. P. and O'MUIRCHEARTAIGH, I. G. (1987). Robust empirical Bayes analyses of event rates. Technometrics 29 1–15. MR0876882
14. GOOD-JACOBSON, K. L., CHEN, Y., VOSS, A. K., SMYTH, G. K., THOMAS, T. and TARLIN-TON, D. (2014). Regulation of germinal center responses and B-cell memory by the chromatin modifier MOZ. Proc. Natl. Acad. Sci. USA 111 9585–9590.
15. GOTTARDO, R., RAFTERY, A. E., YEUNG, K. Y. and BUMGARNER, R. E. (2006). Bayesian robust inference for differential gene expression in microarrays with multiple samples. Biometrics 62 10–18. MR2226551
16. INSUA, D. R. and RUGGERI, F., eds. (2000). Robust Bayesian Analysis. Lecture Notes in Statistics 152. Springer, New York. MR1795206
17. JEANMOUGIN, M., DE REYNIES, A., MARISA, L., PACCARD, C., NUEL, G. and GUEDJ, M. (2010). Should we abandon the t-test in the analysis of gene expression microarray data: A comparison of variance modeling strategies. PLoS ONE 5 e12336.
18. JI, H. and LIU, X. S. (2010). Analyzing ’omics data using hierarchical models. Nature Biotechnology 28 337.
19. KOOPERBERG, C., ARAGAKI, A., STRand, A. D. and OLSON, J. M. (2005). Significance testing for small microarray experiments. Stat. Med. 24 2281–2298. MR2151706
20. LAW, C. W., CHEN, Y., SHI, W. and SMYTH, G. K. (2014). Voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biology 15 R29.
21. LIAO, J. G., MCMURRY, T. and BERG, A. (2014). Prior robust empirical Bayes inference for large-scale data by conditioning on rank with application to microarray data. Biostatistics 15 60–73.
22. LUN, A. T. L., CHEN, Y. and SMYTH, G. K. (2016). It’s DE-licious: A recipe for differential expression analyses of RNA-seq experiments using quasi-likelihood methods in edgeR. Methods in Molecular Biology 1418 391–416.
23. LUN, A. T. L. and SMYTH, G. K. (2015a). diffHic: A bioconductor package to detect differential genomic interactions in Hi-C data. BMC Bioinformatics 16 258.
24. LUN, A. T. L. and SMYTH, G. K. (2015b). From reads to regions: A Bioconductor workflow to detect differential binding in ChIP-seq data. F1000Research 4 1080.
25. LUN, A. T. L. and SMYTH, G. K. (2016). csaw: A bioconductor package for differential binding analysis of ChIP-seq data using sliding windows. Nucleic Acids Res. 44 e45.
26. MAJEWSKI, I. J., BLEWITT, M. E., DE GRAAF, C. A., MCMANUS, E. J., BAHLO, M., HILTON, A. A., HYLand, C. D., SMYTH, G. K., CORBIN, J. E., METCALF, D. et al. (2008). Polycomb repressive complex 2 (PRC2) restricts hematopoietic stem cell activity. PLOS Biology 6 e93.
27. MAJEWSKI, I. J., RITCHIE, M. E., PHIPSON, B., CORBIN, J., PAKUSCH, M., EBERT, A., BUS-SLINGER, M., KOSEKI, H., HU, Y., SMYTH, G. K. et al. (2010). Opposing roles of polycomb repressive complexes in hematopoietic stem and progenitor cells. Blood 116 731–739.
28. MCCARTHY, D. J. and SMYTH, G. K. (2009). Testing significance relative to a fold-change thresh-old is a TREAT. Bioinformatics 25 765–771.
29. MORRIS, C. N. (1983). Parametric empirical Bayes inference: Theory and applications. J. Amer. Statist. Assoc. 78 47–65. With discussion. MR0696849
30. MURIE, C., WOODY, O., LEE, A. Y. and NADON, R. (2009). Comparison of small n statistical tests of differential expression applied to microarrays. BMC Bioinformatics 10 45.
31. PHIPSON, B., LEE, S., MAJEWSKI, I. J., ALEXandER, W. S. and SMYTH, G. K. (2016). Supplement to “Robust hyperparameter estimation protects against hypervariable genes and improves power to detect differential expression.” DOI:10.1214/16-AOAS920SUPP.
32. PICKRELL, J. K., MARIONI, J. C., PAI, A. A., DEGNER, J. F., ENGELHARDT, B. E., NKADORI, E., VEYRIERAS, J.-B., STEPHENS, M., GILAD, Y. and PRITCHARD, J. K. (2010a). Understanding mechanisms underlying human gene expression variation with RNA sequencing. Nature 464 768–772.
33. PICKRELL, J. K., PAI, A. A., GILAD, Y. and PRITCHARD, J. K. (2010b). Noisy splicing drives MRNA isoform diversity in human cells. PLoS Genet. 6 e1001236.
34. RITCHIE, M. E., PHIPSON, B., WU, D., HU, Y., LAW, C. W., SHI, W. and SMYTH, G. K. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nu-cleic Acids Res. 43 e47.
35. ROBINSON, M. D., MCCARTHY, D. J. and SMYTH, G. K. (2010). edgeR: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26 139–140.
36. SARTOR, M. A., TOMLINSON, C. R., WESSELKAMPER, S. C., SIVAGANESAN, S., LEIKAUF, G. D. and MEDVEDOVIC, M. (2006). Intensity-based hierarchical Bayes method improves testing for differentially expressed genes in microarray experiments. BMC Bioinformatics 7 538.
37. SHEIKH, B. N., DOWNER, N. L., PHIPSON, B., VANYAI, H. K., KUEH, A. J., MCCARTHY, D. J., SMYTH, G. K., THOMAS, T. and VOSS, A. K. (2015). MOZ and BMI1 play opposing roles during Hox gene activation in ES cells and in body segment identity specification in vivo. Proc. Natl. Acad. Sci. USA 112 5437–5442.
38. SHI, W., OSHLACK, A. and SMYTH, G. K. (2010). Optimizing the noise versus bias trade-off for Illumina whole genome expression BeadChips. Nucleic Acids Res. 38 e204.
39. SMYTH, G. K. (2004). Linear models and empirical Bayes methods for assessing differential ex-pression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 3 Article 3. MR2101454
40. TUKEY, J. W. (1962). The future of data analysis. Ann. Math. Stat. 33 1–67. MR0133937
41. WRIGHT, G. W. and SIMON, R. M. (2003). A random variance model for detection of differential gene expression in small microarray experiments. Bioinformatics 19 2448–2455.
42. ZHOU, X., LINDSAY, H. and ROBINSON, M. D. (2014). Robustly detecting differential expression in RNA sequencing data using observation weights. Nucleic Acids Res. 42 e91.