在计算xts加权和时如何避免for循环？

Question

我如何避免循环计算xts加权和，如下所示：

library(xts)

thetaSum <- function(theta, w=c(1, 1, 1)) {
    sum(coredata(theta)*rev(w))
}

n <- 10
tmpVec <- rep(1, n)
tmpDates <- seq(as.Date("2000-01-01"), length = n, by = "day")
theta <- xts(tmpVec, order.by=tmpDates)

N <- 3
thetaSummed <- xts(rep(NA, n), order.by=tmpDates)

for (i in N:n) {
    thetaTemp <- theta[(i-N+1):i, ]
    thetaSummed[i] <- thetaSum(thetaTemp, w=rep(1, N))
}

thetaSummed

N是小于n的回顾期。

for循环有哪些快速替代方案？

Answer 1

您可以使用rollapply。

rollapplyr(theta, width=3, FUN=thetaSum, fill=NA)