Pandas_分组(三) 暂时停更

分组

import numpy as np
import pandas as pd
df = pd.read_csv('data/table.csv',index_col='ID')
df.head()

v2-4fd975473f6107b2f044c426d9ed6b21_b.jpg

一、SAC过程

1. 内涵

SAC指的是分组操作中的split-apply-combine过程

其中split指基于某一些规则,将数据拆成若干组,apply是指对每一组独立地使用函数,combine指将每一组的结果组合成某一类数据结构

2. apply过程

在该过程中,我们实际往往会遇到四类问题:

  • 整合(Aggregation)——即分组计算统计量(如求均值、求每组元素个数)
  • 变换(Transformation)——即分组对每个单元的数据进行操作(如元素标准化)
  • 过滤(Filtration)——即按照某些规则筛选出一些组(如选出组内某一指标小于50的组)

综合问题——即前面提及的三种问题的混合

二、groupby函数

通过使用 GroupBy 可以做到涉及一个或多个以下步骤的过程:

  • 根据一些标准将数据拆分成组
  • 独立应用功能给每个组
  • 将结果合并成数据结构

1. 分组函数的基本内容:

(a)根据某一列分组

grouped_single = df.groupby('School')

经过groupby后会生成一个groupby对象,该对象本身不会返回任何东西,只有当相应的方法被调用才会起作用

例如取出某一个组:

grouped_single.get_group('S_1').head()

v2-13ece782b6b917b8813315bca7864862_b.jpg

(b)根据某几列分组

grouped_mul = df.groupby(['School','Class'])
grouped_mul.get_group(('S_2','C_4'))

v2-61ea3024275ce628bc773f42a838c8c0_b.jpg

(c)组容量与组数

grouped_single.size()

v2-27da9b2bdc1086fafa1a712159e08178_b.png
grouped_mul.size()

v2-cbf4ef1f2556e407a6534b3a6e89e566_b.png
grouped_single.ngroups

2

grouped_mul.ngroups

7

(d)组的遍历

for name,group in grouped_single:
    print(name)
    display(group.head())

v2-fdf8cd990121d3f9e75ccdb846aa2fd3_b.jpg

(e)level参数(用于多级索引)和axis参数

df.set_index(['Gender','School']).groupby(level=1,axis=0).get_group('S_1').head()

v2-ec007f5a60d5eb33c4ba6063823bba33_b.jpg

2. groupby对象的特点

(a)查看所有可调用的方法

由此可见,groupby对象可以使用相当多的函数,灵活程度很高

print([attr for attr in dir(grouped_single) if not attr.startswith('_')])

['Address', 'Class', 'Gender', 'Height', 'Math', 'Physics', 'School', 'Weight', 'agg', 'aggregate', 'all', 'any', 'apply', 'backfill', 'bfill', 'boxplot', 'corr', 'corrwith', 'count', 'cov', 'cumcount', 'cummax', 'cummin', 'cumprod', 'cumsum', 'describe', 'diff', 'dtypes', 'expanding', 'ffill', 'fillna', 'filter', 'first', 'get_group', 'groups', 'head', 'hist', 'idxmax', 'idxmin', 'indices', 'last', 'mad', 'max', 'mean', 'median', 'min', 'ndim', 'ngroup', 'ngroups', 'nth', 'nunique', 'ohlc', 'pad', 'pct_change', 'pipe', 'plot', 'prod', 'quantile', 'rank', 'resample', 'rolling', 'sem', 'shift', 'size', 'skew', 'std', 'sum', 'tail', 'take', 'transform', 'tshift', 'var']

(b)分组对象的head和first

对分组对象使用head函数,返回的是每个组的前几行,而不是数据集前几行

grouped_single.head(2)

v2-dd8d87a19f48d3b2c3bc3710751fdf08_b.jpg

first显示的是以分组为索引的每组的第一个分组信息

grouped_single.first()

v2-e73c5548f94e7782d80b6c30d0a22ca1_b.jpg

(c)分组依据

对于groupby函数而言,分组的依据是非常自由的,只要是与数据框长度相同的列表即可,同时支持函数型分组

