Question

我有一个TRUE和FALSE的逻辑向量。现在我想找出10个最长序列的TRUE或FALSE和（最重要的）它们的位置。与

set.seed(123)
data <- sample(seq(from = 1, to = 100, by = 1), size = 3000, replace = TRUE)
reference=3
logical.vec <- data > reference
logical.vec.rle <- rle(logical.vec)
MaxiSequence <- max(logical.vec.rle$length[logical.vec.rle$values=="TRUE"])

它会返回类似MaxiSequence = 55L的内容。我的问题是那个

它只获得最长的序列
位置丢失（我想知道数据在哪里或者是logical.vec）
它似乎并不正确，因为在logical.vec.rle中，$ length至少是60（见于＆＃34;环境＆＃34;在RStudio中）。 BTW typeof（MaxiSequence）返回整数。

Answer 1

使用data.table及其rleid - 功能，您可以执行以下操作：

DT <- data.table(logical.vec, runid = rleid(logical.vec))
DT_stat <- DT[,.(length = .N, position = .I[1], type=logical.vec[1]), by = runid]
# order by length
setorder(DT_stat, -length)
# Get top 3 (or n) by type
DT_stat[, .SD[1:3], type]

所以这导致：

    type runid length position
1:  TRUE     7    125      144
2:  TRUE   103    121     1623
3:  TRUE   115     96     1847
4: FALSE    34      2      717
5: FALSE    80      2     1351
6: FALSE    88      2     1491

R找到最大序列长度和位置

1 个答案: