Re: [Bioc-devel] very slow to use intronsByTranscript in GenomicFeatures

2013-12-20 Thread Ou, Jianhong
Thanks Herve, Robert, Now it works good. Yours sincerely, Jianhong Ou LRB 670A Program in Gene Function and Expression 364 Plantation Street Worcester, MA 01605 On 12/20/13 2:25 PM, "Hervé Pagès" wrote: >Hi Robert, Jianhong, > >This could be related to some changes to the relist() and sp

Re: [Bioc-devel] very slow to use intronsByTranscript in GenomicFeatures

2013-12-20 Thread Hervé Pagès
Hi Robert, Jianhong, This could be related to some changes to the relist() and split() code that I made a few days ago in IRanges. I didn't immediately make the corresponding changes to GenomicRanges and GenomicAlignments so relisting or splitting GRanges and GAlignments objects was broken for a

Re: [Bioc-devel] very slow to use intronsByTranscript in GenomicFeatures

2013-12-20 Thread Robert Castelo
hi, i can reproduce what Jianhong says, i noticed it earlier this week but didn't mention because we all know devel is a moving target and so on, but since this has been raised now i'll report what i'm getting. so, this is for Jianhong, if you downgrade the following packages to these partic

Re: [Bioc-devel] very slow to use intronsByTranscript in GenomicFeatures

2013-12-20 Thread Ou, Jianhong
In my case, looks like never end. I need to check my R first. Yours sincerely, Jianhong Ou LRB 670A Program in Gene Function and Expression 364 Plantation Street Worcester, MA 01605 On 12/20/13 12:05 PM, "Hervé Pagès" wrote: >Hi Jianhong, > >According to my timings, it's a little bit slow

Re: [Bioc-devel] very slow to use intronsByTranscript in GenomicFeatures

2013-12-20 Thread Hervé Pagès
Hi Jianhong, According to my timings, it's a little bit slower than exonsBy() but not that much. It has to do a little bit more work too as the introns are not explicitly stored in the SQLite db (the exons are) but are inferred from the exons and transcript boundaries. So intronsByTranscript() ha

Re: [Bioc-devel] very slow to use intronsByTranscript in GenomicFeatures

2013-12-20 Thread James W. MacDonald
Depends on what you mean by unacceptably slow. Took me almost 20 seconds, during which time I got pretty fidgety. > system.time(introns <- intronsByTranscript(TxDb.Hsapiens.UCSC.hg19.knownGene)) user system elapsed 18.741 0.204 19.340 How long did it take you, and how much RAM do you

[Bioc-devel] very slow to use intronsByTranscript in GenomicFeatures

2013-12-20 Thread Ou, Jianhong
Dear all, When I try to use intronsByTranscript to get introns for hg19 known genes, I found it is unacceptable slow. Does any body has the same problem? My code: library(GenomicFeatures) library(TxDb.Hsapiens.UCSC.hg19.knownGene) introns <- intronsByTranscript(TxDb.Hsapiens.UCSC.hg19.knownGene)