df.groupby(np.random.choice(['a','b','c'],df.shape[0])).get_group('a').head()
#相当于将np.random.choice(['a','b','c'],df.shape[0])当做新的一列进行分组

v2-06387b680e38c6a02ae02a0535e7a0c4_b.jpg

从原理上说,我们可以看到利用函数时,传入的对象就是索引,因此根据这一特性可以做一些复杂的操作

df[:5].groupby(lambda x:print(x)).head(0)

v2-d171fe793e05ceede2eba293492ac916_b.jpg

根据奇偶行分组

df.groupby(lambda x:'奇数行' if not df.index.get_loc(x)%2==1 else '偶数行').groups

v2-bdb946003160061e85439e15d46ffd3f_b.jpg

如果是多层索引,那么lambda表达式中的输入就是元组,下面实现的功能为查看两所学校中男女生分别均分是否及格

注意:此处只是演示groupby的用法,实际操作不会这样写

math_score = df.set_index(['Gender','School'])['Math'].sort_index()
grouped_score = df.set_index(['Gender','School']).sort_index().\
            groupby(lambda x:(x,'均分及格' if math_score[x].mean()>=60 else '均分不及格'))
for name,_ in grouped_score:print(name)

v2-21c1c8eaafb757ca485cd5bedc876076_b.jpg

(d)groupby的[]操作

可以用[]选出groupby对象的某个或者某几个列,上面的均分比较可以如下简洁地写出:

df.groupby(['Gender','School'])['Math'].mean()>=60

v2-a98df79b11642f94272501a0c8d2aa71_b.png

用列表可选出多个属性列:

df.groupby(['Gender','School'])[['Math','Height']].mean()

v2-8a9acb368a3dfaa4a0c5d07a3f136d41_b.jpg

(e)连续型变量分组

例如利用cut函数对数学成绩分组:

bins = [0,40,60,80,90,100]
cuts = pd.cut(df['Math'],bins=bins) #可选label添加自定义标签
df.groupby(cuts)['Math'].count()

v2-e2df2c4dbc225344b6245812db975374_b.png

三、聚合、过滤和变换

1. 聚合(Aggregation)

(a)常用聚合函数

所谓聚合就是把一堆数,变成一个标量,因此mean/sum/size/count/std/var/sem/describe/first/last/nth/min/max都是聚合函数

为了熟悉操作,不妨验证标准误sem函数,它的计算公式是:

\frac{组内标准差}{\sqrt{组容量}} ,下面进行验证:
group_m = grouped_single['Math']
group_m.std().values/np.sqrt(group_m.count().values)== group_m.sem().values

array([ True, True])

(b)同时使用多个聚合函数

group_m.agg(['sum','mean','std'])

v2-504462e6fefa6154095db3c8cd988d58_b.jpg

利用元组进行重命名

group_m.agg([('rename_sum','sum'),('rename_mean','mean')])

v2-49c1ad757f680927fb8fd00e39bb8fa3_b.jpg

指定哪些函数作用哪些列

grouped_mul.agg({'Math':['mean','max'],'Height':'var'})

v2-bb071be9367d97af7f794a8625a395af_b.jpg

(c)使用自定义函数

grouped_single['Math'].agg(lambda x:print(x.head(),'间隔'))
#可以发现,agg函数的传入是分组逐列进行的,有了这个特性就可以做许多事情

v2-24f02b7584461d8d735d6f013d91f474_b.jpg

官方没有提供极差计算的函数,但通过agg可以容易地实现组内极差计算

grouped_single['Math'].agg(lambda x:x.max()-x.min())

v2-b5f0dcf16927ed79d9e8cbb7cad6258f_b.jpg

(d)利用NamedAgg函数进行多个聚合

注意:不支持lambda函数,但是可以使用外置的def函数

def R1(x):
    return x.max()-x.min()
def R2(x):
    return x.max()-x.median()
grouped_single['Math'].agg(min_score1=pd.NamedAgg(column='col1', aggfunc=R1),
                           max_score1=pd.NamedAgg(column='col2', aggfunc='max'),
                           range_score2=pd.NamedAgg(column='col3', aggfunc=R2)).head()

