Rxivist logo

Rxivist combines preprints from bioRxiv with data from Twitter to help you find the papers being discussed in your field. Currently indexing 65,624 bioRxiv papers from 290,715 authors.

Nonparametric expression analysis using inferential replicate counts

By Anqi Zhu, Avi Srivastava, Joseph G. Ibrahim, Rob Patro, Michael I Love

Posted 25 Feb 2019
bioRxiv DOI: 10.1101/561084 (published DOI: 10.1093/nar/gkz622)

A primary challenge in the analysis of RNA-seq data is to identify differentially expressed genes or transcripts while controlling for technical biases. Ideally, a statistical testing procedure should incorporate the inherent uncertainty of the abundance estimates arising from the quantification step. Most popular methods for RNA-seq differential expression analysis fit a parametric model to the counts for each gene or transcript, and a subset of methods can incorporate uncertainty. Previous work has shown that nonparametric models for RNA-seq differential expression may have better control of the false discovery rate, and adapt well to new data types without requiring reformulation of a parametric model. Existing nonparametric models do not take into account inferential uncertainty, leading to an inflated false discovery rate, in particular at the transcript level. We propose a nonparametric model for differential expression analysis using inferential replicate counts, extending the existing SAMseq method to account for inferential uncertainty. We compare our method, Swish, with popular differential expression analysis methods. Swish has improved control of the false discovery rate, in particular for transcripts with high inferential uncertainty. We apply Swish to a single-cell RNA-seq dataset, assessing differential expression between sub-populations of cells, and compare its performance to the Wilcoxon test.

Download data

  • Downloaded 1,150 times
  • Download rankings, all-time:
    • Site-wide: 5,992 out of 65,624
    • In bioinformatics: 1,132 out of 6,492
  • Year to date:
    • Site-wide: 1,399 out of 65,624
  • Since beginning of last month:
    • Site-wide: 7,707 out of 65,624

Altmetric data

Downloads over time

Distribution of downloads per paper, site-wide

Sign up for the Rxivist weekly newsletter! (Click here for more details.)