[SciPy-User] log pdf, cdf, etc

[Chris Strickland]

Hi,

When using any of the distributions of scipy.stats there does not seem to be

the ability (or at least I cannot figure out how) to have the function
return
the log of the pdf, cdf, sf, etc. For statistical analysis this is
essential.
For instance suppose we are interested in an exponential distribution for a
random variable x with a hyperparameter lambda there needs to be an option
that returns -log(lambda)-x/lambda. It is not sufficient (numerically) to
calculate log(scipy.stats.expon.pdf(x,lambda)).

Is there a way to do this using the distributions in scipy.stats?

If there is not is it possible for me to suggest that this feature is added.

There is such an excellent range of distributions, each with such an
impressive range of options, it seems ashame to have to mostly manually code

up the log of pdfs and often call the log of CDFs from R.

Thanks,
Chris.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20100528/1df2721c/attachment.html>