R中的加权滞后计算

时间:2016-07-10 17:10:32

标签: r

我有几个数据集(下面的简短示例)

> 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-> 10])
  3. 增加+1
  4. 计算(partOne [12] * partTwo [12])/ sum(partTwo [2-> 11])
  5. 等。等

1 个答案:

答案 0 :(得分:2)

试试这个

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

数据

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