[BioC] Faster way to apply a function than endoapply on 40K GrangesList

Lakshmanan Iyer laxvid at gmail.com
Wed Aug 10 21:36:54 CEST 2011


Hi,
I am trying to add 10K to the end of the intervals in a GRangesList
object using endoapply and flank using the comand:
p3UTRs10K <- endoapply(p3UTRs,function (z) flank (z, 10000,start=FALSE));

I am writing to find out if someone has a suggestion for a faster way
of doing this?  With 40000 elements it is taking a long time..

Here is what I am doing:

class (p3UTRs)
[1] "GRangesList"
attr(,"package")
[1] "GenomicRanges"

length (p3UTRs)
[1] 41053

p3UTRs10K <- endoapply(p3UTRs,function (z) flank (z, 10000,start=FALSE));

> sessionInfo()
R version 2.13.0 (2011-04-13)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
 [1] BSgenome.Mmusculus.UCSC.mm9_1.3.17 GenomicFeatures_1.4.0
 [3] leeBamViews_0.99.11                BSgenome_1.20.0
 [5] Biobase_2.12.1                     Rsamtools_1.4.1
 [7] Biostrings_2.20.0                  GenomicRanges_1.4.3
 [9] IRanges_1.10.0                     rtracklayer_1.12.0
[11] RCurl_1.5-0                        bitops_1.0-4.1

loaded via a namespace (and not attached):
[1] biomaRt_2.8.0 DBI_0.2-5     RSQLite_0.9-4 tools_2.13.0  XML_3.2-0
>



More information about the Bioconductor mailing list