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