v2-17f5f37c1a08ffa19d28ab26ad24392c_b.jpg

(e)带参数的聚合函数

判断是否组内数学分数至少有一个值在50-52之间:

def f(s,low,high):
    return s.between(low,high).max()
grouped_single['Math'].agg(f,50,52)

v2-69f356a24242cf2da2003134f465d1b1_b.png

如果需要使用多个函数,并且其中至少有一个带参数,则使用wrap技巧:

def f_test(s,low,high):
    return s.between(low,high).max()
def agg_f(f_mul,name,*args,**kwargs):
    def wrapper(x):
        return f_mul(x,*args,**kwargs)
    wrapper.__name__ = name
    return wrapper
new_f = agg_f(f_test,'at_least_one_in_50_52',50,52)
grouped_single['Math'].agg([new_f,'mean']).head()

v2-6f29667766f7417a5901cb9fdaa21c6b_b.jpg

2. 过滤(Filteration)

filter函数是用来筛选某些组的(务必记住结果是组的全体),因此传入的值应当是布尔标量

grouped_single[['Math','Physics']].filter(lambda x:(x['Math']>32).all()).head()

v2-02339e2d630f97d87e21a0ad6c5a04ee_b.png

3. 变换(Transformation)

(a)传入对象

transform函数中传入的对象是组内的列,并且返回值需要与列长完全一致

grouped_single[['Math','Height']].transform(lambda x:x-x.min()).head()

v2-cbcc35aa39021db817af86c6c48696e9_b.png

如果返回了标量值,那么组内的所有元素会被广播为这个值

grouped_single[['Math','Height']].transform(lambda x:x.mean()).head()

v2-2d3d9adc41302d968a11ac1a52272f70_b.jpg

(b)利用变换方法进行组内标准化

grouped_single[['Math','Height']].transform(lambda x:(x-x.mean())/x.std()).head()

v2-00d33a8dfa7af69c9841d2d767c5ba3f_b.png

(c)利用变换方法进行组内缺失值的均值填充

df_nan = df[['Math','School']].copy().reset_index()
df_nan.loc[np.random.randint(0,df.shape[0],25),['Math']]=np.nan
df_nan.head()

v2-9abd1d5bbaf11235c889d5c2cd37bfd7_b.png
df_nan.groupby('School').transform(lambda x: x.fillna(x.mean())).join(df.reset_index()['School']).head()

v2-fa32906a9f1c86f41f8dc6a52fcd5031_b.png

四、apply函数

1. apply函数的灵活性

可能在所有的分组函数中,apply是应用最为广泛的,这得益于它的灵活性:

对于传入值而言,从下面的打印内容可以看到是以分组的表传入apply中:

df.groupby('School').apply(lambda x:print(x.head(1)))

v2-e6336f4b9977d900497f345c1b6076d8_b.jpg

apply函数的灵活性很大程度来源于其返回值的多样性:

① 标量返回值

df[['School','Math','Height']].groupby('School').apply(lambda x:x.max())

v2-3559ffc1072aff3b4bfda745687e2b6b_b.jpg

② 列表返回值

df[['School','Math','Height']].groupby('School').apply(lambda x:x-x.min()).head()

v2-1eb456a4c68850c6f9b29d835ece6496_b.png

③ 数据框返回值

df[['School','Math','Height']].groupby('School')\
    .apply(lambda x:pd.DataFrame({'col1':x['Math']-x['Math'].max(),
                                  'col2':x['Math']-x['Math'].min(),
                                  'col3':x['Height']-x['Height'].max(),
                                  'col4':x['Height']-x['Height'].min()})).head()

v2-7b31c45b7420c354eaf398fa605aa91a_b.jpg

2. 用apply同时统计多个指标

此处可以借助OrderedDict工具进行快捷的统计:

from collections import OrderedDict
def f(df):
    data = OrderedDict()
    data['M_sum'] = df['Math'].sum()
    data['W_var'] = df['Weight'].var()
    data['H_mean'] = df['Height'].mean()
    return pd.Series(data)
grouped_single.apply(f)

v2-65d0d722dd974b0b573893c32852bc1d_b.jpg

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值