我有在每个小区1x35000阵列(或NaN)一个16x10熊猫数据帧。 我想借此逐元素平均超过行对每一列。
1 2 3 ... 10
1 1x35000 1x35000 1x35000 1x35000
2 1x35000 NaN 1x35000 1x35000
3 1x35000 NaN 1x35000 NaN
...
16 1x35000 1x35000 NaN 1x35000
为了避免误解:取每个阵列的第一元件中的第一列和取均值。 然后采取在第一列中的每个阵列的第二元件和再取平均值。 最后,我希望有与每个每列一个1x35000阵列的1×10数据帧。 数组应该是我前阵的元素方面的意思。
1 2 3 ... 10
1 1x35000 1x35000 1x35000 1x35000
你有一个想法,到那里优雅最好不用for循环?
设定
np.random.seed([3,14159])
df = pd.DataFrame(
np.random.randint(10, size=(3, 3, 5)).tolist(),
list('XYZ'), list('ABC')
).applymap(np.array)
df.loc['X', 'B'] = np.nan
df.loc['Z', 'A'] = np.nan
df
A B C
X [4, 8, 1, 1, 9] NaN [8, 2, 8, 4, 9]
Y [4, 3, 4, 1, 5] [1, 2, 6, 2, 7] [7, 1, 1, 7, 8]
Z NaN [9, 3, 8, 7, 7] [2, 6, 3, 1, 9]
解
g = df.stack().groupby(level=1)
g.apply(np.sum, axis=0) / g.size()
A [4.0, 5.5, 2.5, 1.0, 7.0]
B [5.0, 2.5, 7.0, 4.5, 7.0]
C [5.66666666667, 3.0, 4.0, 4.0, 8.66666666667]
dtype: object
如果你坚持你的形状呈现
g = df.stack().groupby(level=1)
(g.apply(np.sum, axis=0) / g.size()).to_frame().T
A B C
0 [4.0, 5.5, 2.5, 1.0, 7.0] [5.0, 2.5, 7.0, 4.5, 7.0] [5.66666666667, 3.0, 4.0, 4.0, 8.66666666667]
方法一:糊涂
鉴于混合D型输入数据,我们可以通过对性能效率要循环。 因此,带有显式循环或在引擎盖的用途循环.apply/.applymap
将可能被建议的解决方案。
这里有一种方法,通过列循环 -
mask = ~df.isnull().values
n = df.shape[1]
out = np.empty((1,n),dtype=object)
for i in range(n):
out[0,i] = df.iloc[mask[:,i],i].mean()
df_out = pd.DataFrame(out)
样品输入,输出 -
In [326]: df
Out[326]:
0 1 2
0 [4, 0, 1, 6] [4, 2, 2, 2] [5, 3, 5, 4]
1 NaN [0, 5, 6, 8] NaN
2 NaN NaN NaN
3 NaN NaN NaN
In [327]: df_out
Out[327]:
0 1 2
0 [4.0, 0.0, 1.0, 6.0] [2.0, 3.5, 4.0, 5.0] [5.0, 3.0, 5.0, 4.0]
方法2:矢量化
如果有向量化,这是一个使用单程matrix-multiplication
来代替mean-reductions
和可能带来改善大数据-
mask = ~df.isnull().values
v = np.vstack(df.values[mask])
r,c = np.where(mask)
n = df.shape[1]
pos_mask = c == np.arange(n)[:,None]
out = pos_mask.dot(v)/np.bincount(c).astype(float)[:,None]
df_out1 = pd.DataFrame(out)
示例输出 -
In [328]: df_out1
Out[328]:
0 1 2 3
0 4.0 0.0 1.0 6.0
1 2.0 3.5 4.0 5.0
2 5.0 3.0 5.0 4.0
方法3:矢量一个更
利用的np.add.reduceat
得到这些mean-reductions
-
mask = ~df.T.isnull().values
v = np.vstack(df.values.T[mask])
count = mask.sum(1)
out0 = np.add.reduceat(v, np.r_[0,count.cumsum()[:-1]])
out = out0/count[:,None].astype(float)
df_out2 = pd.DataFrame(out)