[Scipy-svn] r2937 - trunk/Lib/stats

scipy-svn@scip... scipy-svn@scip...
Tue Apr 24 17:02:10 CDT 2007


Author: oliphant
Date: 2007-04-24 17:02:08 -0500 (Tue, 24 Apr 2007)
New Revision: 2937

Modified:
   trunk/Lib/stats/stats.py
Log:
Use Pearson's definition of kurtosis in kurtosis-test.

Modified: trunk/Lib/stats/stats.py
===================================================================
--- trunk/Lib/stats/stats.py	2007-04-24 19:39:22 UTC (rev 2936)
+++ trunk/Lib/stats/stats.py	2007-04-24 22:02:08 UTC (rev 2937)
@@ -828,7 +828,7 @@
         warnings.warn(
             "kurtosistest only valid for n>=20 ... continuing anyway, n=%i" % 
             int(n))
-    b2 = kurtosis(a, axis)
+    b2 = kurtosis(a, axis, fisher=False)
     E = 3.0*(n-1) /(n+1)
     varb2 = 24.0*n*(n-2)*(n-3) / ((n+1)*(n+1)*(n+3)*(n+5))
     x = (b2-E)/np.sqrt(varb2)
@@ -848,8 +848,6 @@
     # fixme: find reference
     """Tests whether skew and/or kurtosis of dataset differs from normal curve.
 
-    This is the omnibus test of D'Agostino and Pearson, 1973
-
     Parameters
     ----------
     a : array



More information about the Scipy-svn mailing list