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
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
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
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
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
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
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)