Question: use only lib size as normalizing factor in edgeR
0
gravatar for Udi Landau
10 months ago by
Udi Landau30
Udi Landau30 wrote:

Hello, 

I would like to use only the lib size as normalizing factor in edgeR, so I am doing this:

cds <- DGEList(counts = counts, group = group, norm.factors = MyNormFactors)

I skip the 

cds <- calcNormFactors(cds, method = "none")

Because It will bring the norm factor back to 1..

and I continue with calculating dispersion and doing the exact test.

What is seems to me strange is that the pseudo count look the same as the original count.

Isn't it supposed to change?

When I use no normalization (factor = 1)  the pseudo count do change a little..

Just want to know if this is the right way so I could trust my results.. 

edger • 233 views
ADD COMMENTlink modified 10 months ago by Aaron Lun24k • written 10 months ago by Udi Landau30
Answer: use only lib size as normalizing factor in edgeR
2
gravatar for Aaron Lun
10 months ago by
Aaron Lun24k
Cambridge, United Kingdom
Aaron Lun24k wrote:

There's several issues in your post.

  • If you want to do library size normalization, then don't specify norm.factors in the DGEList call. Having all normalization factors equal to 1 means that edgeR will use the log-library sizes directly as offsets, without calculating effective library sizes.
  • It doesn't matter if you skip calcNormFactors or not, your normalization factors should already be 1.
  • All currently recommended edgeR pipelines use glmQLFit and glmQLFTest instead of exactTest. The quasi-likelihood GLM framework is more accurate, more flexible and faster than "classic" edgeR - better in every way.
ADD COMMENTlink modified 10 months ago • written 10 months ago by Aaron Lun24k

Thank you!

Can you tell me what would you recommend to use for plotting the data?

I understand that the pseudo count is not suitable? 

ADD REPLYlink written 10 months ago by Udi Landau30
1

See Section 2.15 of the edgeR User's Guide ("Clustering, heatmaps etc").

ADD REPLYlink modified 10 months ago • written 10 months ago by Gordon Smyth37k
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 16.09
Traffic: 113 users visited in the last hour