[Scipy-svn] r6430 - trunk/scipy/stats

[scipy-svn at scipy.org]

Author: oliphant
Date: 2010-05-29 18:29:26 -0500 (Sat, 29 May 2010)
New Revision: 6430

Modified:
   trunk/scipy/stats/distributions.py
Log:
Fix typo in F distribution.

Modified: trunk/scipy/stats/distributions.py
===================================================================
--- trunk/scipy/stats/distributions.py	2010-05-29 23:22:56 UTC (rev 6429)
+++ trunk/scipy/stats/distributions.py	2010-05-29 23:29:26 UTC (rev 6430)
@@ -2206,8 +2206,8 @@
     def _logpdf(self, x, dfn, dfd):
         n = 1.0*dfn
         m = 1.0*dfd
-        lPx = m/2 * log(m) + n/2 * log(n) + (n/2-1)*log(x)
-        lPx -= (n+m)/2*log(m+n*x) - special.betaln(n/2,m/2)
+        lPx = m/2*log(m) + n/2*log(n) + (n/2-1)*log(x)
+        lPx -= ((n+m)/2)*log(m+n*x) + special.betaln(n/2,m/2)
         return lPx
     def _cdf(self, x, dfn, dfd):
         return special.fdtr(dfn, dfd, x)