[Shootout-list] fasta data error
Simon Geard
simon@whiteowl.co.uk
Wed, 19 Jan 2005 21:25:14 +0000
The second set of fasta data, iub, seems to be wrong. The cumulative
probability is > 1
iub = (/ pair('a', 0.270d0), &
pair('c', 0.125d0), &
pair('g', 0.125d0), &
pair('t', 0.270d0), &
pair('B', 0.02d0), &
pair('D', 0.02d0), &
pair('H', 0.02d0), &
pair('K', 0.02d0), &
pair('M', 0.02d0), &
pair('N', 0.02d0), &
pair('R', 0.02d0), &
pair('S', 0.02d0), &
pair('V', 0.02d0), &
pair('W', 0.02d0), &
pair('Y', 0.02d0) /)
call makeCumulative(iub)
print *,iub(size(iub))%p
outputs 1.01. What are the correct values?
Simon