Question: Deseq Sizefactors Help
0
gravatar for Sara
7.6 years ago by
Sara130
Sara130 wrote:

Hello,

I try to do the diff analysis by DESeq on two samples coming from two different condition without replicate

I run the script of DESeq and when I give the library size of my sample, the size factor is 1 for both conditions

libsizes <- c( Cond1 =80040653, Cond2 =360265740)
cds=estimateSizeFactors(cds)
sizeFactors(cds)

    Cond1 Cond2 
       1    1

I don't understand this estimation?

when I don't give the libsizes, it become

 sizeFactors(cds)
     Cond1 Cond2 
1.7320508 0.5773503

Could you please explain for me this issue?

Thanks in advance Sara

deseq • 12k views
ADD COMMENTlink modified 7.6 years ago by Damian Kao15k • written 7.6 years ago by Sara130

Probably you assigned sizes first and then did the estimate, which may just scale your assigned sizes back to 1. I think the right way is to estimate the factor directly from the 'cds' object without assigning any libsizes.

ADD REPLYlink written 7.6 years ago by Vitis2.1k

Hi Jeremy thanks for your comment could u please tell me what did this sizeFactors function? is it for normalization of data (scale normalization)?

ADD REPLYlink written 7.6 years ago by Sara130
9
gravatar for Damian Kao
7.4 years ago by
Damian Kao15k
USA
Damian Kao15k wrote:

The author of DESeq wrote a post on SeqAnswers a while back about how the size factors thing work in DESeq. I'll try to find the post, but basically it normalizes the datasets by:

-Take the geometric mean of each condition for a gene and use that as the reference expression data set.

-For each condition, get a list of quotients of each gene expression value to its reference expression.

-The median of each condition quotient list is the normalization factor for that data set.

edit*

Here is the post: http://seqanswers.com/forums/showpost.php?p=16468&postcount=13

ADD COMMENTlink modified 7.4 years ago • written 7.4 years ago by Damian Kao15k
0
gravatar for Jeremy Leipzig
7.6 years ago by
Philadelphia, PA
Jeremy Leipzig18k wrote:

I don't see where in the documentation it says libsizes is a magic word.

Either use the estimateSizeFactors to use your count data as the source of the estimate, or set it manually using sizeFactors.

Now, why are there so many Sara's?

alt text

ADD COMMENTlink written 7.6 years ago by Jeremy Leipzig18k
3

Because Sara is a very popular name amongst bioinformatician LOL

ADD REPLYlink written 7.5 years ago by Pasta1.3k
1

could be some OpenID flaking - good catch - merged the Sara-s - There Can Be Only One

ADD REPLYlink written 7.6 years ago by Istvan Albert ♦♦ 80k
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 2.3.0
Traffic: 1007 users visited in the last hour