[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