Pandas“Group By”查询HDFStore中的大数据?

时间:2013-04-03 21:11:04

标签: python pandas pytables

我在HDFStore中有大约700万行,超过60列。数据超出了我的记忆能力。我希望根据列“A”的值将数据聚合成组。 pandas splitting/aggregating/combining的文档假定我已将所有数据都放在DataFrame中,但我无法将整个商店读入内存中DataFrame。在HDFStore中分组数据的正确方法是什么?

1 个答案:

答案 0 :(得分:18)

这是一个完整的例子。

import numpy as np
import pandas as pd
import os

fname = 'groupby.h5'

# create a frame
df = pd.DataFrame({'A': ['foo', 'foo', 'foo', 'foo',
                         'bar', 'bar', 'bar', 'bar',
                         'foo', 'foo', 'foo'],
                   'B': ['one', 'one', 'one', 'two',
                         'one', 'one', 'one', 'two',
                         'two', 'two', 'one'],
                   'C': ['dull', 'dull', 'shiny', 'dull',
                         'dull', 'shiny', 'shiny', 'dull',
                         'shiny', 'shiny', 'shiny'],
                   'D': np.random.randn(11),
                   'E': np.random.randn(11),
                   'F': np.random.randn(11)})


# create the store and append, using data_columns where I possibily
# could aggregate
with pd.get_store(fname) as store:
    store.append('df',df,data_columns=['A','B','C'])
    print "store:\n%s" % store

    print "\ndf:\n%s" % store['df']

    # get the groups
    groups = store.select_column('df','A').unique()
    print "\ngroups:%s" % groups

    # iterate over the groups and apply my operations
    l = []
    for g in groups:

        grp = store.select('df',where = [ 'A=%s' % g ])

        # this is a regular frame, aggregate however you would like
        l.append(grp[['D','E','F']].sum())


    print "\nresult:\n%s" % pd.concat(l, keys = groups)

os.remove(fname)

输出

store:
<class 'pandas.io.pytables.HDFStore'>
File path: groupby.h5
/df            frame_table  (typ->appendable,nrows->11,ncols->6,indexers->[index],dc->[A,B,C])

df:
      A    B      C         D         E         F
0   foo  one   dull -0.815212 -1.195488 -1.346980
1   foo  one   dull -1.111686 -1.814385 -0.974327
2   foo  one  shiny -1.069152 -1.926265  0.360318
3   foo  two   dull -0.472180  0.698369 -1.007010
4   bar  one   dull  1.329867  0.709621  1.877898
5   bar  one  shiny -0.962906  0.489594 -0.663068
6   bar  one  shiny -0.657922 -0.377705  0.065790
7   bar  two   dull -0.172245  1.694245  1.374189
8   foo  two  shiny -0.780877 -2.334895 -2.747404
9   foo  two  shiny -0.257413  0.577804 -0.159316
10  foo  one  shiny  0.737597  1.979373 -0.236070

groups:Index([bar, foo], dtype=object)

result:
bar  D   -0.463206
     E    2.515754
     F    2.654810
foo  D   -3.768923
     E   -4.015488
     F   -6.110789
dtype: float64

一些警告:

1)如果你的团体密度相对较低,这种方法是有意义的。按数百或数千个组的顺序。如果你得到更多,那么效率更高(但更复杂的方法),你正在应用的函数(在这种情况下sum)变得更加严格。

基本上你会在整个商店中通过块进行迭代,在你去的时候进行分组,但保持组只是半崩溃(想象做一个意思,所以你需要保持一个运行总数加上一个运行计数,然后除以结束)。所以一些操作会有点棘手,但可能会处理许多组(并且非常快)。

2)通过保存坐标(例如组位置,但这有点复杂)可以提高效率。

3)这种方案不可能实现多分组(这是可能的,但需要一种更像2的方法)

4)要分组的列,必须是data_column!

5)你可以在select btw中组合你想要的任何其他过滤器(这是一种多组合btw的一种偷偷摸摸的方式,你只需要在它们的产品上形成2个唯一的组和迭代器列表,如果不是非常有效的话你有很多小组,但可以工作)

HTH

让我知道这是否适合您