This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DESEQ2 size factors change when number of samples change in the design

For example, when you analyse Untreated condition vs Day1Treatment condition with only these samples in your design you will get some size factors.

If you have also set up a Day6Treatment condition in your design (but you don't use it yet), and want to to do the same comparison Untreated vs Day1Treatment, the size factors will change. (taking into account Day6)

I thought size factors was only relative to library size of its sample , so why do they change dependent on the design file even if you don't use a set of samples.

I was doing analyse Day1Treatment vs Untreated & Day6Treatment vs Untreated with two separate design files. But now I am wondering if it's better to have one design with all the samples, and do the two comparison to get same sizefactors because at the end of the day you finish with different differential genes detected.

deseq2
• 3,384 views • link
updated by erwan.scaon erwan.scaon • written by ZheFrench ZheFrench

2 answers

SizeFactors are calculated based on all the samples in your dds object. So if you have 1 large experiment its better to put everything together and then perform comparisons for example using contrasts, rather than making new dds objects with subsets of your original dataset.

• 1 views • link
written by grant.hovhannisyan grant.hovhannisyan

Agreed. The idea is that you estimate a size factor for each column that best scales the datasets based on a large set of genes that do not change upon conditions. Given that you do not have samples with extreme global changes, it is probably the best to have as many samples in the matrix as possible. This probably produces more robust size factors than with only two or three samples.

• 1 views • link
written by ATpoint ATpoint

Ok I got it , but do you know how is it computed ?

• 1 views • link
written by ZheFrench ZheFrench

Yes, it is actually pretty simple but powerful from the concept. Check out StatQuest for a nice explanation.

• 1 views • link
written by ATpoint ATpoint

So if you have 1 large experiment its better to put everything together and then perform comparisons

It is still true if said large experiment was done on multiple Illumina runs (compute size factor for all samples on all runs) ? Or shall we compute size factor for all samples per Illumina run ?

• 0 views • link
written by erwan.scaon erwan.scaon

You should put everything together and also ideally add batch effect to your formula design.

• 0 views • link
written by grant.hovhannisyan grant.hovhannisyan

Log in to answer this question.

More posts like this