2016-07-10 24 views
-1

Ich habe eine Reihe von Datensätzen (Kurz Beispiele unten)Weighted verzögerte Berechnung in 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 

Ich brauche jede PartOne Wert durch den aktuellen Part2 Wert durch die Summe von N vorhergehenden geteilt gewichten Teil2 Werte.

Also, für das kurzes Beispiel oben (die 20 Werte in jeder Array), ein Pseudo-Code wäre:

  1. Punkte überspringen N + 1 (zB annehmen, N = 10 in diesem Beispiel)
  2. berechnen (PartOne [11] * PartTwo [11])/Summe (PartTwo [1-> 10])
  3. Inkrement +1
  4. berechnen (PartOne [12] * PartTwo [12])/Summe (partTwo [2-> 11])
  5. usw. usw.
+0

Versuchen Sie, eine [reproduzierbaren] zu schreiben (http://stackoverflow.com/questions/5963269/how-to-make-a- groß-r-reproduzierbar-Beispiel). – user2100721

Antwort

2

Versuchen Sie, diese

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

Daten

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