Question: Computing Normalized Counts For Exons In Edger
0
gravatar for Nick
6.3 years ago by
Nick270
Spain
Nick270 wrote:

I've done standard rna-seq differential expression analysis on a gene level using edgeR. For downstream analysis I would like to extract the normalized counts for a particular exon. What is the best way to go about it? Can I derive the normalized counts for that exon straight from my analysis (e.g. by using, say, cpm())?

The first workaround that pops to my mind is that I could simply edit either the GTF or the raw count data to discard any counts for other exons of that gene, than simply run the standard count normalisation and I would get the normalized counts for that exon masquerading as the normalized count for the corresponding gene. But this seems too hacky and dirty solution. Can anyone give advice on how to do this properly?

edger rna-seq • 2.2k views
ADD COMMENTlink modified 6.3 years ago by Devon Ryan96k • written 6.3 years ago by Nick270
1
gravatar for Devon Ryan
6.3 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

You can't derive exon-level metrics from gene-level metrics. DEXSeq comes with some scripts meant to count reads mapping to exonic bins, so just use that. I suspect that featureCounts can do this as well, but have never tried.

ADD COMMENTlink written 6.3 years ago by Devon Ryan96k
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: 1082 users visited in the last hour