0

我想在 r 中构建一个固定的滚动窗口。假设我有一个这样的数据集：

      Apple Microsoft     Amazon      Tesla
2000 0.91903890 0.5454254 0.22267026 0.41857974
2001 0.13638076 0.7640585 0.56990815 0.04490846
2002 0.19977390 0.9170585 0.04391331 0.72987722
2003 0.70627675 0.2583974 0.03485766 0.35594681
2004 0.08069703 0.2085965 0.19865918 0.30194120
2005 0.03334810 0.7955248 0.75876036 0.28129544
2006 0.94652374 0.6095904 0.98855677 0.36792881
2007 0.90060660 0.3144541 0.78201742 0.02731390

我知道以下功能会给我一个扩展窗口：

 all.cov.matrix <- lapply(1:nrow(stocks), function(y) cov(stocks[1:y,]))

为了获得一个内部有四个句点的固定窗口，我尝试了以下功能：

library(zoo)     
all.cov.matrix <- apply(rollapply(1:nrow(stocks), 4, c), 1, function(ix) cov(stocks[ix, ]))
    # However, this returns one big matrix, that is not what I am looking for. 
    # Ideally I want to get the following results:
     cov(stocks[1:4,]) # 1 period
     cov(stocks[2:5,]) # 2 period and so on

我希望将每个矩阵单独存储在 中all.cov.matrix，因此在 thix 示例中all.cov.matrix应该存储 5 个不同的矩阵。

4

2 回答 2

3

您需要使用apply(..., margin = 2, ...)来获得柱状协方差。但是，我不推荐apply。您可以lapply改用，即

library(zoo)
lapply(as.data.frame(t(rollapply(1:nrow(stocks), 4, c))), function(i) cov(stocks[i,]))

另一种写法（根据@G.Grothendieck 评论），

r <- rollapplyr(stocks, 4, function(x) c(cov(x)), by.column = FALSE)
lapply(split(r, row(r)), matrix, 4)
于 2018-06-15T13:57:10.723 回答
3

如果您的扩展窗口示例接近您想要的，您可以通过一些细微的调整使其工作，以使“窗口”以您想要的方式滚动：

all.cov.matrix <- lapply( 2:(nrow(stocks)-2),
                          function(y) {
                              cov(stocks[(y-1):(y+2),])
                          } )

您的方式1:y意味着窗口始终从 开始1，而y每次迭代都会增加，这意味着窗口会扩展。将其更改为(y-1):(y+2)意味着您的窗口的起点和终点都会随之移动y，始终制作一个大小为的窗口4。

为了注意不要超出数据的边界，我们还需要更改为X向量 in lapply（第一个参数），以使窗口的边缘保持在数据内。

于 2018-06-15T14:01:01.943 回答