Roll применяется для неполных временных рядов в R

Данный

z <- zoo(c(1:10))

Я хочу иметь возможность агрегировать следующее:

> z
 4  8  10 
 10 26 19

Я пробовал следующее, используя rollapply, но безрезультатно:

> rollapply(zoo(c(1:10)), width = 4, FUN = "sum", by = 4, partial = TRUE, align = "right")
 1  5  9 
 1 14 30 
> rollapply(zoo(c(1:10)), width = 4, FUN = "sum", by = 4, partial = TRUE, align = "left")
 1  5  9 
10 26 19 
> rollapply(zoo(c(1:10)), width = 4, FUN = "sum", by = 4, partial = TRUE, align = "center")
 1  5  9 
 6 22 27

Любая помощь будет принята с благодарностью. Второй выглядит наиболее многообещающе, но мне придется настроить лаг?


r rollapply
person bob    schedule 17.05.2015    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Аргумент partial всегда применяется к обоим концам; однако можно указать ширину для каждого элемента отдельно, используя вектор для аргумента width, а затем подмножить его самостоятельно вместо использования by:

library(zoo)

# inputs
z <- zoo(1:10)
k <- 4    

n <- length(z)
w <- rep(1:k, length = n)  # 1 2 3 4 1 2 3 4 1 2 
ok <- w == k | seq(n) == n  # F F F T F F F T F T

rollapplyr(z, w, sum)[ok]

давая:

 4  8 10 
10 26 19

2) Мы могли бы использовать align = "left", а затем исправить время (используя ok сверху):

r <- rollapply(z, k, by = k, sum, partial = TRUE, align = "left")
time(r) <- time(z)[ok]

3) Это можно сделать с помощью aggregate.zoo (используя ok сверху):

tt <- na.locf(replace(time(z), !ok, NA), fromLast = TRUE)  # 4 4 4 4 8 8 8 8 10 10
aggregate(z, tt, sum)
person G. Grothendieck    schedule 17.05.2015