r - 对数据框中具有相同值的列元素进行计数并以矩阵格式显示

Question

我在互联网上搜索类似的解决方案，但我找不到适合我的案例的具体解决方案。假设有以下数据框：

a = c(1, 1, 1, 2, 2)
b = c(2, 1, 1, 1, 2)
c = c(2, 2, 1, 1, 1)
d = c(1, 2, 2, 1, 1)
df <- data.frame(a = a, b = b, c = c, d = d)

df看起来像这样：

注意：在此示例中，我使用[1,2]了一对值，但它可以是一组不同的值：[-1,1]甚至两个以上的可能值：[-1,1,2].

现在我想要一个矩阵，其中每个[i,j]元素将代表行数，其值为1columni和j。对于这种特殊情况，我们有（显示上对角线，因为它是对称的）：

对角线应计算1给定列中具有值的行数。在这种情况下，所有列都具有相同数量的 value 1。格式应该类似于cor()函数（相关矩阵）。

我试图使用table()（也crosstab来自descr包），但它通过成对的列显示信息。

可以通过手动计算1每对列的出现（即：）nrow(df[df$a==1 & df$b==1,])=2然后放入一个矩阵来完成，但我想知道是否有一个内置函数可以简化这个过程。

score 3 · Accepted Answer

我们可以使用crossprodamatrix来计算1问题示例的值的出现：

m1 <- as.matrix(df == 1) # see Note[1]
out <- crossprod(m1)

注意[1]由@imo 指出（请参阅下面的评论），用于解决一般情况（具有值的矩阵：[x,y]）。对于具有[0,1]值的矩阵df==1可以替换为df. 要计算2问题示例中的值，请使用：df == 2。

如果下对角线应该是0或NA

out[lower.tri(out)] <- NA
out
#   a  b  c d
#a  3  2  1 1
#b NA  3  2 1
#c NA NA  3 2
#d NA NA NA 3

1 回答 1