powered by
collapseMat(m1, m2, collapse_matrix = matrix(c(0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 3L, 0L, 1L, 3L, 3L), 4, 4))