AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • 主页
  • 系统&网络
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • 主页
  • 系统&网络
    • 最新
    • 热门
    • 标签
  • Ubuntu
    • 最新
    • 热门
    • 标签
  • Unix
    • 最新
    • 标签
  • DBA
    • 最新
    • 标签
  • Computer
    • 最新
    • 标签
  • Coding
    • 最新
    • 标签
主页 / coding / 问题 / 79285449
Accepted
amance
amance
Asked: 2024-12-17 00:55:33 +0800 CST2024-12-17 00:55:33 +0800 CST 2024-12-17 00:55:33 +0800 CST

使用 groupby 变换查找特定年份每组的平均比率

  • 772

我正在尝试寻找一种更好/更快的方法来做到这一点。我有一个相当大的数据集(约 2 亿行),每行都有单独的日期。我想找出 2018 年至 2019 年每组的平均年利率。我知道我可以用结果创建一个小的 df 并将其合并回去,但我试图找到一种使用转换的方法。不确定合并是否会更快。一行代码可以加分。

示例数据

rng = np.random.default_rng(seed=123)
df = pd.DataFrame({'group':rng.choice(list('ABCD'), 100),
                   'date':[(pd.to_datetime('2018')+pd.Timedelta(days=x)).normalize() for x in rng.integers(0, 365*5, 100)],
                   'foo':rng.integers(1, 100, 100),
                   'bar':rng.integers(50, 200, 100)})
df['year'] = df['date'].dt.year

这有效

#find average 2018 and 2019 'foo' and 'bar'
for col in ['foo', 'bar']:
    for y in [2018, 2019]:
        df[col+'_'+str(y)+'_total'] = df.groupby('group')['year'].transform(lambda x: df.loc[x.where(x==y).dropna().index, col].sum())

#find 2018 and 2019 rates
for y in [2018, 2019]:
    df['rate_'+str(y)] =  df['foo_'+str(y)+'_total'].div(df['bar_'+str(y)+'_total'])

#find average rate
df['2018_2019_avg_rate'] = df[['rate_2018', 'rate_2019']].mean(axis=1)

我尝试过但不太管用的东西(我apply在切换到之前用它来测试它是否有效transform)

#gives yearly totals for each year and each column, but further 'apply'ing to find rates then averaging isn't working after I switch to transform
df.groupby(['group', 'year'])['year'].apply(lambda x: df.loc[x.where(x.between(2018, 2019)).dropna().index, ['foo', 'bar']].sum())

#close but is averaging too early
df.groupby(['group', 'year'])['year'].apply(lambda x: df.loc[i, 'foo'].sum()/denom if (denom:=df.loc[i:=x.where(x.between(2018, 2019)).dropna().index, 'bar'].sum())>0 else np.nan)
python
  • 2 2 个回答
  • 74 Views

2 个回答

  • Voted
  1. Best Answer
    mozway
    2024-12-17T01:57:31+08:002024-12-17T01:57:31+08:00

    您无法使用 有效地执行多个过滤/聚合groupby.transform。您必须循环。

    更有效的方法是组合pivot_table+ merge:

    cols = ['foo', 'bar']
    years = [2018, 2019]
    
    tmp = (df[df['year'].isin(years)]
           .pivot_table(index='group', columns='year',
                        values=cols, aggfunc='sum')
           [cols]
           .pipe(lambda x: x.join(pd.concat({'rate': x['foo'].div(x['bar'])}, axis=1)))
          )
    
    avg_rate = tmp['rate'].mean(axis=1)
    
    tmp.columns = tmp.columns.map(lambda x: f'{x[0]}_{x[1]}_total')
    
    tmp[f'{"_".join(map(str, years))}_avg_rate'] = avg_rate
    
    out = df.merge(tmp, left_on='group', right_index=True)
    

    输出:

       group       date  foo  bar  year  foo_2018_total  foo_2019_total  bar_2018_total  bar_2019_total  rate_2018_total  rate_2019_total  2018_2019_avg_rate
    0      A 2022-03-11   59   91  2022             343             270             972             875         0.352881         0.308571            0.330726
    1      C 2018-08-22   56   52  2018             175             325             331             902         0.528701         0.360310            0.444506
    2      C 2019-04-24   47   89  2019             175             325             331             902         0.528701         0.360310            0.444506
    3      A 2019-04-16   43  102  2019             343             270             972             875         0.352881         0.308571            0.330726
    4      D 2019-11-25    3   56  2019             126             222             224             696         0.562500         0.318966            0.440733
    5      A 2018-01-06   86  148  2018             343             270             972             875         0.352881         0.308571            0.330726
    ...
    99     B 2018-02-25   32   90  2018             253             204             703             400         0.359886         0.510000            0.434943
    
    • 2
  2. Ouroboroski
    2024-12-17T17:15:42+08:002024-12-17T17:15:42+08:00

    使用变换的方法:

    mask_year = df['year'].isin([2018, 2019])
    
    df.loc[mask_year, 'foo_total'] = df[mask_year].groupby(['group', 'year'])['foo'].transform('sum')
    df.loc[mask_year, 'bar_total'] = df[mask_year].groupby(['group', 'year'])['bar'].transform('sum')
    
    df['rate'] = df['foo_total'] / df['bar_total']
    df.loc[:, '2018_2019_avg_rate'] = df.groupby(['group'])['rate'].transform(lambda x: np.nanmean(x.unique()))
    

    从我有限的测试来看,这比@mozway 的答案稍微慢一些。

    对于 100 万行样本数据,在我的笔记本电脑上运行 100 次的平均时间:

    • OP 的方法:0.76 秒
    • Mozway 的方法:0.23 秒
    • 这种方法:0.29 秒
    • 1

