Question: Deseq2: RNA-seq and Negative binomial distribution
2
gravatar for tony_cybercloud
3.9 years ago by
United States
tony_cybercloud20 wrote:

Hello,

I would like to know. Why is a Negative Binomial distribution assumed for RNA-seq data?

How is this distribution assumption used in Deseq2?

Related to the questions above: what's the difference with microarray data and limma?

 

Thanks a lot in advance

Any suggestions would be much appreciated.

Tony.

rnaseq deseq2 • 8.5k views
ADD COMMENTlink modified 3.9 years ago by Michael Love25k • written 3.9 years ago by tony_cybercloud20
Answer: Deseq2: RNA-seq and Negative binomial distribution
7
gravatar for Michael Love
3.9 years ago by
Michael Love25k
United States
Michael Love25k wrote:

A distributional assumption is needed because we want to estimate the probability of extreme events (large fold change just appearing by chance) from limited replicates. The negative binomial (a.k.a. Gamma-Poisson) is a good choice for RNA-seq experiments because

  1. the observed data at gene level is inherently counts or estimated counts of fragments for each feature and
  2. the spread of values among biological replicates is more than given by a simpler, one parameter distribution, the Poisson; and it seems to be captured by the NB sufficiently well

(The NB may not necessarily be a sufficient for single cell experiments, because here it's been observed that more parameters may be needed to capture the probability of "drop out" of fragments.) For more details, check out the paper for DESeq2, as well as the DESeq and edgeR papers. Simon gave a lecture at the Bioconductor summer course in Brixen which also discusses distributions: http://www.bioconductor.org/help/course-materials/2015/CSAMA2015/lect/L05-deseq2-anders.pdf For your question on limma, I recommend you read the voom paper which helps explain what extra is needed to use linear models on count data.

ADD COMMENTlink modified 3.9 years ago by Wolfgang Huber13k • written 3.9 years ago by Michael Love25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 197 users visited in the last hour