[BioRuby] Consensus sequence

Tomoaki NISHIYAMA tomoakin at kenroku.kanazawa-u.ac.jp
Thu Aug 5 00:28:04 UTC 2010


Hi,

> 2) How can the count of the greatest residue (b[0][1]) be larger  
> than or equal to the total number of residues?


It is obvious that the count of the greatest residue is equal to the  
total number of residues
if all the residues are identical.

Presumably, the parameter threshold should be
0 <= threshold <= 1.0

> 1) Why is it sorting? Shouldn't it use a linear search?


I really don't know.  So this is just my feeling,
but this could be simplicity, readability, and extensibility...
Did you compare the performance with linear search?
-- 
Tomoaki NISHIYAMA

Advanced Science Research Center,
Kanazawa University,
13-1 Takara-machi,
Kanazawa, 920-0934, Japan




More information about the BioRuby mailing list