【问题标题】:Pandas - group by consecutive rangesPandas - 按连续范围分组
【发布时间】:2016-08-18 13:29:12
【问题描述】:

我有一个具有以下结构的数据框 - 开始、结束和高度。

数据框的一些属性:

  • 数据帧中的一行总是从前一行结束的地方开始,即如果第 n 行的结束是 100,那么第 n+1 行的开始是 101。
  • 第 n+1 行的高度始终与第 n+1 行的高度不同(这就是数据位于不同行的原因)。

我想对数据框进行分组,将高度分组为 5 个长的桶,即桶是 0、1-5、6-10、11-15 和 >15

请参阅下面的代码示例,我正在寻找的是 group_by_bucket 函数的实现。

我尝试查看其他问题,但无法得到我所寻找的确切答案。

提前致谢!

>>> d = pd.DataFrame([[1,3,5], [4,10,7], [11,17,6], [18,26, 12], [27,30, 15], [31,40,6], [41, 42, 7]], columns=['start','end', 'height'])
>>> d
   start  end  height
0      1    3       8
1      4   10       7
2     11   17       6
3     18   26      12
4     27   30      15
5     31   40       6
6     41   42       7
>>> d_gb = group_by_bucket(d)
>>> d_gb
   start  end height_grouped
0      1   17           6_10
1     18   30          11_15
2     31   42           6_10

【问题讨论】:

  • 这不是骗人的吗:stackoverflow.com/questions/21441259/…?在您的情况下,您想在“高度”上调用 cut 并传递一个范围,例如 pd.cut(d['height'], bins=np.arange(1, d['height'].max()+1, 5))
  • 在您的示例中,height_grouped 值不是唯一的。 d_gb 有三个组,但其中两个是相同的,所以你应该有一个额外的分组标准,你没有提到。

标签: python pandas group-by intervals


【解决方案1】:

您可以使用cutgroupby by cutSeriescumsum 用于生成组并通过aggfirstlast 聚合:

bins = [-1,0,1,5,10,15,100]
print bins
[-1, 0, 1, 5, 10, 15, 100]

cut_ser = pd.cut(d['height'], bins=bins)
print cut_ser
0     (5, 10]
1     (5, 10]
2     (5, 10]
3    (10, 15]
4    (10, 15]
5     (5, 10]
6     (5, 10]
Name: height, dtype: category
Categories (6, object): [(-1, 0] < (0, 1] < (1, 5] < (5, 10] < (10, 15] < (15, 100]]

print (cut_ser.shift() != cut_ser).cumsum()
0    0
1    0
2    0
3    1
4    1
5    2
6    2
Name: height, dtype: int32

print d.groupby([(cut_ser.shift() != cut_ser).cumsum(), cut_ser])
       .agg({'start' : 'first','end' : 'last'})
       .reset_index(level=1).reset_index(drop=True)
       .rename(columns={'height':'height_grouped'})

  height_grouped  start  end
0        (5, 10]      1   17
1       (10, 15]     18   30
2        (5, 10]     31   42

编辑:

时间安排

In [307]: %timeit a(df)
100 loops, best of 3: 5.45 ms per loop

In [308]: %timeit b(d)
The slowest run took 4.45 times longer than the fastest. This could mean that an intermediate result is being cached 
100 loops, best of 3: 3.28 ms per loop

代码

d = pd.DataFrame([[1,3,5], [4,10,7], [11,17,6], [18,26, 12], [27,30, 15], [31,40,6], [41, 42, 7]], columns=['start','end', 'height'])
print d

df = d.copy()


def a(df):
    df['groups']=pd.cut(df.height,[-1,0,5,10,15,1000])
    df['categories']=(df.groups!=df.groups.shift()).cumsum()
    f = {'start':['first'],'end':['last'], 'groups':['first']}
    return df.groupby('categories').agg(f)

def b(d):
    bins = [-1,0,1,5,10,15,100]
    cut_ser = pd.cut(d['height'], bins=bins)
    return d.groupby([(cut_ser.shift() != cut_ser).cumsum(), cut_ser]).agg({'start' : 'first','end' : 'last'}).reset_index(level=1).reset_index(drop=True).rename(columns={'height':'height_grouped'})


print a(df)    
print b(d)

【讨论】:

  • 谢谢。这似乎是与 B.M. 类似的解决方案。两者似乎都可以解决问题
  • 是的,看起来很相似,但我的解决方案要快 1.6 倍。
【解决方案2】:

一种方法:

df = pd.DataFrame([[1,3,10], [4,10,7], [11,17,6], [18,26, 12],
[27,30, 15], [31,40,6], [41, 42, 6]], columns=['start','end', 'height'])

使用cut 进行分组:

df['groups']=pd.cut(df.height,[-1,0,5,10,15,1000])

寻找断点:

df['categories']=(df.groups!=df.groups.shift()).cumsum()

那么df就是:

"""
   start  end  height    groups  categories
0      1    3      10   (5, 10]           0
1      4   10       7   (5, 10]           0
2     11   17       6   (5, 10]           0
3     18   26      12  (10, 15]           1
4     27   30      15  (10, 15]           1
5     31   40       6   (5, 10]           2
6     41   42       6   (5, 10]           2
"""

定义有趣的数据:

f = {'start':['first'],'end':['last'], 'groups':['first']}

并使用groupby.agg 函数:

df.groupby('categories').agg(f)
"""
              groups  end start
               first last first
categories                     
0            (5, 10]   17     1
1           (10, 15]   30    18
2            (5, 10]   42    31
"""

【讨论】:

    猜你喜欢
    • 2017-01-30
    • 1970-01-01
    • 2021-10-02
    • 1970-01-01
    • 2011-08-05
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2022-06-29
    相关资源
    最近更新 更多