[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