Bioperl: relative-majority consensus, fast code sought

Gustavo Glusman bmgustav@bioinformatics.weizmann.ac.il
Tue, 6 Apr 1999 14:12:49 +0300


At 12:02 +0100 6/4/99, James Gilbert wrote:
>I agree, this is more maintainable.  To get it to
>work, however, you need to do:
>
>foreach $letter (qw/A T C G/) { $count{$letter} = eval
>"tr/$letter/$letter/"; }

Oops, my original error. Make it s/$letter/$letter/g, then...

     ____________________________________________________________
     Gustavo Glusman              Founder/administrator of BioMOO
     Gustavo@bioinfo.weizmann.ac.il    (public PGP key available)
     http://bioinfo.weizmann.ac.il/Gustavo
        Visit BioMOO, the biologists' virtual meeting place, at
     ___________ http://bioinfo.weizmann.ac.il/BioMOO ___________


=========== Bioperl Project Mailing List Message Footer =======
Project URL: http://bio.perl.org/
For info about how to (un)subscribe, where messages are archived, etc:
http://www.techfak.uni-bielefeld.de/bcd/Perl/Bio/vsns-bcd-perl.html
====================================================================