cheetahfly 发表于 2017-8-31 16:23 
谢谢,亲测之后,举了一个特殊的例子,觉得不行呢。比方说,
a <- matrix(c(1:16),4,4)
A <- which(upper.tri(a),arr.ind = TRUE)
b <- c(10,5,9, 13, 14,15)
G1 <- b %% ncol(a)
G2 <- floor(b / ncol(a)) + 1
B <- interactionOrder <- cbind(G1,G2)
colnames(A) <- NULL
colnames(B) <- NULL
apply(apply(B, 1, function(y) {apply(A, 1, function(x) identical(x, y))}), 2, which)
会报错。