Entering edit mode
AntonS • 0@antons-13533
Last seen 5.7 years ago
I tried to find the solution on my own, but I think the whole tdx, annotationdb, ... are circuitous.
I have a chromosomal location (chr1:48902063-48902085, chr2:202340342-202340364, ....) and I just want to know, whether this region is an exon. Is there a simple command to do it?
Best regards and thank you in advanced
If you want to add a comment, please use the ADD COMMENT link. The Add your answer box below is for people to add answers, not more questions or comments.
I have no idea what you mean by 'nested in at least 3 levels'. Are you complaining that it took three lines of code to generate an answer?
If so, there are always tradeoffs to be made - you can make things really simple and straightforward, but the cost to that is you force people to do what you think they should do, and make it difficult to do other things that they might want to do. The alternative is to make things very powerful, but the cost to that is complexity. All of the Bioconductor infrastructure for dealing with genomic position data is very powerful, but at the same time very complex. What people gain from the complexity is the ability to do lots of things that a simpler API would likely prevent.
There is extensive documentation for all of the objects that I have generated, so I would point you to the help pages, as well as the vignettes for the GenomicRanges package. As a hint towards what you want, do note that
Gives you a GRanges object, where the names of the GRanges object are the Entrez Gene ID, which is, I presume, what you wanted for the Gene column.
as.data.frame()to get the simple data.frame that the user desires.