Genic annotations of whole hg38
1
0
Entering edit mode
bundo • 0
@6c4206de
Last seen 2.5 years ago
Japan

Hello,

I'm performing genic annotations with my data using annotatr, and I have a question. How can I know the proportion of genic annotations (exon, intron, intergenic and so on) of whole hg38 genome determined in annotatr?

Thank you.

hg38 genicannotations • 665 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 13 minutes ago
United States

What exactly do you mean by 'proportion of genic annotations'? One possibility would be to see what proportion of all the exons are covered by the regions you are annotating. If so, you end up with a GRanges object from annotatr, and you presumably have a GRanges object for the genomic sites you were annotating against. Let's say your object from annotatr is called 'annot' and the regions GRanges object is called 'exons'.

We don't want to count the width of exons or annotations that are duplicated, so we want to use reduce first. Then we get the widths, of those regions, and then sum it up.

## total width of annotated regions
totannot <- sum(width(reduce(annot)))
## total width of exons
totexon <- sum(width(reduce(exons)))
totannot/totexon
ADD COMMENT

Login before adding your answer.

Traffic: 978 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6