rollapply与函数rle（x）

Question

我有时间序列数据作为data.table类，每列（观察点）都有我想在滑动窗口（30宽度）内计算它们的值。 我试图使用rle（sort（x））来计算rollapply中的每个值，但它不起作用。

例如，如果我有如下表，

dt <- data.frame(v1=c(1,0,1,4,4,4,4,4),v2=c(1,1,1,4,3,3,3,3),
          v3=c(0,1,1,3,3,3,3,2),v4=c(1,1,0,3,3,3,3,3),
       v5=c(1,1,1,5,5,5,5,5))

我试过这个;

rollapply(dt, 3, function(x) {rle(sort(x))$values; rle(sort(x))$length})

但结果却没有意义。 请给我一些方向......

Answer 1

解决方案1假设目标是获得3个值的滚动计数，请尝试以下操作：

m <- as.matrix(dt)
levs <- sort(unique(c(m)))
f <- function(x) table(factor(x, levs))
r <- rollapply(m, 3, f)

这里levs是0,1，...，5因此对于函数的每个应用程序，我们将得到一个向量6长度为0，1，...，5的计数。 有5个输入列，因此对每列应用这样的函数给出5 * 6 = 30列输出。

请注意， rollapply适用于矩阵或动物园对象，而不是数据框，因此我们对其进行了转换。 另外，为了确保每个函数应用程序输出相同长度的向量，我们将每个输入转换为具有相同级别的因子。

注意：

ra <- array(r, c(6, 6, 5))

给出一个3d数组，其中ra [i]是由rollapply(dt[, i], 3, f)形成的矩阵。 也就是说，在矩阵ra[,,i] ，对于列i上的f每个应用，存在一行，并且该行中的列计数0，1，...，5的数量。

另一种可能性是，它将相同的5个矩阵（每个输入列一个）作为结果列表的组件：

lapply(dt, rollapply, 3, f)

例如，请考虑以下内容。 输出的第1行表示ft在dt[,1]上的第一次应用有一个0，两个1，没有其他值。 这也可以从r[,,1]或lapply(dt, rollapply, 3, f)[[1]] ：

> rollapply(dt[, 1], 3, f)
     0 1 2 3 4 5
[1,] 1 2 0 0 0 0  <- dt[1:3,1] has 1 zero and 2 ones
[2,] 1 1 0 0 1 0  <- dt[2:4,1] has 1 zero and 1 one and 1 four, etc.
[3,] 0 1 0 0 2 0
[4,] 0 0 0 0 3 0
[5,] 0 0 0 0 3 0
[6,] 0 0 0 0 3 0

解决方案2

这说看输出的单元格1,1在dt[1:3,1]有一个0和两个1。 查看输出的单元格2,1我们看到dt[2:4,1]中有一个0,1个1和1个等等。

> g <- function(x) { tab <- table(x); toString(paste(names(tab), tab, sep = ":")) }
> sapply(dt, rollapply, 3, g) # or rollapply(m, 3, g) where m was defined in solution 1
     v1              v2              v3         v4              v5        
[1,] "0:1, 1:2"      "1:3"           "0:1, 1:2" "0:1, 1:2"      "1:3"     
[2,] "0:1, 1:1, 4:1" "1:2, 4:1"      "1:2, 3:1" "0:1, 1:1, 3:1" "1:2, 5:1"
[3,] "1:1, 4:2"      "1:1, 3:1, 4:1" "1:1, 3:2" "0:1, 3:2"      "1:1, 5:2"
[4,] "4:3"           "3:2, 4:1"      "3:3"      "3:3"           "5:3"     
[5,] "4:3"           "3:3"           "3:3"      "3:3"           "5:3"     
[6,] "4:3"           "3:3"           "2:1, 3:2" "3:3"           "5:3"     

补充：补充讨论和解决方案2。