pandas 滚动 window 有效添加新行

pandas rolling window effective add new rows

我有一个相当大的数据集(大约 500 万行),其中包含多个计算列,例如滞后(1 和 7)和滚动 windows(7、30、90,每个列都有多个值,例如平均值、标准差、最小值、最大值等)。现在我需要向 df 添加一行或多行,我想知道(重新)计算这些特征的最有效方法。重新计算整个 df 会花费太多时间,但我不能简单地将函数应用于新添加的行。

对于滞后这不是一个大问题,我可以简单地说,例如对于滞后 1,

df.iloc[-1, -2] = df.iloc[-2, -2]

这应该可以解决问题(滞后 7 也是如此),但是滚动 windows 呢?同样的方法行不通,我负担不起 (re)运行 整个数据帧上的滚动 windows。我应该手工计算吗?还有其他方法吗?

Python 3.7.7 和 pandas 1.0.3

如果我正确地解释了这个问题,那么您有一个大型 DataFrame,其中包含一个或多个源数据列,然后是多个列,其中包含基于源列的 windowed 汇总统计信息。在将新行附加到源数据列后,您正在尝试更新 windowed 摘要列的底部,而没有 re-calculating 整个摘要列。

实现这一点的方式取决于很多因素,包括您是否使用居中 windows。但希望这能让你入门。

我将从你的问题的玩具版本开始,只有一个 source 列和两个 windowed 方法:

In [2]: df = pd.DataFrame({'source': np.arange(0, 20, 2)})

In [3]: for window in [3, 5]:
   ...:     df[f'rolling_mean_{window}'] = (
   ...:         df.source.rolling(window, center=True).mean())
   ...:

然后我们在底部追加一个新行:

In [4]: df = df.append(pd.Series({'source': 100}), ignore_index=True)

In [5]: df
Out[5]:
    source  rolling_mean_3  rolling_mean_5
0      0.0             NaN             NaN
1      2.0             2.0             NaN
2      4.0             4.0             4.0
3      6.0             6.0             6.0
4      8.0             8.0             8.0
5     10.0            10.0            10.0
6     12.0            12.0            12.0
7     14.0            14.0            14.0
8     16.0            16.0             NaN
9     18.0             NaN             NaN
10   100.0             NaN             NaN

我们必须更新的数据量取决于 window 的长度。例如,要更新 rolling_mean_3,我们需要使用最后五行的信息更新最后两行。为了安全起见,我们可以 re-calculate 最后 2*window 行加上您添加的行数:

In [6]: df.source.iloc[-(2*window+1):].rolling(window, center=True).mean()
Out[6]:
4           NaN
5     10.000000
6     12.000000
7     14.000000
8     16.000000
9     44.666667
10          NaN
Name: source, dtype: float64

第 5-10 行的数据正确。请注意,第 4 行在此版本中不正确(现在为 NaN),但我们可以使用此结果仅更新最后 [-(window+1):] 行。这是完整的解决方案:

In [7]: updated_rows = 1
In [8]: for window in [3, 5]:
   ...:     update_column_name = f'rolling_mean_{window}'
   ...:     update_column_index = df.columns.get_loc(update_column_name)
   ...:     df.iloc[-(window+updated_rows):, update_column_index] = (
   ...:         df.source
   ...:         .iloc[-(window*2+updated_rows):]
   ...:         .rolling(window, center=True).mean()
   ...:         .iloc[-(window+updated_rows):]
   ...:     )

In [9]: df
Out[9]:
    source  rolling_mean_3  rolling_mean_5
0      0.0             NaN             NaN
1      2.0        2.000000             NaN
2      4.0        4.000000             4.0
3      6.0        6.000000             6.0
4      8.0        8.000000             8.0
5     10.0       10.000000            10.0
6     12.0       12.000000            12.0
7     14.0       14.000000            14.0
8     16.0       16.000000            32.0
9     18.0       44.666667             NaN
10   100.0             NaN             NaN

现在已更新为具有正确计算的尾巴。

从技术上讲,对于居中滚动操作,您只需要更新最后 floor(window/2)+updated_rows 行,从数据帧的最后 window+updated_rows 行绘制。所以你可以这样做来真正优化事情。

如果您要生成不居中的滚动统计数据,方法应该相同,但不包括居中标志。

您只需要估计最后一组项目的平均值。见下文

updated_rows = 1
for window in [3, 5]:
    update_column_name = f'rolling_mean_{window}'
    update_column_index = df.columns.get_loc(update_column_name)
    df.iloc[-(updated_rows):, update_column_index] = df.source.iloc[-(window):].mean()