[Scipy-tickets] [SciPy] #1802: Infinite time taken by scipy hypergeom function for some particular values
SciPy Trac
scipy-tickets@scipy....
Tue Jan 1 08:17:11 CST 2013
#1802: Infinite time taken by scipy hypergeom function for some particular values
--------------------+-------------------------------------------------------
Reporter: imsc | Owner: somebody
Type: defect | Status: new
Priority: normal | Milestone: Unscheduled
Component: Other | Version: 0.11.0
Keywords: |
--------------------+-------------------------------------------------------
Comment(by josefpktd):
Looks like a bug to me. It's strange that it works (althoug slowly) for
0.94 and 0.955, but not for 0.945 and 0.95.
The problem is that the calculations are done by the generic methods, and
they are difficult to get right for most cases.
_sf got a special implementation in recent scipy, but it is not used for
the ``interval`` and ``ppf``, and also not for the isf, calculations.
one other strange part: .a is negative (which might throw of the root
finding)
{{{
>>> stats.hypergeom.b
array(757)
>>> stats.hypergeom.a
-143652
}}}
--
Ticket URL: <http://projects.scipy.org/scipy/ticket/1802#comment:1>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.
More information about the Scipy-tickets
mailing list