通过块选择R矩阵值

Question

我有一个形式的矩阵：

[,1] [,2]
1    0
100  0
200  1
300  1
400  1
500  0
600  0
700  1
800  1
900  1

我正在尝试制作一个R脚本，该脚本在第2列上选择至少3个连续的“ 1”值，例如，检索

200  1
300  1
400  1
700  1
800  1
900  1

但找不到简单的方法。 有没有人遇到过这样的问题？ 我真的很感激！！！ 谢谢！！！

Answer 1

我将使用rle和inverse.rle如下：

r <- rle(mat[, 2])
r$values[r$values != 1 | r$lengths < 3] <- 0
keep <- as.logical(inverse.rle(r))
mat[keep, ]

Answer 2

您可以使用rle 。 假设您的矩阵是M ，感兴趣的列是2 ，那么您可以使用以下方法创建一个过滤器：

filter <- with(rle(M[,2]), rep(lengths>=3 & values==1, lengths))

然后使用过滤矩阵

M[filter,]