Question: EdgeR: different values in expression and differential expression tables
0
gravatar for Floydian_slip
5.7 years ago by
Floydian_slip130
United States
Floydian_slip130 wrote:

Hi, 
I am running EdgeR as a part of a bigger pipeline (CAP-miRSeq to be precise) where I quantify the expression and differential expression of miRNAs. However, the values in differential_expression tables does not match (though they are usually close) with values in expression tables (miRNA_expression_norm.xls, mature_miRNA_expression.xls and miRNA_expression_raw.xls). I am not sure which read count value were used in the differential expression tables. 

 

Does EdgeR also calculate the read count first before doing differential expression calculations? There is another tool in the big pipeline which has the capabilities to calculate read count. But I am not sure if the read counts come from that tool or from EdgeR. 

 

Can you please tell me where what could be the reason for different numbers in expression and differential expression tables? 

ADD COMMENTlink modified 5.7 years ago by Gordon Smyth2.0k • written 5.7 years ago by Floydian_slip130
0
gravatar for Sam
5.7 years ago by
Sam3.2k
New York
Sam3.2k wrote:

EdgeR take read count as an input instead of performing the counting itself. What EdgeR will do is that given the read count, it will calcuate the CPM for the statistical analysis.(Thanks for Devon's correction, EdgeR will not use CPM for statistic calculation) Maybe the value you see in the EdgeR output is the CPM instead of the raw input?

ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by Sam3.2k

Just for clarification, edgeR doesn't use CPM for statistics. You can filter genes by CPM metrics, which is likely what you're thinking of.

ADD REPLYlink written 5.7 years ago by Devon Ryan97k

opps, my bad. 

ADD REPLYlink written 5.7 years ago by Sam3.2k
0
gravatar for Gordon Smyth
5.7 years ago by
Gordon Smyth2.0k
Australia
Gordon Smyth2.0k wrote:

The differential expression (topTags) tables from edgeR don't contain expression values, so it would seem that you are comparing different things. The topTags tables do contain an Average Log CPM value, but this is an average for the miR across all samples, not an individual expression value.

edgeR analyses the counts but doesn't compute them. Many people use subread-featureCounts or htseq to compute read counts.

 

ADD COMMENTlink written 5.7 years ago by Gordon Smyth2.0k
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 2.3.0
Traffic: 2183 users visited in the last hour