y4 <- matrix(sample(0:1, 40, replace = T), nrow = 10)
pair <- as.matrix(expand.grid(col1 = 1:ncol(y4), col2 = 1:ncol(y4)))
pair <- pair[pair[, 1] < pair[, 2], ]
comp <- apply(pair, 1, function(x) y4[, x[1]] & y4[, x[2]])
colnames(comp) <- apply(pair, 1, function(x) paste0("m", x, collapse=""))
cbind(y4, comp)