統計變量或變量組合的實例數量爲真

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