[Shootout-list] Re: OCaml k-nucleotide

[Joel Hoffman]

Daniel South wrote:

>Joel, are you using a correctly generated fasta file?
>Change the output line in the Tcl program from
>
> puts [format "%s %0.2f" [lindex $item 0] $percent]
>
>to
>
> puts "[format "%s %0.2f %0.12f" [lindex $item 0] $percent 
>$percent]\t$sum\t[lindex $item 1]"
>
>this will give you, after the usual output, percent out to 12 decimal 
>places, then the sum, then the count.
>
>When the modified Tcl program is run against your input file, do the 
>sums and counts match?
>
>Daniel South 
>  
>

Yes, here's the output:

T 30.41 30.408000000000 100000  30408
A 30.30 30.305000000000 100000  30305
C 19.65 19.652000000000 100000  19652
G 19.63 19.635000000000 100000  19635


I'm using the file knucleotide-input20000.txt, from the shootout CVS. 
(in /shootout/website/desc/).

joel