aggreate forward instead of backward.

project-aero · Feb 13, 2018 · 026cb56 · 026cb56
1 parent dd35898
commit 026cb56
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/eric-figs/fig1_new_variance.R b/eric-figs/fig1_new_variance.R
@@ -105,8 +105,12 @@ tychoL1.measles.CA.cases.imp.zoo.statswindow <- window(tychoL1.measles.CA.cases.
 tmpfcn <- function(x,agg.interval=1) {
   index <- index(x)
   interval <- mean(diff(index))
-  agg.groups <- rep(seq(index[length(index)], index[1], -interval*agg.interval), each=agg.interval)
-  agg.groups <- rev(agg.groups[1:length(x)])
+  # forwards:
+  agg.groups <- rep(seq(index[1], index[length(index)], interval*agg.interval), each=agg.interval)
+  agg.groups <- agg.groups[1:length(x)]
+  # reverse:
+  # agg.groups <- rep(seq(index[length(index)], index[1], -interval*agg.interval), each=agg.interval)
+  # agg.groups <- rev(agg.groups[1:length(x)])
   out <- aggregate(x, by = agg.groups, sum)
   return(out)
 }