Search
Question: DESeq2 Which gene used to normalize?
0
14 months ago by
michi.s0
michi.s0 wrote:

Is it possible to get a function which gives the name of the gene used to normalize?

modified 14 months ago by Steve Lianoglou12k • written 14 months ago by michi.s0
0
14 months ago by
Michael Love20k
United States
Michael Love20k wrote:

hi,

Take a look at the original DESeq paper. DESeq2 uses the median ratio method described here for normalization. The default normalization is the estimation of size factors presented in Equation 5. The default size factor estimate uses all of the genes that have positive counts for all samples.

https://genomebiology.biomedcentral.com/articles/10.1186/gb-2010-11-10-r106

Hi Michael,

thank you for the fast reply. I now the method, but my point is, if I take the median of ratios to the geometric mean, then there is in fact for each sample one gene I normalize to (or two genes, if the total number of genes is even). This information could be really helpful for explaining biologists what happened to their data...

0
14 months ago by
Denali
Steve Lianoglou12k wrote:

Unless you specify a set of control genes to use, every gene which has all-non-zero-observations is used.

See the help in ?estimateSizeFactors and ?estimateSizeFactorsForMatrix for more information.