相关问题

  • 如何将 for 循环拆分为 3 个单独的数据框?

  • 如何检查 Pandas DataFrame 中的所有浮点列是否近似相等或接近

  • “load_dataset”如何工作,因为它没有检测示例文件?

  • 为什么 pandas.eval() 字符串比较返回 False

  • Python tkinter/ ttkboostrap dateentry 在只读状态下不起作用

Sidebar

Stats

  • 问题 205573
  • 回答 270741
  • 最佳答案 135370
  • 用户 68524
  • 热门
  • 回答
  • Marko Smith

    Vue 3:创建时出错“预期标识符但发现‘导入’”[重复]

    • 1 个回答
  • Marko Smith

    为什么这个简单而小的 Java 代码在所有 Graal JVM 上的运行速度都快 30 倍,但在任何 Oracle JVM 上却不行?

    • 1 个回答
  • Marko Smith

    具有指定基础类型但没有枚举器的“枚举类”的用途是什么?

    • 1 个回答
  • Marko Smith

    如何修复未手动导入的模块的 MODULE_NOT_FOUND 错误?

    • 6 个回答
  • Marko Smith

    `(表达式,左值) = 右值` 在 C 或 C++ 中是有效的赋值吗?为什么有些编译器会接受/拒绝它?

    • 3 个回答
  • Marko Smith

    何时应使用 std::inplace_vector 而不是 std::vector?

    • 3 个回答
  • Marko Smith

    在 C++ 中,一个不执行任何操作的空程序需要 204KB 的堆,但在 C 中则不需要

    • 1 个回答
  • Marko Smith

    PowerBI 目前与 BigQuery 不兼容:Simba 驱动程序与 Windows 更新有关

    • 2 个回答
  • Marko Smith

    AdMob:MobileAds.initialize() - 对于某些设备,“java.lang.Integer 无法转换为 java.lang.String”

    • 1 个回答
  • Marko Smith

    我正在尝试仅使用海龟随机和数学模块来制作吃豆人游戏

    • 1 个回答
  • Martin Hope
    Aleksandr Dubinsky 为什么 InetAddress 上的 switch 模式匹配会失败,并出现“未涵盖所有可能的输入值”? 2024-12-23 06:56:21 +0800 CST
  • Martin Hope
    Phillip Borge 为什么这个简单而小的 Java 代码在所有 Graal JVM 上的运行速度都快 30 倍,但在任何 Oracle JVM 上却不行? 2024-12-12 20:46:46 +0800 CST
  • Martin Hope
    Oodini 具有指定基础类型但没有枚举器的“枚举类”的用途是什么? 2024-12-12 06:27:11 +0800 CST
  • Martin Hope
    sleeptightAnsiC `(表达式,左值) = 右值` 在 C 或 C++ 中是有效的赋值吗?为什么有些编译器会接受/拒绝它? 2024-11-09 07:18:53 +0800 CST
  • Martin Hope
    The Mad Gamer 何时应使用 std::inplace_vector 而不是 std::vector? 2024-10-29 23:01:00 +0800 CST
  • Martin Hope
    Chad Feller 在 5.2 版中,bash 条件语句中的 [[ .. ]] 中的分号现在是可选的吗? 2024-10-21 05:50:33 +0800 CST
  • Martin Hope
    Wrench 为什么双破折号 (--) 会导致此 MariaDB 子句评估为 true? 2024-05-05 13:37:20 +0800 CST
  • Martin Hope
    Waket Zheng 为什么 `dict(id=1, **{'id': 2})` 有时会引发 `KeyError: 'id'` 而不是 TypeError? 2024-05-04 14:19:19 +0800 CST
  • Martin Hope
    user924 AdMob:MobileAds.initialize() - 对于某些设备,“java.lang.Integer 无法转换为 java.lang.String” 2024-03-20 03:12:31 +0800 CST
  • Martin Hope
    MarkB 为什么 GCC 生成有条件执行 SIMD 实现的代码? 2024-02-17 06:17:14 +0800 CST

热门标签

python javascript c++ c# java typescript sql reactjs html

Explore

  • 主页
  • 问题
    • 最新
    • 热门
  • 标签
  • 帮助

Footer

AskOverflow.Dev

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve