R中的重叠矩阵

时间:2012-07-31 18:50:25

标签: r

我有以下数据框

id channel
1  a
1  b
1  c
2  a
2  c
3  a

我想创建并重叠矩阵。它基本上是一个带行和列的方阵。列标签为a,b,c,该表中的每个条目显示每个通道共有多少个id。例如,在上面的例子中,矩阵看起来像

  a b c
a 3 1 2
b 1 1 1
c 2 1 2

非常感谢。

2 个答案:

答案 0 :(得分:8)

这应该可以解决问题:

df <- data.frame(id=c(1,1,1,2,2,3), channel=letters[c(1,2,3,1,3,1)]) # your data

m <- table(df[[1]], df[[2]])   ## Alternatively:  m <- do.call(table, df)
t(m) %*% m
#   a b c
# a 3 1 2
# b 1 1 1
# c 2 1 2

答案 1 :(得分:5)

library(plyr)
df
  id channel
1  1       a
2  1       b
3  1       c
4  2       a
5  2       c
6  3       a
tb <- table(ddply(df, .(id), function(x) {x$id <- x$channel; expand.grid(x)}))
tb
   channel
id  a b c
  a 3 1 2
  b 1 1 1
  c 2 1 2
names(dimnames(tb)) <- NULL
tb
  a b c
a 3 1 2
b 1 1 1
c 2 1 2

现在有一些关于矩阵表的解释和一些内容作为table()的输出。 ?table

中有一个例子
a <- letters[1:3]
(b <- sample(a))
[1] "b" "c" "a"
table(a, b)
   b
a   a b c
  a 0 1 0
  b 0 0 1
  c 1 0 0

所以它按位置匹配元素..现在如果我们有

  id channel
  1       a
  1       b
  1       c
  2       a
  ...

然后,可以通过id拆分数据框来显示共享相同的id,创建channel列的副本并获取这两列的所有组合:

tbl <- expand.grid(data.frame(x = c("a","b","c"), y = c("a", "b", "c")))
tbl
  x y
1 a a
2 b a
3 c a
4 a b
5 b b
6 c b
7 a c
8 b c
9 c c
table(tbl$x, tbl$y)

    a b c
  a 1 1 1
  b 1 1 1
  c 1 1 1