<p dir="ltr">Why don't you use CDF? </p>
<div class="gmail_quote">On Feb 4, 2016 3:22 PM, "Matthew Pitkin" <<a href="mailto:matthew.pitkin@glasgow.ac.uk">matthew.pitkin@glasgow.ac.uk</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear SciPy-dev list,<br>
<br>
I often require using trapezoidal integration of probability distribution functions (pdfs). For numerical reasons it is generally best to compute the logarithm of the pdfs, so integrating them whilst staying in log-space is generally required. Provided the steps between points are equal the scipy.misc.logsumexp function can be used for this, but I was wondering if a function specifically for this purpose (and following the usage of scipy.integrate.trapz function) might be more generally useful (called e.g. logtrapz)? I have my own code to do this, so if the list thinks this could be useful I can add it into scipy and create a pull request.<br>
<br>
Regards,<br>
<br>
Matt<br>
<br>
-- <br>
-----------------------------------------------------------<br>
Matthew Pitkin<br>
School of Physics and Astronomy<br>
Kelvin Building<br>
University of Glasgow<br>
Glasgow, G12 8QQ<br>
<br>
Tel. no.  <a href="tel:%2B44%20%280%29141%20330%207168" value="+441413307168" target="_blank">+44 (0)141 330 7168</a><br>
Mob. no.  <a href="tel:%2B44%20%280%29774%20050%207792" value="+447740507792" target="_blank">+44 (0)774 050 7792</a><br>
e-mail    <a href="mailto:matthew.pitkin@glasgow.ac.uk" target="_blank">matthew.pitkin@glasgow.ac.uk</a><br>
<a href="http://www.astro.gla.ac.uk/~matthew" rel="noreferrer" target="_blank">http://www.astro.gla.ac.uk/~matthew</a><br>
-----------------------------------------------------------<br>
<br>
_______________________________________________<br>
SciPy-Dev mailing list<br>
<a href="mailto:SciPy-Dev@scipy.org" target="_blank">SciPy-Dev@scipy.org</a><br>
<a href="https://mail.scipy.org/mailman/listinfo/scipy-dev" rel="noreferrer" target="_blank">https://mail.scipy.org/mailman/listinfo/scipy-dev</a><br>
</blockquote></div>