Question: What Is The Effect Of Using Ordering In The Exonby Function In R (Bioconductor)
gravatar for Jetse
6.6 years ago by
Jetse0 wrote:

Hello everyone,

I am using the programming language R, with this language I first want to download exons from the UCSC website. Now I'm using this R code:

txdb <- makeTranscriptDbFromUCSC( genome='hg19', tablename='refGene' )
tx_by_exon <- exonsBy(txdb, 'tx')

When I'm changing the 'tx' to 'gene', the data is totally differend. When converting to a data.frame (with, the by gene ordered list is only half of the unordered list... In the manual this is described as: One of ‘"gene"’, ‘"exon"’, ‘"cds"’ or ‘"tx"’. Determines the grouping. And I thought grouping doesn't change the data...

I use this data to check which mapping tool to use, tophat or bwa. The difference between the ordering, makes or tophat much better (when not ordering), or bwa just a little bit better (when ordering by gene)... The largest difference between those mapping tools is the RNA splicing. So this change in ordering has something to do with the RNA splicing...

Anyone knows what this difference is?

R exon • 1.7k views
ADD COMMENTlink modified 6.6 years ago by Jeremy Leipzig18k • written 6.6 years ago by Jetse0
gravatar for Jeremy Leipzig
6.6 years ago by
Philadelphia, PA
Jeremy Leipzig18k wrote:

grouping changes the data

you can expect a lot more exons by transcript because multiple transcripts exist for the same gene, so you will see that same exon appear multiple times under different transcript names

ADD COMMENTlink written 6.6 years ago by Jeremy Leipzig18k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1070 users visited in the last hour