[Numpy-discussion] Medians that ignore values
Thu Sep 18 13:23:25 CDT 2008
On Thursday 18 September 2008 13:31:18 Peter Saffrey wrote:
> The version in the Ubuntu package repository. It says 1:1.0.4-6ubuntu3.
So it's 1.0 ? It's fairly old, that would explain.
> > if you don't give an axis
> > parameter, you should get the median of the flattened array, therefore a
> > scalar, not an array.
> Not for my version.
Indeed. Looks like the default axis changed from 0 in 1.0 to None in the
incoming 1.2. But that's a detail at this point.
> > Anyway: you should use ma.median for masked arrays. Else, you're just
> > keeping the NaNs where they were.
> That will be the problem. My version does not have median or mean methods
> for masked arrays, only the average() method.
The method mean has always been around for masked arrays, so has the
corresponding function. But I'm surprised, median has been in numpy.ma.extras
for a while. Maybe not 1.0...
> According to this page:
> 1.1.0 is the latest release.
You need to update your internet ;) 1.1.1 was released 6 weeks ago.
> Do I need to use an SVN build to get the
> ma.median functionality?
No, you can install 1.1.1, that should work.
Note that I just fixed a bug in median in SVN (it would fail when trying to
get the median of a 2D array with axis=1), so you may want to check this one
instead if you feel like it. You can still use 1.1.1 : as a quick workaround
the forementioned bug, use ma.median(a.T, axis=0) instead of
ma.median(a,axis=1) when working w/ 2D arrays.
More information about the Numpy-discussion