Shrinkage estimation of dispersion in Negative Binomial models for RNA-seq experiments with small sample size
Open Access
- 14 April 2013
- journal article
- research article
- Published by Oxford University Press (OUP) in Bioinformatics
- Vol. 29 (10) , 1275-1282
Abstract
Motivation: RNA-seq experiments produce digital counts of reads that are affected by both biological and technical variation. To distinguish the systematic changes in expression between conditions from noise, the counts are frequently modeled by the Negative Binomial distribution. However, in experiments with small sample size, the per-gene estimates of the dispersion parameter are unreliable. Method: We propose a simple and effective approach for estimating the dispersions. First, we obtain the initial estimates for each gene using the method of moments. Second, the estimates are regularized, i.e. shrunk towards a common value that minimizes the average squared difference between the initial estimates and the shrinkage estimates. The approach does not require extra modeling assumptions, is easy to compute and is compatible with the exact test of differential expression. Results: We evaluated the proposed approach using 10 simulated and experimental datasets and compared its performance with that of currently popular packages edgeR, DESeq, baySeq, BBSeq and SAMseq. For these datasets, sSeq performed favorably for experiments with small sample size in sensitivity, specificity and computational time. Availability:http://www.stat.purdue.edu/∼ovitek/Software.html and Bioconductor. Contact:ovitek@purdue.edu Supplementary information: Supplementary data are available at Bioinformatics online.Keywords
This publication has 46 references indexed in Scilit:
- Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variationNucleic Acids Research, 2012
- Normalization, testing, and false discovery rate estimation for RNA-sequencing dataBiostatistics, 2011
- A powerful and flexible approach to the analysis of RNA sequence count dataBioinformatics, 2011
- Computational methods for transcriptome annotation and quantification using RNA-seqNature Methods, 2011
- Conservation of an RNA regulatory map between Drosophila and mammalsGenome Research, 2010
- edgeR: a Bioconductor package for differential expression analysis of digital gene expression dataBioinformatics, 2009
- DEGseq: an R package for identifying differentially expressed genes from RNA-seq dataBioinformatics, 2009
- Computation for ChIP-seq and RNA-seq studiesNature Methods, 2009
- RNA-Seq: a revolutionary tool for transcriptomicsNature Reviews Genetics, 2009
- The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurementsNature Biotechnology, 2006