如果
如何使用rollapply
(或其他R函数)随着函数在数据中的进展而增加窗口大小?换句话说,第一个应用程序适用于第一个元素,第二个适用于前两个元素,第三个适用于前三个元素,等等。
如果您希望应用min
, max
, sum
或prod
,这些函数已经有它们的累积对应项:
cummin
, cummax
, cumsum
, cumprod
要在增长/扩展的窗口上应用更多奇特的函数,您可以简单地使用sapply
如
# your vector of interest
x <- c(1,2,3,4,5)
sapply(seq_along(x), function(y,n) yourfunction(y[seq_len(n)]), y = x)
对于基本的动物园对象
x.Date <- as.Date("2003-02-01") + c(1, 3, 7, 9, 14) - 1
x <- zoo(rnorm(5), x.Date)
# cumsum etc will work and return a zoo object
cs.zoo <- cumsum(x)
# convert back to zoo for the `sapply` solution
# here `sum`
foo.zoo <- zoo(sapply(seq_along(x), function(n,y) sum(y[seq_len(n)]), y= x), index(x))
identical(cs.zoo, foo.zoo)
## [1] TRUE
从窥视文档在?zooapply
我认为这将做你想要的,其中a
是你的矩阵和sum
可以是任何功能:
a <- cbind(1:5,1:5)
# [,1] [,2]
# [1,] 1 1
# [2,] 2 2
# [3,] 3 3
# [4,] 4 4
# [5,] 5 5
rollapply(a,width=seq_len(nrow(a)),sum,align="right")
# [,1] [,2]
# [1,] 1 1
# [2,] 3 3
# [3,] 6 6
# [4,] 10 10
# [5,] 15 15
除了@mnel的回答:
对于更奇特的函数,您可以简单地使用sapply
sapply
方法花费的时间太长,那么您最好迭代地表述您的函数。