4 回答
TA贡献1783条经验 获得超4个赞
我不会说这是重复的,但您提到的相关问题是一个很好的起点。您链接的大多数答案都需要对数组进行排序,提取组开始的索引,然后调用np.split它。这里不是这种情况,因为它会返回一个大小不平衡的组列表。
相反,您可以使用np.bincount方法。它计算每个加权值出现的次数,这实际上与 groupby sum 相同,只是输出中缺少组键。
def group_by_sum(x):
u, idx = np.unique(x[:,0], return_inverse=True)
s = np.bincount(idx, weights = x[:,1])
return np.c_[u, s]
奖金。它实际上是numpy_indexed包装中的一个单行:
np.transpose(npi.group_by(x[:, 0]).sum(x[:, 1]))
标杆管理
import numpy as np
import perfplot
import matplotlib.pyplot as plt
def bincount(x):
u, idx = np.unique(x[:,0], return_inverse=True)
s = np.bincount(idx, weights = x[:,1])
return np.c_[u, s]
def reduceat(x):
x = x[np.argsort(x[:, 0])]
i = np.flatnonzero(np.diff(x[:, 0]))
i = np.r_[0, i + 1]
s = np.add.reduceat(x[:, 1], i)
return np.stack((x[i, 0], s), axis=-1)
def setup(N, s):
x = np.linspace(0,1,N+1)[np.random.randint(N, size = s)]
return np.c_[x, (x**2)%1]
def build_args(k):
return {'setup': lambda x: setup(k, x),
'kernels': [bincount, reduceat],
'n_range': [2**k for k in range(1, 20)],
'title': f'Testing for x samples in [0, 1] with no more than {k} groups',
'show_progress': True,
'equality_check': False}
outs = [perfplot.bench(**build_args(n)) for n in (10, 100, 1000, 10000)]
fig = plt.figure(figsize=(20, 20))
for i in range(len(outs)):
ax = fig.add_subplot(2, 2, i + 1)
ax.grid(True, which="both")
outs[i].plot()
plt.show()
TA贡献1777条经验 获得超10个赞
Numpy 提供了无需显式循环即可完成此操作的工具。
首先对行进行排序:
a = a[np.argsort(a[:, 0])]
然后找到值发生变化的索引:
i = np.flatnonzero(np.diff(a[:, 0])) i = np.r_[0, i + 1]
然后将元素相加:
s = np.add.reduceat(a[:, 1], i)
索引只是a
每次运行中的第一个元素,因此结果是
result = np.stack((a[i, 0], s), axis=-1)
TA贡献1860条经验 获得超9个赞
这是一个使用唯一值来计算每个元素的重复次数并将其乘以其值来计算 groupby 总和的解决方案(您可以通过实现仅计算重复和唯一值的哈希图来更快地实现它)O(n):
编辑原始问题已编辑:
keys2, idx, count = np.unique(x[:,0], return_counts=True, return_index=True)
values2 = x[:,1][idx]*count
另一种方法是使用 pandas groupby:
df = pd.DataFrame({'keys':x[:,0], 'values':x[:,1]})
df2 = df.groupby(keys)['values'].agg('sum')
keys2, values2 = df2.index.to_numpy(), df2.values
输出:
[1.2 2.3]
[20 30]
TA贡献1812条经验 获得超5个赞
这是一个方法
d = {}
for k,v in x:
d[k] = d.get(k,0) + v
x = np.array(list(d.items()))
请记住,这是测试浮动相等性......您可能不应该这样做
添加回答
举报