[SciPy-dev] RFR: Proposed fixes in scipy.stats functions for calculation of variance/error/etc.
Sun Oct 25 22:49:16 CDT 2009
Hi Josef and all,
thank for looking. Concerning the z-score functions - I am also
confused by those and I would suggest unifying them under one
function. In particular, I can't imagine what the function 'z' is for.
However, I don't want to just remove these without discussion. What do
you think about this?
Another, more general thing, concerning the axis - I am wondering: why
is the default axis for scipy is 0, while the default for numpy (in
np.mean, for example) is None? I think that it would be good to have
one convention for both libraries. I think that the more parsimonious
one is the one using "None" as the default value. This doesn't favor
any of the dimensions of an array over others, by default. I don't
know - how wide-spread is this convention within scipy?
On Sun, Oct 25, 2009 at 8:16 PM, <email@example.com> wrote:
> On Sun, Oct 25, 2009 at 10:50 PM, Ariel Rokem <firstname.lastname@example.org> wrote:
>> Hi everyone,
>> I have been working on some fixes to the functions in scipy.stats
>> which calculate variance/error and related quantities. In particular,
>> in order to comply with the deprecation warnings that appear in use of
>> scipy.stats.samplevar/scipy.stats.samplestd, I have replaced use of
>> these functions with calls to np.std/np.var. I have also cleaned up
>> the documentation a bit.
>> This can all be found here: http://codereview.appspot.com/141051
> I just gave it a quick look, looks good so far
> in def zs looks like a shape error for axis>0
> "return (a-mu)/sigma"
> def zs changes definition, before it normalized with raveled mean,
> std not by axis
> - mu = np.mean(a,None)
> - sigma = samplestd(a)
> - return (array(a)-mu)/sigma
> + a,axis = _chk_asarray(a,axis)
> + mu = np.mean(a,axis)
> + sigma = np.std(a,axis)
> + return (a-mu)/sigma
> I never looked closely at these,
> zmap has a description I don't understand.
> z, zs, zm ???
> Which is which? they look a bit inconsistent, population might refer
> to dof correction in z ?
> Is there a standard terminology for z scores?
> I think for axis, I have seen more "int or None" ?
>> Ariel Rokem
>> Helen Wills Neuroscience Institute
>> University of California, Berkeley
>> Scipy-dev mailing list
> Scipy-dev mailing list
Helen Wills Neuroscience Institute
University of California, Berkeley
More information about the Scipy-dev