2016-07-10 31 views
-1

滞后计算我有几个数据集(以下简称实施例)加权中的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 

我需要由电流2部分值通过的N个先前总和除以加权每个partOne值第2部分值。

因此,对于上面的短的例子(其具有在每个阵列中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. 等等
+0

尝试发布一个[再现的](http://stackoverflow.com/questions/5963269/how-to-make-a-伟大的 - 可重现的例子)的例子。 – user2100721

回答

2

尝试此

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

数据

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