R соответствует / сравнивает две матрицы в R, но цикл, кажется, не работает хорошо

Question

N = 4
comN <- gtools::permutations(2, N, v = c(0, 1), repeats = TRUE)
comstate <- matrix (NA, N, N)
for (i in 1:2^N) {
  for (j in 1:N) {
    comstate[1:N, ] = comN[i, ] # spoint<-comN[i,]
    tcomstate = t(comstate)     # need transpose "comstate" to ensure col stands for the component state
    tcomstate[j, j] = (1 - comN[i, j])
    y <- which(colSums(t(comN) == tcomstate[j, ]) == ncol(comN))
  }
}
y

Я пытался получить y вручную, (i = 3) comN [3,], это работает хорошо, y = 1,4,7,15, но когда я возвращаюсь к циклу, это только 15должно быть четыре значения 1,4,7,15