[BioC] DESeq/DESeq normalization on different experiments

Gilgi Friedlander gilgi.friedlander at weizmann.ac.il
Thu May 22 20:13:37 CEST 2014


Thanks a lot!

-----Original Message-----
From: Michael Love [mailto:michaelisaiahlove at gmail.com] 
Sent: Thursday, May 22, 2014 6:10 PM
To: Ryan
Cc: Gilgi [guest]; Gilgi Friedlander; bioconductor at r-project.org
Subject: Re: [BioC] DESeq/DESeq normalization on different experiments

hi Gilgi,

Yes I would second Ryan's recommendation to not normalize all the experiments together.

Mike

On Thu, May 22, 2014 at 11:06 AM, Ryan <rct at thompsonclan.org> wrote:
> Hi Gilgi,
>
> If you are not going to test for differential expression between 
> experiments, then there is no purpose in normalizing them together. 
> The more worrying problem with analyzing all your experiments as a 
> single data set is that a single dispersion value will be estimated 
> for each gene across all experiments. This is only ok if you believe 
> that every gene has equal biological variability in all your 
> experiments, which is unlikely to be the case.
>
> -Ryan
>
>
> On Thu May 22 01:26:55 2014, Gilgi [guest] wrote:
>>
>> Hello,
>>
>> I have ~100 RNA seq-samples from the same organism, but they include 
>> different experiments (each experiment of 6-16 samples). I would like 
>> to put all together into our RNA seq pipeline, that includes  
>> mapping, htseq count and DESeq. Differential expression comparisons 
>> will be done of course only between samples of a single experiment (due to possible batch effects).
>> However, will it be a good idea to normalize all samples together?
>> The assumption behind the normalization is that most genes are not 
>> differentially expressed, but between different experiments the 
>> variability might be higher. Therefore I wanted to ask your opinion 
>> on doing the normalization on such a large combined data-set.
>>
>> Thank you,
>> Gilgi
>>
>>   -- output of sessionInfo():
>>
>> none
>>
>> --
>> Sent via the guest posting facility at bioconductor.org.
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor


More information about the Bioconductor mailing list