2

假设我有一个具有以下结构的数据框:

> DF <- data.frame(x=1:5, y=6:10)
> DF
  x  y
1 1  6
2 2  7
3 3  8
4 4  9
5 5 10

我需要构建一个新的数据框,其中包含来自第一个数据框的重叠观察结果,以用作构建 Rglpk 优化库的 A 矩阵的输入。我将使用 n 长度的观察窗口,这样如果 n=2,则生成的数据框将连接第 1&2、2&3、3&4 等行。结果数据帧的长度为

(numberOfObservations-windowSize+1)*windowSize

这个例子的结果 windowSize=2 将是一个类似的结构

  x  y
1 1  6
2 2  7
3 2  7
4 3  8
5 3  8
6 4  9
7 4  9
8 5 10

我可以做一个循环

DFResult <- NULL
numBlocks <- nrow(DF)-windowSize+1
for (i in 1:numBlocks) {
    DFResult <- rbind(DFResult, DF[i:(i+horizon-1), ])
}

但这似乎效率很低,尤其是对于非常大的数据帧。

我也试过

rollapply(data=DF, width=windowSize, FUN=function(x) x, by.column=FALSE, by=1)
     x y
[1,] 1 6
[2,] 2 7
[3,] 2 7
[4,] 3 8

我试图在不应用任何聚合函数的情况下重复一行行。这不起作用,因为我缺少一些行

我对此感到有些困惑,并四处寻找类似的问题,但找不到任何问题。有没有人有更好的想法?

4

1 回答 1

2

我们可以做一个矢量化的方法

i1 <- seq_len(nrow(DF))
res <- DF[c(rbind(i1[-length(i1)], i1[-1])),]
row.names(res) <- NULL   
res
#  x  y
#1 1  6
#2 2  7
#3 2  7
#4 3  8
#5 3  8
#6 4  9
#7 4  9
#8 5 10
于 2017-09-04T11:08:24.567 回答