6

我有以下数据框。

date id value
2012-01-01 1 0.3
2012-01-01 2 0.5
2012-01-01 3 0.2
2012-01-01 4 0.8
2012-01-01 5 0.2
2012-01-01 6 0.8
2012-01-01 7 0.1
2012-01-01 8 0.4
2012-01-01 9 0.3
2012-01-01 10 0.2

有几个日期，对于每个日期，我有 10 个 id 值，如上所示和一个值字段。我想做的是为每个 id 找到“值”字段中的前 n 个值。例如，如果 n = 3，那么我希望输出如下。

date id value value1 value2 value3
2012-01-01 1 0.3 NA NA NA
2012-01-01 2 0.5 NA NA NA
2012-01-01 3 0.2 NA NA NA
2012-01-01 4 0.8 0.2 0.5 0.3
2012-01-01 5 0.2 0.8 0.2 0.5
...

有没有一种简单的方法可以通过 plyr 或使用 mapply 来实现？提前非常感谢。

4

3 回答 3

6

您可以使用基本函数很容易地做到这一点：

id <- 1:10
value <- c(0.3,0.5,0.2,0.8,0.2,0.8,0.1,0.4,0.3,0.2)
test <- data.frame(id,value)

test$valprev1 <- c(rep(NA,1),head(test$value,-1))
test$valprev2 <- c(rep(NA,2),head(test$value,-2))
test$valprev3 <- c(rep(NA,3),head(test$value,-3))

结果

   id value valprev1 valprev2 valprev3
1   1   0.3       NA       NA       NA
2   2   0.5      0.3       NA       NA
3   3   0.2      0.5      0.3       NA
4   4   0.8      0.2      0.5      0.3
5   5   0.2      0.8      0.2      0.5
6   6   0.8      0.2      0.8      0.2
7   7   0.1      0.8      0.2      0.8
8   8   0.4      0.1      0.8      0.2
9   9   0.3      0.4      0.1      0.8
10 10   0.2      0.3      0.4      0.1

之前在这里犯了一个错误 - 这是sapply一个函数中的一个版本：

prevrows <- function(data,n) {sapply(1:n,function(x) c(rep(NA,x),head(data,-x)))}
prevrows(test$value,3)

这给出了这个：

      [,1] [,2] [,3]
 [1,]   NA   NA   NA
 [2,]  0.3   NA   NA
 [3,]  0.5  0.3   NA
 [4,]  0.2  0.5  0.3
 [5,]  0.8  0.2  0.5
 [6,]  0.2  0.8  0.2
 [7,]  0.8  0.2  0.8
 [8,]  0.1  0.8  0.2
 [9,]  0.4  0.1  0.8
[10,]  0.3  0.4  0.1

然后，您可以将其应用于数据中的每组日期，如下所示：

result <- tapply(test$value,test$date,prevrows,3)

它为每个日期集提供了一堆列表。您可以使用以下方法将这些绑定起来以添加回您的数据集：

data.frame(test,do.call(rbind,result))
于 2012-05-29T06:25:54.780 回答
3

使用data.tablev1.9.5+ 这很简单：

library(data.table)
setDT(dt)

lags <- dt[, shift(value, n = c(1,2,3))]

或将它们作为附加列附加到同一 data.table 中：

dt[, c("lag1", "lag2", "lag3") := shift(value, n = c(1,2,3))]
于 2016-08-03T21:39:06.127 回答
0

只想添加到@thelatemail 的答案（我无法直接评论我的声誉）：

prevrows2 <- function(data,n) {
if (length(data) >= 10){
sapply(1:n,function(x) c(rep(NA,x),head(data,-x)))
} else {
cbind(sapply(1:length(data),function(x) c(rep(NA,x),head(data,-x))),
matrix(NA,nrow = length(data),ncol= n - length(data)))}
}

此添加可防止组中的行数小于您希望选择的行数 (n) 的情况

于 2019-12-18T18:38:13.360 回答