Numpy mean and std over every terms of arrays

Question

I have a list of 2 dimensional arrays (same shape), and would like to get the mean and deviation for all terms, in a result array of the same shape as the inputs. I have trouble understanding from the doc whether this is possible. All my attempts with axis and keepdims parameters produce results of different shapes.

I would like for example to have: mean([x, x]) equal to x, and std([x, x]) zeroes shaped like x.

Is this possible without reshaping the arrays ? If not, how to do it with reshaping ?

Example:

>> x= np.array([[1,2],[3,4]])
>>> y= np.array([[2,3],[4,5]])
>>> np.mean([x,y])
3.0

I want [[1.5,2.5],[3.5,4.5]] instead.

Answer 1

正如Divikar指出的那样，您可以将数组列表传递给，np.mean并指定axis=0对列表中每个数组的对应值求平均值：

In [13]: np.mean([x,y], axis=0)
Out[13]: 
array([[ 1.5,  2.5],
       [ 3.5,  4.5]])

这适用于任意长度的列表。对于仅两个数组，(x+y)/2.0速度更快：

In [20]: %timeit (x+y)/2.0
100000 loops, best of 3: 1.96 µs per loop

In [21]: %timeit np.mean([x,y], axis=0)
10000 loops, best of 3: 21.6 µs per loop