id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
5095,incorrect results from hypergeometric pdf,David Koes <dkoes@…>,John Maddock,"#include <boost/math/distributions/hypergeometric.hpp>
#include <boost/math/policies/policy.hpp>
#include <iostream>

using namespace std;
using namespace boost;

int main()
{
  unsigned N = 16086184;
  unsigned n = 256004;
  unsigned Q = 251138;
  math::hypergeometric_distribution<double> hyper(n, Q, N);
  cout <<  math::pdf<double>(hyper, 4000) << "" "" <<  math::pdf<double>(hyper, 4001) << "" "" <<  math::pdf<double>(hyper, 4002) << ""\n"";

  return 0;

}

Output:
0.00640003 1.11519e-09 0.00638443
The value for 4001 is incorrect (according to http://stattrek.com/Tables/Hypergeometric.aspx).  In fact, every value where k is odd appears to be incorrect.  ",Bugs,closed,To Be Determined,math,Boost 1.45.0,Problem,fixed,,