在另一个向量中使用分组值的平均值(numpy / Python)

时间:2015-03-24 21:49:37

标签: python numpy

我想基于另一个向量中的分组信息来取一个向量的平均值。两个向量长度相同。我根据每个用户的平均预测创建了一个最小的例子。我如何在NumPy中做到这一点?

       >>> pred
           [ 0.99  0.23  0.11  0.64  0.45  0.55 0.76  0.72  0.97 ] 
       >>> users
           ['User2' 'User3' 'User2' 'User3' 'User0' 'User1' 'User4' 'User4' 'User4']

3 个答案:

答案 0 :(得分:2)

'纯粹的numpy'解决方案可能会使用np.uniquenp.bincount的组合:

import numpy as np

pred = [0.99,  0.23,  0.11,  0.64,  0.45,  0.55, 0.76,  0.72,  0.97]
users = ['User2', 'User3', 'User2', 'User3', 'User0', 'User1', 'User4',
         'User4', 'User4']

# assign integer indices to each unique user name, and get the total
# number of occurrences for each name
unames, idx, counts = np.unique(users, return_inverse=True, return_counts=True)

# now sum the values of pred corresponding to each index value
sum_pred = np.bincount(idx, weights=pred)

# finally, divide by the number of occurrences for each user name
mean_pred = sum_pred / counts

print(unames)
# ['User0' 'User1' 'User2' 'User3' 'User4']

print(mean_pred)
# [ 0.45        0.55        0.55        0.435       0.81666667]

如果您安装了pandas,则DataFrame拥有some very nice methods for grouping and summarizing data

import pandas as pd

df = pd.DataFrame({'name':users, 'pred':pred})

print(df.groupby('name').mean())
#            pred
# name           
# User0  0.450000
# User1  0.550000
# User2  0.550000
# User3  0.435000
# User4  0.816667

答案 1 :(得分:1)

如果您想坚持使用numpy,最简单的方法是使用np.uniquenp.bincount

>>> pred = np.array([0.99, 0.23, 0.11, 0.64, 0.45, 0.55, 0.76, 0.72, 0.97])
>>> users = np.array(['User2', 'User3', 'User2', 'User3', 'User0', 'User1',
...                   'User4', 'User4', 'User4'])
>>> unq, idx, cnt = np.unique(users, return_inverse=True, return_counts=True)
>>> avg = np.bincount(idx, weights=pred) / cnt
>>> unq
array(['User0', 'User1', 'User2', 'User3', 'User4'],
      dtype='|S5')
>>> avg
array([ 0.45      ,  0.55      ,  0.55      ,  0.435     ,  0.81666667])

答案 2 :(得分:0)

一个紧凑的解决方案是使用numpy_indexed(声明:我是它的作者),它实现了类似于Jaime提出的矢量化解决方案的解决方案;但是使用更清晰的界面和更多测试:

import numpy_indexed as npi
npi.group_by(users).mean(pred)