SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA seq data normalization question slny Bioinformatics 35 10-19-2016 05:32 AM
edgeR ECHo Bioinformatics 13 04-18-2013 01:01 PM
edgeR Puva Bioinformatics 2 05-19-2011 09:04 AM
Normalization: DESeq vs. EdgeR with method="RLE" aquila RNA Sequencing 3 04-26-2011 09:00 AM
DEGseq or edgeR mmanrique Bioinformatics 10 02-12-2010 02:13 PM

Reply
 
Thread Tools
Old 03-18-2012, 11:05 AM   #1
polsum
Member
 
Location: Texas

Join Date: May 2009
Posts: 32
Default edgeR Normalization question

The edgeR manual gives us three examples for DE analysis with each example normalizing the data in different manner.

In the case study 1, the Normalization was done by calcNormFactors. In the case study 2 however, only counts per million or cpm was calculated and no calcNormFactors was used. In the case study 3, both cpm and calcNormFactors was used. In this 3rd example, we are essentially modifying the data two times before performing DE analysis which is strange to me.

I basically want to perform miRNA DE analysis between two samples (each with 3 replicates). Which of the above methods is applicable for me?

Can any one please educate me? thanks in advance.
polsum is offline   Reply With Quote
Old 05-07-2012, 12:21 AM   #2
Gordon Smyth
Member
 
Location: Melbourne, Australia

Join Date: Apr 2011
Posts: 91
Default

Dear polsum,

The edgeR User's Guide contains five case studies, and all are normalized exactly the same way, except for the second which deals with deepSAGE data for which simple library size normalization was considered sufficient.

You are mis-understanding the role of cpm, which are used only as descriptive quantities for data exploration in edgeR, not as part of the formal normalization or differential expression process.

We obviously recommend, unless there are good reasons otherwise, that you use calcNormFactors() for an RNA-Seq analysis.

Note that edgeR is a Bioconductor package, so the way to get a prompt answer is to post your question to the Bioconductor mailing list. See the section "how to get help" in the edgeR User's Guide.

It also seems that your version of edgeR must be very old. Could I suggest that you update?

Gordon
Gordon Smyth is offline   Reply With Quote
Reply

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 11:20 PM.


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