统计变量或变量组合的实例数量为真

Question

我是一位热情的R新手，需要一些帮助！ :)

我有看起来像这样的数据帧：

id<-c(100,200,300,400) 
a<-c(1,1,0,1) 
b<-c(1,0,1,0) 
c<-c(0,0,1,1) 

y=data.frame(id=id,a=a,b=b,c=c) 

凡ID是一个唯一标识符（例如人）和一个，b和Ç是虚设变量是否具有此功能（如始终为1 = TRUE）。

我希望R创建一个矩阵或数据框，其中变量a，b和c都作为列和行的名称。对于矩阵R的值将必须计算具有此特征的标识符的数量或特征的组合。

因此，例如，ID为100,200和400的矩阵的对角线上有一个特征a，其中a和十字形R将输入3.只有ID 100具有特征a和b，因此R将输入1 a和b交叉，等等。

产生的数据帧将不得不像这样：

l<-c("","a","b","c") 
m<-c("a",3,1,1) 
n<-c("b",1,2,1) 
o<-c("c",1,1,2) 
result<-matrix(c(l,m,n,o),nrow=4,ncol=4) 

由于我的数据集有10个变量和数以百计的观察，我会自动的全过程。

您的帮助将不胜感激。 非常感谢！

Answer 1

随着基础R：

crossprod(as.matrix(y[,-1])) 
# a b c 
# a 3 1 1 
# b 1 2 1 
# c 1 1 2 

Answer 2

这被称为邻接矩阵。你可以用qdap包做到这一点很容易地：

library(qdap) 
adjmat(y[,-1])$adjacency 

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

它，因为你喂养它一个数据帧抛出一个警告。没什么大不了的，可以忽略。也注意到我放弃了负号索引y[, -1]的第一列（ID）。

需要注意的是，因为你有一个布尔矩阵开始时你可能会得到与有：

Y <- as.matrix(y[,-1]) 
t(Y) %*% Y