SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq problem NicoBxl Bioinformatics 21 08-11-2012 10:57 AM
DESeq 1.5.30 - estimateDispersions horizon Bioinformatics 5 03-29-2012 12:08 PM
DESeq questions hanshart RNA Sequencing 10 03-07-2012 12:42 PM
DESeq: question about baseMean Azazel Bioinformatics 2 05-22-2011 06:40 PM
DESeq: question about baseMean. Also, replicates. Azazel Bioinformatics 5 05-18-2011 11:51 PM

Reply
 
Thread Tools
Old 09-09-2012, 05:41 PM   #1
Azazel
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 52
Question [DESeq] baseMean

I would like to ask how exacly baseMean is calculated in DESeq.

According to the manual from 2012-07-24 (http://tinyurl.com/8krvbny) the meaning of the three baseMean columns is:

Quote:
baseMean mean normalised counts, averaged over all samples from both conditions
baseMeanA mean normalised counts from condition A
baseMeanB mean normalised counts from condition B
I interpret this as:
Code:
baseMean = (baseMeanA + baseMeanB) / 2
However I have this result from a DE test which I did with DESeq (only first entry shown):
Code:
id	baseMean	baseMeanA	baseMeanB	foldChange	log2FoldChange	pval	padj
geneA	4.6870306180134	3.72258134619134	6.1337045257465	1.64770194532405	0.720455295377725	0.137300747770432	0.461541670327437
Code:
(3.72258134619134+6.1337045257465)/2 == 4.928143
4.928143 != 4.6870306180134
There is obviously a discrepancy. Maybe I am interpreting it wrong.

I'd be thankful if someone could explain how I can replicate the calculation of baseMean.

P.S.: I asked a similar question in May this year. Simon nicely answered me that:
Quote:
Yes, it is the normalized count, i.e., simply the number of reads divided by the size factor (normalization constant) of the sample
However, I do not know how I can get the "size factor" from the results table. If at all possible, I want to recompute baseMean given only baseMeanA and baseMeanB (is this even possible?). Because this way I wouldn't need to determine the sizeFactor.

Last edited by Azazel; 09-09-2012 at 05:45 PM. Reason: Added a PS for clarification.
Azazel is offline   Reply With Quote
Old 09-12-2012, 02:43 AM   #2
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

1. Maybe you have different numbers of samples in your two conditions. Then, the base mean will not be the average of the two per-condition means.

2. You real problem seems to be that you do not know how to get the size factors. If so, why do you not ask about this? Simply use the accessor function 'sizeFactors' on your CountDataSet to get a vector of size factors.
Simon Anders is offline   Reply With Quote
Reply

Tags
basemean, deseq, normalization

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:29 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO