[Seqinr-forum] Bug in comp()

[Darren Obbard]

Hi,

There is a bug in comp

# y and r are missing from the list of all bases, and return NA
#---------------------------
require("seqinr")
example(comp)
allbases
comp(allbases)
comp(allbases,ambiguous=TRUE)
allbases<-c(allbases,"y","r")
comp(allbases,ambiguous=TRUE)
#---------------------------
#

Darren
-- 

Darren Obbard
Institute of Evolutionary Biology
Room 123, Ash 2, Ashworth Labs
Kings Buildings
University of Edinburgh, UK
darren.obbard at ed.ac.uk


The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.