scipy-svn@scip... scipy-svn@scip...
Wed Jul 18 22:52:16 CDT 2007

```Author: pierregm
Date: 2007-07-18 22:52:07 -0500 (Wed, 18 Jul 2007)
New Revision: 3174

Modified:
Log:
MA:
morestats      : fixed hdquantiles for empty dataset
TS:
tseries        : fixed compressed w/ 2D series with one date
added support for optional arguments/paremeters in convert

===================================================================
--- trunk/Lib/sandbox/maskedarray/morestats.py	2007-07-18 13:12:12 UTC (rev 3173)
+++ trunk/Lib/sandbox/maskedarray/morestats.py	2007-07-19 03:52:07 UTC (rev 3174)
@@ -58,13 +58,15 @@
"""
def _hd_1D(data,prob,var):
"Computes the HD quantiles for a 1D array."
-        xsorted = numpy.sort(data.compressed().view(ndarray))
+        xsorted = numpy.squeeze(numpy.sort(data.compressed().view(ndarray)))
n = len(xsorted)
#.........
hd = empty((2,len(prob)), float_)
if n < 2:
hd.flat = numpy.nan
-            return hd
+            if var:
+                return hd
+            return hd[0]
#.........
v = arange(n+1) / float(n)
betacdf = beta.cdf
@@ -87,7 +89,7 @@
p = numpy.array(prob, copy=False, ndmin=1)
# Computes quantiles along axis (or globally)
if (axis is None):
-        result = _hd_1D(data.compressed(), p, var)
+        result = _hd_1D(data, p, var)
else:
assert data.ndim <= 2, "Array should be 2D at most !"
result = apply_along_axis(_hd_1D, axis, data, p, var)

```