Re: [Bioc-devel] equivalent of seqselect,vector,Rle

2013-12-13 Thread Hervé Pagès
Hi Michael, On 12/13/2013 06:39 PM, Michael Lawrence wrote: Coercion might suffice. I do remember Patrick optimizing these selections with e.g. memcpy(), so they are pretty fast. The memcpy() trick was used (and is still used in extractROWS) when seqselect'ing by a Ranges object. For subsettin

Re: [Bioc-devel] equivalent of seqselect,vector,Rle

2013-12-13 Thread Michael Lawrence
Coercion might suffice. I do remember Patrick optimizing these selections with e.g. memcpy(), so they are pretty fast. No profiling data though. I do have some performance critical code that has relied on the Rle-based extraction. Would be nice to avoid re-evaluating the performance. On Fri, Dec

Re: [Bioc-devel] equivalent of seqselect,vector,Rle

2013-12-13 Thread Hervé Pagès
On 12/13/2013 01:49 PM, Michael Lawrence wrote: Thanks, makes sense. Didn't realize we could dispatch on the 'i' parameter. I sort of recall the perception that we couldn't, and that was one of the main motivations behind seqselect. But it does appear possible. Well I was hoping I could do this

Re: [Bioc-devel] equivalent of seqselect,vector,Rle

2013-12-13 Thread Michael Lawrence
Thanks, makes sense. Didn't realize we could dispatch on the 'i' parameter. I sort of recall the perception that we couldn't, and that was one of the main motivations behind seqselect. But it does appear possible. Michael On Fri, Dec 13, 2013 at 1:10 PM, Hervé Pagès wrote: > Hi Michael, > > >

Re: [Bioc-devel] equivalent of seqselect,vector,Rle

2013-12-13 Thread Hervé Pagès
Hi Michael, On 12/13/2013 01:03 PM, Michael Lawrence wrote: I used to use seqselect for subsetting ordinary R vectors by Ranges and Rle. IRanges:::extractROWS does this, but it's hidden behind the namespace. What is the public way of doing this? Maybe we just need to export extractROWS()? Or so

[Bioc-devel] equivalent of seqselect,vector,Rle

2013-12-13 Thread Michael Lawrence
I used to use seqselect for subsetting ordinary R vectors by Ranges and Rle. IRanges:::extractROWS does this, but it's hidden behind the namespace. What is the public way of doing this? Maybe we just need to export extractROWS()? Or something with a better name? Michael [[alternative HTM

Re: [Bioc-devel] error when install ŒRsamtools¹, ŒGenomicAlignments¹ and Œrtracklayer¹

2013-12-13 Thread Ou, Jianhong
Hi Herve, Thanks for your reply. I updated my R to yesterday's building and it works great. Looks like it is the problem of my old version of R 3.10. Yours sincerely, Jianhong Ou LRB 670A Program in Gene Function and Expression 364 Plantation Street Worcester, MA 01605 On 12/13/13 2:50 PM

Re: [Bioc-devel] error when install ŒRsamtools¹, ŒGenomicAlignments¹ and Œrtracklayer¹

2013-12-13 Thread Hervé Pagès
Hi Jianhong, On 12/13/2013 11:45 AM, Ou, Jianhong wrote: Dear all, When I try to renew Rsamtools, GenomicAlignments and rtracklayer, I got error Error : objects ‘GAlignments’, ‘GAlignmentPairs’ are not exported by 'namespace:GenomicRanges' You didn't show us what you did exactly. Try to do

[Bioc-devel] error when install ŒRsamtools¹, ŒGenomicAlignments¹ and Œrtracklayer¹

2013-12-13 Thread Ou, Jianhong
Dear all, When I try to renew Rsamtools, GenomicAlignments and rtracklayer, I got error Error : objects ‘GAlignments’, ‘GAlignmentPairs’ are not exported by 'namespace:GenomicRanges' Does anybody know how to figure out this problem? > biocValid() * sessionInfo() R Under development (unstable

Re: [Bioc-devel] Consistency between miRNATarget and MiRaGE

2013-12-13 Thread Hervé Pagès
Hi tag, The data-experiment repository doesn't use a trunk/branch organization like the software repo, so the devel and release builds both point to trunk. You can see this in the "svn info" box ("URL" field) at the top of the HTML reports: http://bioconductor.org/checkResults/2.14/data-experi