计算列表中的出现次数，对吗？

debugcn 发表于 Dev

Richpiana

我有一个1和0的矩阵。有关此表的规则如下。

我想计算1，1系列的出现次数（其中1不被0隔开！），并为1，1，1系列做同样的事情。我已经尝试过colSums，但似乎不太合适。

矩阵最终为

      t1 t2 t3 t4 t5 t6 t7
 [1,]  0  0  0  0  1  1  0
 [2,]  0  0  1  1  0  0  1
 [3,]  1  1  0  0  0  0  0
 [4,]  0  0  1  1  1  0  0
 [5,]  0  0  1  1  0  0  0
 [6,]  1  1  0  0  0  0  0
 [7,]  0  0  0  0  0  0  1
 [8,]  0  0  0  0  1  1  0
 [9,]  1  1  0  0  1  1  0
[10,]  0  0  0  0  0  1  1
[11,]  1  1  0  0  0  0  0
[12,]  0  0  1  1  0  0  0
[13,]  0  0  0  0  0  0  0
[14,]  0  0  0  0  0  0  1
[15,]  0  0  0  0  0  0  0

因此，对于第一行，我想要1的序列为1,1，0为1的序列，1,1。对于第4行，我希望0表示1，1，但1等于1，1，1。

谁能告诉我下面的1,1系列的代码出了什么问题？

occ <- matrix()
occ_temp <- matrix
for (j in 1:nrow(final)){
  for (i in 2:7){
    if (sum(final[j,i-1:i])==2){occ_temp[j,i-1]=1}
  }
occ[j] <- sum(occ_temp)
}

阿克伦

我们可以使用遍历所有行apply，使用来获取run-length-type rle，提取lengths出valuesare为1的位置，检查是否等于'n1'和'n2'，然后获取sum。

n1 <- 2
n2 <- 3
res <- t(apply(m1, 1, FUN=function(x) {
       x1 <- with(rle(x), lengths[!!values])
       c(sum(x1==n1), sum(x1==n2))
}))
colnames(res) <- paste0("count", c(11, 111))
res
#      count11 count111
# [1,]       1        0
# [2,]       1        0
# [3,]       1        0
# [4,]       0        1
# [5,]       1        0
# [6,]       1        0
# [7,]       0        0
# [8,]       1        0
# [9,]       2        0
#[10,]       1        0
#[11,]       1        0
#[12,]       1        0
#[13,]       0        0
#[14,]       0        0
#[15,]       0        0

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。