On Thursday 18 September 2008 10:59:12 Peter Saffrey wrote: > I had looked at masked arrays, but couldn't quite get them to work.
That's unfortunate. > >>> from numeric import * Mmh, typo? > >>> from pylab import rand > >>> a = rand(10,3) > >>> a[a > 0.8] = nan > >>> m = ma.masked_array(a, isnan(a)) > >>> m Another way would be m = ma.masked_where(a>0.8,a) > Remember I want medians of each triple, so I need to median the > transposed matrix: > >>> median(m.T) Ohoh. What version of numpy are you using ? if you don't give an axis parameter, you should get the median of the flattened array, therefore a scalar, not an array. Anyway: you should use ma.median for masked arrays. Else, you're just keeping the NaNs where they were. _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion
