Re: [Numpy-discussion] Using matplotlib's prctile on masked arrays

2009-10-28 Thread josef . pktd
On Wed, Oct 28, 2009 at 9:52 AM, Gökhan Sever wrote: > > > On Tue, Oct 27, 2009 at 12:23 PM, Pierre GM wrote: >> >> On Oct 27, 2009, at 7:56 AM, Gökhan Sever wrote: >> > >> > >> > Unfortunately, matplotlib.mlab's prctile cannot handle this division: >> >> Actually, the division's OK, it's mlab.pr

Re: [Numpy-discussion] Using matplotlib's prctile on masked arrays

2009-10-28 Thread Gökhan Sever
On Tue, Oct 27, 2009 at 12:23 PM, Pierre GM wrote: > > On Oct 27, 2009, at 7:56 AM, Gökhan Sever wrote: > > > > > > Unfortunately, matplotlib.mlab's prctile cannot handle this division: > > Actually, the division's OK, it's mlab.prctile which is borked. It > uses the length of the input array ins

Re: [Numpy-discussion] Using matplotlib's prctile on masked arrays

2009-10-28 Thread Gökhan Sever
On Tue, Oct 27, 2009 at 8:25 AM, wrote: > This should not be the correct results if you use > scipy.stats.scoreatpercentile, > it doesn't have correct missing value handling, it treats nans or > mask/fill values as regular numbers sorted to the end. > > stats.mstats.scoreatpercentile is the corr

Re: [Numpy-discussion] Using matplotlib's prctile on masked arrays

2009-10-27 Thread Pierre GM
On Oct 27, 2009, at 7:56 AM, Gökhan Sever wrote: > > > Unfortunately, matplotlib.mlab's prctile cannot handle this division: Actually, the division's OK, it's mlab.prctile which is borked. It uses the length of the input array instead of its count to compute the nb of valid data. The easiest

Re: [Numpy-discussion] Using matplotlib's prctile on masked arrays

2009-10-27 Thread josef . pktd
On Tue, Oct 27, 2009 at 7:56 AM, Gökhan Sever wrote: > Hello, > > Consider this sample two columns of data: > >  99. 99. >  99. 99. >  99. 99. >  99.   1693.9069 >  99.   1676.1059 >  99.   1621.5875 >     651.8040   1542.

[Numpy-discussion] Using matplotlib's prctile on masked arrays

2009-10-27 Thread Gökhan Sever
Hello, Consider this sample two columns of data: 99. 99. 99. 99. 99. 99. 99. 1693.9069 99. 1676.1059 99. 1621.5875 651.8040 1542.1373 691.0138 1650.4214 678.5558 1710.7311 621.577