[Scipy-tickets] [SciPy] #488: scipy.stats.expon is broken for x=0

SciPy scipy-tickets@scipy....
Sun Aug 26 08:53:14 CDT 2007


#488: scipy.stats.expon is broken for x=0
-------------------------+--------------------------------------------------
 Reporter:  vomjom       |       Owner:  somebody
     Type:  defect       |      Status:  new     
 Priority:  normal       |   Milestone:  0.6     
Component:  scipy.stats  |     Version:  0.5.2   
 Severity:  normal       |    Keywords:          
-------------------------+--------------------------------------------------
 >>> scipy.stats.expon.pdf(0)
 array(0.0)

 The pdf for the exponential distribution is:
 (1/scale)*(e**(-(1/scale)*x)

 For x=0 this should return 1/scale, not 0.0.

 Strange enough, expon._pdf works, so it's probably a bug in rv_continuous.

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


More information about the Scipy-tickets mailing list