Hi all,
I am now using raw reads to calculate the differential gene expression. I see DESeq gives SizeFactors for each sample and can use that for normalization. But it seems DESeq only use that normalization for visualization, not for calculate the differential pvalue. Is that correct?
I've read DESeq paper and still not very clear about that. I know DESeq use NB distribution to calculate the pvalue. But for raw reads, it must be normalized first, right?? Thanks,
Che
2 answers
I'm curious what makes you think that the sizeFactors are only used for visualization?
If you download the source code to for the package and poke around a bit, you'll see that the sizeFactors are used everywhere from dispersion estimation, nbinomTest-ing, fitting GLMs, getting variance stabilized data, etc.
Have you read the DESeq vignette? Fold changes are calculated for you after a call to nbinomTest, for example. To just get something like a matrix of expression that you can calculate fold changes for data exploration, look at the example in section 6, where they shoot their data through a variance stabilizing transformation.
Log in to answer this question.
-
Exporting normalized counts data from DeSeq2
written by mohammedtoufiq91
-
Best normalization method to use with ChIP-seq counts
written by science_lizard •
-
p-value in Limma vs Graphpad
written by LuisNagano •
-
Differential Gene Expression Analysis From TMM normalized counts
written by boyinepallyki •
-
GFOLD diff input
written by tonja.r
-
RNA-seq counts analysis with DESeq
written by iside •
-
Z-Score On Raw Counts
written by camelbbs
-
Get Conservation Score Of Chromosome Region
written by camelbbs
-
Draw Heatmap Or Do Pca Analysis With Raw Read Counts?
written by camelbbs
-
Normalization Method In Comparing Differential Gene Expression
written by camelbbs
