[CHANGESET]: Statistics function incorrectly computing median

[Ben Abbott]

On Mar 5, 2008, at 4:50 PM, John W. Eaton wrote:

> On 28-Feb-2008, Ben Abbott wrote:
>
> | changeset is attached.
>
> | +2008-02-28  Ben Abbott <bpabbott at mac.com>
> | +
> | +	* statistics/base/statistics.m: Modified to calculate median and
> | +	  quantiles in a manner consistent with method #7 used by GNU's R.
> | +	* statistics/base/__quantile__.m: New function.
> | +	* statistics/base/quantile.m: New function. Matlab compatible.
> | +	* statistics/base/prctile.m: New function. Matlab compatible.
> | +	* miscellaneous/dimfunc.m: New function. Operate on a specific
> | +	  dimension of an N-d array.
>
> The part of this patch that I'm not sure about is dimfunc.  Is that
> really necessary?  If I understand the way it works, it seems that it
> will be really slow to have nested loops and calling a function
> repeatedly instead of working on the full array.  Is there no way to
> avoid this using permute/ipermute to rearrange the data before/after
> processing?
>
> jwe

I spent a day looking at various approaches. I don't recall the detail  
at this point.

However, if you have some idea for how this might be improved, I'm  
happy to work on it ... I'll start by looking over permute (which I  
may have done already have forgotten).

Ben