R data.table 查找当前行与上一行之间的滞后
R data.table find lags between current row to previous row
> tempDT <- data.table(colA = c("E","E","A","A","E","A","E")
+ , lags = c(NA,1,1,2,3,1,2))
> tempDT
colA lags
1: E NA
2: E 1
3: A 1
4: A 2
5: E 3
6: A 1
7: E 2
我有第 colA
列,需要找出当前行和前一行 colA == "E"
.
之间的滞后
注意:如果我们能找到 colA == "E"
的前一行的行引用,那么我们就可以计算滞后。但是,我不知道如何实现。
1)定义lastEpos
给定i
returns最后E
在第一个[=中的位置16=] 行并将其应用于每个行号:
lastEpos <- function(i) tail(which(tempDT$colA[1:i] == "E"), 1)
tempDT[, lags := .I - shift(sapply(.I, lastEpos))]
这里有一些变化:
2) i-1 在这个变化中 lastEpos
returns 最后 E
在第一个 i-1
中的位置] 行而不是 i
:
lastEpos <- function(i) tail(c(NA, which(tempDT$colA[seq_len(i-1)] == "E")), 1)
tempDT[, lags := .I - sapply(.I, lastEpos)]
3) Position 类似于 (2) 但使用 Position
:
lastEpos <- function(i) Position(c, tempDT$colA[seq_len(i-1)] == "E", right = TRUE)
tempDT[, lags := .I - sapply(.I, lastEpos)]
4) 滚动应用
library(zoo)
w <- lapply(1:nrow(tempDT), function(i) -rev(seq_len(i-1)))
tempDT[, lags := .I - rollapply(colA == "E", w, Position, f = c, right = TRUE)]
5) sqldf
library(sqldf)
sqldf("select a.colA, a.rowid - b.rowid lags
from tempDT a left join tempDT b
on b.rowid < a.rowid and b.colA = 'E'
group by a.rowid")
> tempDT <- data.table(colA = c("E","E","A","A","E","A","E")
+ , lags = c(NA,1,1,2,3,1,2))
> tempDT
colA lags
1: E NA
2: E 1
3: A 1
4: A 2
5: E 3
6: A 1
7: E 2
我有第 colA
列,需要找出当前行和前一行 colA == "E"
.
注意:如果我们能找到 colA == "E"
的前一行的行引用,那么我们就可以计算滞后。但是,我不知道如何实现。
1)定义lastEpos
给定i
returns最后E
在第一个[=中的位置16=] 行并将其应用于每个行号:
lastEpos <- function(i) tail(which(tempDT$colA[1:i] == "E"), 1)
tempDT[, lags := .I - shift(sapply(.I, lastEpos))]
这里有一些变化:
2) i-1 在这个变化中 lastEpos
returns 最后 E
在第一个 i-1
中的位置] 行而不是 i
:
lastEpos <- function(i) tail(c(NA, which(tempDT$colA[seq_len(i-1)] == "E")), 1)
tempDT[, lags := .I - sapply(.I, lastEpos)]
3) Position 类似于 (2) 但使用 Position
:
lastEpos <- function(i) Position(c, tempDT$colA[seq_len(i-1)] == "E", right = TRUE)
tempDT[, lags := .I - sapply(.I, lastEpos)]
4) 滚动应用
library(zoo)
w <- lapply(1:nrow(tempDT), function(i) -rev(seq_len(i-1)))
tempDT[, lags := .I - rollapply(colA == "E", w, Position, f = c, right = TRUE)]
5) sqldf
library(sqldf)
sqldf("select a.colA, a.rowid - b.rowid lags
from tempDT a left join tempDT b
on b.rowid < a.rowid and b.colA = 'E'
group by a.rowid")