Re: [Bioc-devel] export() to 2bit file

2014-05-12 Thread Hervé Pagès
On 05/12/2014 12:23 PM, Michael Lawrence wrote: On Mon, May 12, 2014 at 11:41 AM, Hervé Pagès mailto:hpa...@fhcrc.org>> wrote: Hi Michael, On 05/09/2014 04:39 PM, Michael Lawrence wrote: What would be the fastest way to do this with a DNAString? Just an alphabetFre

Re: [Bioc-devel] export() to 2bit file

2014-05-12 Thread Michael Lawrence
On Mon, May 12, 2014 at 11:41 AM, Hervé Pagès wrote: > Hi Michael, > > > On 05/09/2014 04:39 PM, Michael Lawrence wrote: > >> What would be the fastest way to do this with a DNAString? Just an >> alphabetFrequency? >> > > That would do it. > > A couple of other issues I ran into with the 2bit

Re: [Bioc-devel] export() to 2bit file

2014-05-12 Thread Hervé Pagès
Hi Michael, On 05/09/2014 04:39 PM, Michael Lawrence wrote: What would be the fastest way to do this with a DNAString? Just an alphabetFrequency? That would do it. A couple of other issues I ran into with the 2bit code: (1) It fails on empty sequences: > export(DNAStringSet(c("AA", "",

Re: [Bioc-devel] export() to 2bit file

2014-05-09 Thread Michael Lawrence
What would be the fastest way to do this with a DNAString? Just an alphabetFrequency? On Fri, May 9, 2014 at 4:07 PM, Hervé Pagès wrote: > Hi Michael, > > library(rtracklayer) > library(Biostrings) > x <- DNAStringSet("AAA-CCC-GGG-TTT-NNN-KKK") > > Then: > > > x > A DNAStringSet

[Bioc-devel] export() to 2bit file

2014-05-09 Thread Hervé Pagès
Hi Michael, library(rtracklayer) library(Biostrings) x <- DNAStringSet("AAA-CCC-GGG-TTT-NNN-KKK") Then: > x A DNAStringSet instance of length 1 width seq [1]23 AAA-CCC-GGG-TTT-NNN-KKK > export(x, "x.2bit") > import("x.2bit") A DNAStringSet instance of length 1