Question

我有以下数据框。

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在“value”字段中找到前面的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来实现这一目标？非常感谢。

Answer 1

您可以使用基本功能轻松完成此操作：

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))

Answer 2

使用data.table v1.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))]

Answer 3

只想添加到@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）时，此附加项可防止出现这种情况

获取数据框中的前一个n行？

3 个答案: