R中的加权滞后计算

Question

我有几个数据集（下面的简短示例）

> print(partOne)
 [1]  0.010966943 -0.006819522 -0.007189830  0.039736714  0.002013070
 [6] -0.043946666  0.003808415  0.199010991 -0.094001478 -0.053006526
[11] -0.051489992  0.019122283 -0.011215761  0.057408738 -0.020809564
[16] -0.041295915  0.010134001 -0.011398076
> print(part2)
 [1] 0.13070012 0.15793754 0.06980192 0.13270089 0.11384798 0.24417631
 [7] 0.10363273 0.09182783 0.12217490 0.47649356 0.33660955 0.23079863
[13] 0.21581061 0.13967763 0.05988797 0.28255164 0.16277804 0.12716975
[19] 0.19299641 0.21452418

我需要通过当前 part2 值除以前N个 part2 值的总和来对每个 partOne 值进行加权。

因此，对于上面的简短示例（每个数组中有20个值），伪代码将是：

1. 跳至项目N + 1（例如，假设此示例为N = 10）
2. 计算（partOne [11] * partTwo [11]）/ sum（partTwo [1-＆gt; 10]）
3. 增加+1
4. 计算（partOne [12] * partTwo [12]）/ sum（partTwo [2-＆gt; 11]）
5. 等。等

Answer 1

试试这个

x[11:20]*y[11:20]/sapply(1:10,function(t) sum(y[t:(t+9)]))

数据

x <- rnorm(20)
y <- rnorm(20)