The Annals of Applied Statistics
- Ann. Appl. Stat.
- Volume 9, Number 2 (2015), 901-925.
hmmSeq: A hidden Markov model for detecting differentially expressed genes from RNA-seq data
We introduce hmmSeq, a model-based hierarchical Bayesian technique for detecting differentially expressed genes from RNA-seq data. Our novel hmmSeq methodology uses hidden Markov models to account for potential co-expression of neighboring genes. In addition, hmmSeq employs an integrated approach to studies with technical or biological replicates, automatically adjusting for any extra-Poisson variability. Moreover, for cases when paired data are available, hmmSeq includes a paired structure between treatments that incoporates subject-specific effects. To perform parameter estimation for the hmmSeq model, we develop an efficient Markov chain Monte Carlo algorithm. Further, we develop a procedure for detection of differentially expressed genes that automatically controls false discovery rate. A simulation study shows that the hmmSeq methodology performs better than competitors in terms of receiver operating characteristic curves. Finally, the analyses of three publicly available RNA-seq data sets demonstrate the power and flexibility of the hmmSeq methodology. An R package implementing the hmmSeq framework will be submitted to CRAN upon publication of the manuscript.
Ann. Appl. Stat., Volume 9, Number 2 (2015), 901-925.
Received: February 2014
Revised: January 2015
First available in Project Euclid: 20 July 2015
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
Cui, Shiqi; Guha, Subharup; Ferreira, Marco A. R.; Tegge, Allison N. hmmSeq: A hidden Markov model for detecting differentially expressed genes from RNA-seq data. Ann. Appl. Stat. 9 (2015), no. 2, 901--925. doi:10.1214/15-AOAS815. https://projecteuclid.org/euclid.aoas/1437397117
- Supplement to “hmmSeq: A hidden Markov model for detecting differentially expressed genes from RNA-seq data.”. The R code for hmmSeq.