[BioC] Computing normalized counts for exons in edgeR
Ryan
rct at thompsonclan.org
Sun Mar 30 08:21:41 CEST 2014
There is no way to derive exon-level counts from gene-level counts. You
will have to go back to the bam files and perform a separate counting
step for exons.
On Sat Mar 29 17:48:02 2014, Nick N 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?
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list