交叉产品的对面：如何从两个矩阵的交集创建一个新的矩阵？

Question

我有两个R表（女性和男性）与存在缺勤数据。我想在它们之间进行两两比较（逐行）以找出每对之间不共享的细胞数量（即细胞总数等于女性中的1，而不是男性中，反之亦然） 。

我知道交叉产品（％*％）与我所需要的相反。它创建了一个新的矩阵，其中包含男性和女性之间共享单元格的总和（即两个单元格的总和等于1）。

下面是一个例子数据集：

females <- as.data.frame(matrix(c(0,0,0,1,1,0,1,0,1,0,1,0,1,0,1,1,1,0,1,1,1,0,1,1,1), nrow=5, byrow=T)) 
males <- as.data.frame(matrix(c(1,0,0,1,1,0,1,0,1,1,1,0,1,0,1,1,1,0,1,1,1,0,1,0,1), nrow=5, byrow=T)) 
rownames(females) <-c ("female_1","female_2","female_3","female_4","female_5") 
rownames(males) <-c ("male_1","male_2","male_3","male_4","male_5") 

所以，如果我做跨产品

as.matrix(females) %*% t(as.matrix(males)) 

我得到这个

  male_1 male_2 male_3 male_4 male_5 
female_1  2  2  1  2  1 
female_2  1  2  0  2  0 
female_3  2  1  3  2  3 
female_4  3  3  2  4  2 
female_5  3  2  3  3  3 

但我需要这个（仅第一行显示）

  male_1 male_2 male_3 male_4 male_5 
female_1  1  1  3  2  3 
. 
. 

实际上，我的数据集不是对称的（我有47位女性和32位男性）。

感谢您的帮助！

Answer 1

建立一个对象，以接收结果：呼吁

xy <- matrix(NA, nrow(females), nrow(males)) 
for (x in 1:nrow(females)){ 
     for(y in 1:nrow(males)){ 
       xy[x,y] <- sum(females[x, 1:ncol(females)] != males[y,1:ncol(males)])}} 

应与嵌套sapply已经做的很好，可能是有点清洁，因为没有需要有一个单独的“设置”，（但只有一点点清洁和流行的神话相反没有任何更快）：

xy <- sapply(1:nrow(females) , 
       function(x) sapply(1:nrow(males) , 
        function(y) sum(females[x, 1:ncol(females)] != males[y,1:ncol(males)]))) 
xy 
#----- 
    [,1] [,2] [,3] [,4] [,5] 
[1,] 1 3 2 1 1 
[2,] 1 1 4 1 3 
[3,] 3 5 0 3 1 
[4,] 2 2 3 0 2 
[5,] 3 5 0 3 1 

dimnames(xy) <- list(rownames(females), rownames(males)) 

Answer 2

inverseCross <- function(females, males){ 
inverse <- data.frame(a=integer(), b=integer(), c=integer(), d=integer(), e=integer()) 
tempRow <- NULL 
for(i in 1:nrow(females)){ 
for(j in 1:nrow(males)){ 
sum <- 0 
for(k in 1: ncol(males)){ 
if(females[i,k] != males[j,k]){ 
sum <- sum + 1 
} 
} 
tempRow <- c(tempRow, sum) 
} 
inverse[i,] <- tempRow 
} 
colnames(inverse) <- rownames(males) 
rownames(inverse) <- rownames(females) 
inverse 
}