runMEME,DNAStringSet fails when the seqs are un-named

[malcook]

Hi,

In my hands, calling runMEME fails unless sequences are named as it generates a temp fasta file lacking sequence identifiers, which breaks MEME.

Here's a test with an obvious "workaround"

library(assertthat)
seq<-readDNAStringSet(file.path(system.file("extdata", package="TFBSTools"), "crp0.s"))
names(seq)<-NULL
assert_that(is.error(try(motifSet <- runMEME(seq, 
                                             binary="meme", 
                                             arguments=list("-nmotifs"=3
                                                            )))))
names(seq)<-seq_along(seq)  # the workaround
assert_that(! is.error(try(motifSet <- runMEME(seq, 
                                             binary="meme", 
                                             arguments=list("-nmotifs"=3
                                                            )))))

[ge11232002]

Thanks. It's fixed in devel branch.