[Scipy-tickets] [SciPy] #1262: scipy.stats.maxwell.ppf(0.25)=0?

SciPy Trac scipy-tickets@scipy....
Fri Aug 13 12:37:31 CDT 2010


#1262: scipy.stats.maxwell.ppf(0.25)=0?
-----------------------+----------------------------------------------------
 Reporter:  pozharski  |       Owner:  somebody
     Type:  defect     |      Status:  new     
 Priority:  normal     |   Milestone:  0.9.0   
Component:  Other      |     Version:  0.7.0   
 Keywords:             |  
-----------------------+----------------------------------------------------

Comment(by josefpktd):

 here are the calls to gammaincinv for the example and I don't see any
 problems:

 {{{
 >>> from scipy import special
 >>> q=0.25; np.sqrt(2*special.gammaincinv(1.5,q))
 1.1011507176793143

 >>> q=np.linspace(0,1,41)
 >>> q
 array([ 0.   ,  0.025,  0.05 ,  0.075,  0.1  ,  0.125,  0.15 ,  0.175,
         0.2  ,  0.225,  0.25 ,  0.275,  0.3  ,  0.325,  0.35 ,  0.375,
         0.4  ,  0.425,  0.45 ,  0.475,  0.5  ,  0.525,  0.55 ,  0.575,
         0.6  ,  0.625,  0.65 ,  0.675,  0.7  ,  0.725,  0.75 ,  0.775,
         0.8  ,  0.825,  0.85 ,  0.875,  0.9  ,  0.925,  0.95 ,  0.975,  1.
 ])
 >>> np.sqrt(2*special.gammaincinv(1.5,q))
 array([ 0.        ,  0.46453771,  0.59316635,  0.68701716,  0.76444383,
         0.83208036,  0.89318052,  0.94961712,  1.00258367,  1.05289766,
         1.10115072,  1.14779054,  1.19316899,  1.23757213,  1.28123987,
         1.32437949,  1.36717534,  1.40979603,  1.45240026,  1.49514148,
         1.53817225,  1.58164821,  1.62573222,  1.67059893,  1.71643994,
         1.76347021,  1.81193611,  1.86212611,  1.91438522,  1.96913543,
         2.02690526,  2.08837439,  2.1544437 ,  2.22635082,  2.30587247,
         2.39570729,  2.50027771,  2.62766894,  2.79548348,  3.05751592,
                Inf])
 }}}

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1262#comment:2>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list