[BioC] consensusString Function
Erik Wright
eswright at wisc.edu
Tue Apr 6 22:31:10 CEST 2010
Hello,
I am trying to get a consensus string for a DNAStringSet, but I am getting an error. The documentation for consensusString says the argument "x" is either a consensus matrix or an XStringSet. So this should work, right?:
> myDNAStringSet <- DNAStringSet(c("NNNN","ACTG"))
> consensusString(myDNAStringSet)
Error in .local(x, ...) :
'threshold' must be a numeric in (0, 1/sum(rowSums(x) > 0)]
Specifying a threshold in the arguments doesn't seem to make a difference.
Thanks!,
Erik
More information about the Bioconductor
mailing list