【发布时间】:2015-08-21 20:08:42
【问题描述】:
我有一个数据框列表,每个数据框对应一个日期,每个数据框的行对应一天中的每小时时段。我需要在每天的等效时间戳上应用滚动功能。
例如,对于滚动 5 个周期的窗口,我希望将一个函数应用于每个 01:00:00、02:00:00 等,前 5 天的时间戳一直滚动到 end = length(list)-5,每个函数应用都会产生一个值。结果将是每个长度为 end 的时间戳的值向量。
我正在考虑创建一个新的数据框列表,其中每个数据框将对应一个时间戳,而行对应于天,在这种情况下,我想知道将我当前的数据框列表转换为此的最佳方法新格式。但是,如果有其他更有效的方法,我愿意接受其他建议。
编辑: 我尝试在下面包含示例数据集和结果,但我不知道如何正确格式化数据
例如,如果一个列表包含 4 个具有以下格式的数据框。 Col1 是一个变量,Col2 是另一个变量
Date Time Col1 Col2
1 2014-08-12 09:00:00 -0.0001079389 0.0002159128
2 2014-08-12 09:30:00 -0.0002158895 0.0002159361
3 2014-08-12 10:00:00 -0.0001079739 0.0008642109
4 2014-08-12 10:30:00 -0.0003241666 0.0003242717
5 2014-08-12 11:00:00 -0.0006482282 0.0005405113
6 2014-08-12 11:30:00 -0.0010794474 0.0002160294
7 2014-08-12 12:00:00 -0.0002157963 0.0002158429
8 2014-08-12 12:30:00 -0.0001079040 0.0009716600
9 2014-08-12 13:00:00 -0.0005397528 0.0002159827
10 2014-08-12 13:30:00 0.0000000000 0.0005398111
11 2014-08-12 14:00:00 0.0000000000 0.0003238517
12 2014-08-12 14:30:00 -0.0002158895 0.0003239216
13 2014-08-12 15:00:00 -0.0002159128 0.0001079739
14 2014-08-12 15:30:00 0.0000000000 0.0008642109
15 2014-08-12 16:00:00 -0.0003240966 0.0004322922
16 2014-08-12 16:30:00 -0.0003241666 0.0003242717
17 2014-08-12 17:00:00 -0.0008642109 0.0001080789
18 2014-08-12 17:30:00 -0.0001079739 0.0004320121
19 2014-08-12 18:00:00 -0.0005396945 0.0001079739
如果函数的滚动窗口为 2,则使用乘法等简单函数,将前两个数据帧中的等效行相乘,向前滚动直到第 3 和第 4 个数据帧中的等效行相乘。
Res1 = data.frame(matrix(NA,ncol=3,nrow=19))
rownames(Res1) = Example[[1]][,2]
colnames(Res1) = c("2014-08-13","2014-08-14","2014-08-15")
Res1[,1]= Example[[1]][,3]*Example[[2]][,3]
Res1[,2]= Example[[2]][,3]*Example[[3]][,3]
Res1[,3]= Example[[3]][,3]*Example[[4]][,3]
这是一个简单的示例,因为我要应用的函数将一次滚动向前滚动 5 个数据帧中的等效行。下面,该函数仅应用于 Col1 变量以进行简化,但我将分别对所有变量执行此操作。
结果应如下所示,其中第一列名称是每个滚动窗口的结尾,因此 4 个数据框,相当于 4 个日期变成 3 列
row.names 2014-08-13 2014-08-14 2014-08-15
1 09:00:00 0.000000e+00 0.000000e+00 4.612407e-08
2 09:30:00 6.987091e-08 0.000000e+00 0.000000e+00
3 10:00:00 6.987846e-08 6.954805e-08 4.614140e-08
4 10:30:00 4.889919e-07 1.621737e-07 2.307567e-08
5 11:00:00 0.000000e+00 0.000000e+00 0.000000e+00
6 11:30:00 6.977682e-07 4.869397e-07 8.083873e-08
7 12:00:00 2.325015e-08 5.794719e-08 1.154280e-07
8 12:30:00 0.000000e+00 0.000000e+00 0.000000e+00
9 13:00:00 5.816615e-08 0.000000e+00 0.000000e+00
10 13:30:00 0.000000e+00 9.270554e-08 6.921964e-08
11 14:00:00 0.000000e+00 0.000000e+00 0.000000e+00
12 14:30:00 0.000000e+00 0.000000e+00 9.224832e-08
13 15:00:00 1.627863e-07 1.621387e-07 2.305836e-08
14 15:30:00 0.000000e+00 0.000000e+00 0.000000e+00
15 16:00:00 6.972068e-08 0.000000e+00 0.000000e+00
16 16:30:00 1.045980e-07 3.122420e-07 4.148280e-07
17 17:00:00 0.000000e+00 0.000000e+00 2.419698e-07
18 17:30:00 0.000000e+00 0.000000e+00 2.303854e-07
19 18:00:00 5.810978e-08 1.040920e-07 6.218404e-07
【问题讨论】:
-
您可能不需要数据帧列表,通过使用任何聚合函数(即聚合、data.table、dplyr)将
01:00:00、2:00:00作为新的分组变量可能更容易,并且在组内,可以使用来自library(zoo)的rollapply。最好提供一个小的可重现示例,该示例具有基于该示例的预期输出以进行测试。 -
请提供可重现的示例和所需的输出。
-
@akrun 感谢您的建议,如果有帮助,我现在提供了一个示例
-
@DavidArenburg 我试图包含一个示例和所需的输出,我希望这会让它更清楚