按行过滤最长重复数序列的数据帧
按行过滤最长重复数序列的数据帧
提问于 2018-06-09 03:14:18
我正在尝试创建一个"filter-by“矩阵，我可以用它来隔离我的数据帧中的数据行，这样每一行只包含对应于相同数字的最长连续序列的值，而其余的则保持为零。在四处搜索后，我认为rle是要使用的函数，但这并没有给我想要的东西。以下是我的代码和结果的示例。我们将非常感谢您的建议和解决方案。谢谢!
示例数据：
a<- c(1,0,1,1,1,1,0,0)
b<- c(0,0,0,1,1,1,0,1)
c<- c(0,0,1,1,0,0,0,1)
d<- c(1,0,0,1,1,1,1,0)
e<- c(1,0,0,1,0,0,1,1)
f<- c(0,0,0,1,1,1,0,1)
g<- c(0,0,1,1,0,0,0,1)
test.data <- data.frame(cbind(a,b,c,d,e,f,g))
# > test.data
# a b c d e f g
# 1 1 0 0 1 1 0 0
# 2 0 0 0 0 0 0 0
# 3 1 0 1 0 0 0 1
# 4 1 1 1 1 1 1 1
# 5 1 1 0 1 0 1 0
# 6 1 1 0 1 0 1 0
# 7 0 0 0 1 1 0 0
# 8 0 1 1 0 1 1 1尝试解决方案的示例代码：
result <- data.frame(lapply(test.data, function(x) {
r <- rle(x)
r$values[r$lengths!=max(r$lengths)]==1
r2=inverse.rle(r)
r2
}))我得到的结果(看起来和放进去的东西完全一样？)：
# > result
# a b c d e f g
# 1 1 0 0 1 1 0 0
# 2 0 0 0 0 0 0 0
# 3 1 0 1 0 0 0 1
# 4 1 1 1 1 1 1 1
# 5 1 1 0 1 0 1 0
# 6 1 1 0 1 0 1 0
# 7 0 0 0 1 1 0 0
# 8 0 1 1 0 1 1 1这是我想要得到的结果(T/F可以用来代替1和0，如果容易的话)：
# > result
# a b c d e f g
# 1 0 0 0 1 1 0 0
# 2 0 0 0 0 0 0 0
# 3 0 0 0 0 0 0 0
# 4 1 1 1 1 1 1 1
# 5 1 1 0 0 0 0 0
# 6 1 1 0 0 0 0 0
# 7 0 0 0 1 1 0 0
# 8 0 0 0 0 1 1 1敬请指教！
回答 2
Stack Overflow用户
回答已采纳
发布于 2018-06-09 03:35:09
library(magrittr)
val <- 1
test.data %>%
apply(1, function(x){
rle(x) %$% {
if(all(values != val)) rep(0, length(x))
else {
m <- max(lengths[values == val])
# Get only longest sequences
values <- (lengths == m & values == val)*values*(m > 1)
# Get only one of them
values[seq_along(values) != which(values == val)[1]] <- 0
rep(values, lengths)
}
}}) %>% t
# [,1] [,2] [,3] [,4] [,5] [,6] [,7]
# [1,] 0 0 0 1 1 0 0
# [2,] 0 0 0 0 0 0 0
# [3,] 0 0 0 0 0 0 0
# [4,] 1 1 1 1 1 1 1
# [5,] 1 1 0 0 0 0 0
# [6,] 1 1 0 0 0 0 0
# [7,] 0 0 0 1 1 0 0
# [8,] 0 0 0 0 1 1 1Stack Overflow用户
发布于 2018-06-09 03:28:32
我想这就是你想要的..。
test.data[] <- t(apply(test.data,1,function(x) {y<-rle(x)
y$values[y$lengths==1] <- 0
y$values[y$lengths!=max(y$lengths[y$values==1])] <- 0
return(inverse.rle(y))}))
test.data
a b c d e f g
1 0 0 0 1 1 0 0
2 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0
4 1 1 1 1 1 1 1
5 1 1 0 0 0 0 0
6 1 1 0 0 0 0 0
7 0 0 0 1 1 0 0
8 0 0 0 0 1 1 1页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/50767114
复制相关文章
相似问题