id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
8308,Poisson Quantile function returns wrong results in certain cases,Caleb.Welton@…,John Maddock,"Sample use case:

{{{
main() {
  double pi = 3.1415927;
  boost::math::poisson poisson(pi);

  //  std::cout << quantile(poisson, 0.00) << "" "";  
  std::cout << quantile(poisson, 0.05) << "" "";
  std::cout << quantile(poisson, 0.30) << "" "";
  std::cout << quantile(poisson, pi/6) << "" "";
  std::cout << quantile(poisson, 0.7) << "" "";
  std::cout << quantile(poisson, 1-1e-10) << "" "";
  // std:: cout << quantile(poisson, 1) << std::endl;  
}
}}}
Produces the results
(error), 0, 1, 3, 4, 20, (error)

Similar code in R or MATLAB results in: (0),1,2,3,4,20,(Inf)

{{{
R> qpois(c(0,0.05,0.3,pi/6,0.7,1-1e-10,1),pi)
[1]   0   1   2   3   4  20 Inf

MATLAB> f := stats::poissonQuantile(PI):
MATLAB> f(0), f(1/20), f(0.3), f(PI/6), f(0.7), f(1-1/10^10), f(1)
0,1,2,3,4,20,Inf
}}}

The errors on the edge cases are less troubling than the seemingly incorrect results at low probabilities.

",Bugs,closed,To Be Determined,math,Boost 1.52.0,Problem,fixed,,