[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