Question: Can someone please explain in simple terms how DESeq2 works?
10
4.3 years ago by
BioStars1993100
BioStars1993100 wrote:

As simply as possible please, and if you don't mind giving a quick definition of any specialist terms you use.

deseq2 • 11k views
modified 2.0 years ago by Biostar ♦♦ 20 • written 4.3 years ago by BioStars1993100

19
4.3 years ago by
Michael Love1.8k
United States
Michael Love1.8k wrote:

We really tried to write the main text of the paper such that it would be understood by non-statisticians. that said, I'll try to do it in a few sentences. before further questions on details please at least *try* to read the paper :)

Let's say we want to compare counts between two groups. We build a model for the observed counts. This model has some parameters: (1) a normalization parameter, for differences in library size at least, or it can be extended by other software; (2) a variance parameter, called dispersion; (3) parameters representing the group differences. Fit (1) using the same method from the original DESeq. Fit (2) in two steps: first find the value of the parameter that makes the likelihood largest, which is called maximum likelihood estimation. Look at all the values from all of the genes and move these values towards a middle value. Bayes theorem guides the amount of movement for each gene: if the information for the gene is low, the value is moved more to the middle, if the information for the gene is high, the value is moved very little. Fit (3) using the same technique as used for (2). The values for (3) are a useful final product, as are sets of genes where the group differences are likely to be above a threshold (zero or otherwise). These sets are defined by their false discovery rate.

10
4.3 years ago by
geek_y9.4k
Barcelona/CRG/London/Imperial
geek_y9.4k wrote:

Look at posts by Simon Anders on Seqanswers. He explains in plain english about differential expression analysis in DESeq and also edgeR.

Few of the interesting discussions are here:

This is very nice link. I have been searching for it, just got it.

1
4.3 years ago by
karl.stamm3.4k
United States
karl.stamm3.4k wrote:

The publication is the most complete explanation you can get.

http://genomebiology.com/2014/15/12/550

In short: it estimates expression values, and calculates differential expression.

6

Not sure that the publication qualifies for "as simply as possible" :)