[Numpy-discussion] Medians that ignore values
Peter Saffrey
pzs at dcs.gla.ac.uk
Thu Sep 18 13:31:18 EDT 2008
Pierre GM <pgmdevlist <at> gmail.com> writes:
> Mmh, typo?
>
Yes, apologies. I was aiming for thorough, but ended up just careless. It's been
a long day.
> Ohoh. What version of numpy are you using ?
The version in the Ubuntu package repository. It says 1:1.0.4-6ubuntu3.
> 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.
>>> a = rand(10,3)
>>> a
array([[ 0.1269796 , 0.43003978, 0.4700416 ],
[ 0.28867077, 0.85265318, 0.35908364],
[ 0.72967127, 0.41856028, 0.54724918],
[ 0.28821876, 0.69684144, 0.54647616],
[ 0.09592476, 0.83704808, 0.52425368],
[ 0.743552 , 0.44433314, 0.7362179 ],
[ 0.4283931 , 0.13305385, 0.68422292],
[ 0.68860674, 0.15057373, 0.99206493],
[ 0.31846329, 0.77237046, 0.986883 ],
[ 0.4578616 , 0.4580833 , 0.97754176]])
>>> median(a.T)
array([ 0.43003978, 0.35908364, 0.54724918, 0.54647616, 0.52425368,
0.7362179 , 0.4283931 , 0.68860674, 0.77237046, 0.4580833 ])
> 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.
According to this page:
http://www.scipy.org/Download
1.1.0 is the latest release. Do I need to use an SVN build to get the ma.median
functionality?
Peter
More information about the NumPy-Discussion
mailing list