bioRxivBioinformaticsPreprint
A revisit of RSEM generative model and its EM algorithm for quantifying transcript abundances.
Abstract
RSEM has been mainly known for its accuracy in transcript abundance quantification. However, its quantification time is extremely high compared to that of recent quantification tools. In this paper, we revised the RSEMs EM algorithm. In particular, we derived accurate M-step updates to eliminate incorrect heuristic updates in RSEM. We also implement some optimizations that reduce the quantification time about a hundred times while still have better accuracy compared to RSEM. In particular, we noticed that different parameters have different convergence rates, therefore we identified and removed early converged parameters to significantly reduce the model complexity in further iterations, and we also use SQUAREM method to further speed up the convergence rate. We implemented these revisions in a packaged named Hera-EM, with source code available at: https://github.com/bioturing/hera/tree/master/hera-EM
- Citations
- 2
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Simultaneous estimation of transcript abundances and transcript specific fragment distributions of RNA-Seq data with the Mix2 model
Andreas Tuerk, Gregor Wiktorin·bioRxiv·Jun 4, 2014
- RNA-seq transcript quantification from reduced-representation data in recount2
Jack M. Fu, Kai Kammers, … Margaret A. Taub·bioRxiv·May 25, 2018
- Salmon provides accurate, fast, and bias-aware transcript expression estimates using dual-phase inference
Rob Patro, Geet Duggal, … Carl Kingsford·Nature MethodsPeer-reviewed·Mar 6, 2017
- Modeling of RNA-seq fragment sequence bias reduces systematic errors in transcript abundance estimation
Michael I Love, John B Hogenesch, Rafael A Irizarry·bioRxiv·Aug 29, 2015
- The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads
Yang Liao, Gordon K. Smyth, Wei Shi·Nucleic Acids ResearchPeer-reviewed·Feb 17, 2019
- Scikit-ribo: Accurate estimation and robust modeling of translation dynamics at codon resolution
Han Fang, Yi-Fei Huang, … Michael C. Schatz·Cell SystemsPeer-reviewed·Jan 17, 2018