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
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
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", "",
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
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