我有一个名为“df”的数据框,如下所示:
ID Value
1 a
1 b
1 c
1 d
3 a
3 b
3 e
3 f
. .
. .
. .
我有一个用这样的零填充的矩阵:
a b c d e f
a x 0 0 0 0 0
b 0 x 0 0 0 0
c 0 0 x 0 0 0
d 0 0 0 x 0 0
e 0 0 0 0 x 0
f 0 0 0 0 0 x
然后我想循环遍历这样的数据框:
for each ID, for each value i, for each value j != i, matrix[i,j] += 1
因此,对于每个ID,对于每个值的组合,我想将矩阵中的值提高1,结果是:
a b c d e f
a x 2 1 1 1 1
b 2 x 1 1 1 1
c 1 1 x 1 0 0
d 1 1 1 x 0 0
e 1 1 0 0 x 1
f 1 1 0 0 1 x
因此,例如,[a,b] = 2,因为这两个值的组合出现了两个不同的ID,而[a,c] = 1,因为这种值的组合仅在ID = 1时发生,而不在ID时发生= 3.
我怎样才能做到这一点?我已经制作了一个包含唯一ID的向量。
提前致谢。
答案 0 :(得分:5)
最简单的方法是获取table
,然后执行crossprod
out <- crossprod(table(df))
diag(out) <- NA #replace the diagonals with NA
names(dimnames(out)) <- NULL #set the names of the dimnames as NULL
out
# a b c d e f
#a NA 2 1 1 1 1
#b 2 NA 1 1 1 1
#c 1 1 NA 1 0 0
#d 1 1 1 NA 0 0
#e 1 1 0 0 NA 1
#f 1 1 0 0 1 NA
df <- structure(list(ID = c(1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L), Value = c("a",
"b", "c", "d", "a", "b", "e", "f")), .Names = c("ID", "Value"
), class = "data.frame", row.names = c(NA, -8L))