2017-05-28 60 views
1

在R中,我试图将指示符变量的数据帧转换为计数表。将指标列转换为计数表?

我的数据看起来像这样,其中A,B,C和D是不同条件下的指示变量,ID1,ID2,ID3和ID4是患者ID。

test <- data.frame(rbind(c(0,0,0,1),c(1,0,0,1),c(1,0,0,0),c(1,1,0,0))) 
colnames(test) <- c("A","B","C","D") 
rownames(test) <- c("ID1","ID2","ID3","ID4"); test 

    A B C D 
ID1 0 0 0 1 
ID2 1 0 0 1 
ID3 1 0 0 0 
ID4 1 1 0 0 

我想用它做的计数A,B,C,和d的摘要以及它们如何重叠的表来结束:

A B C D 
A 3 1 0 1 
B 1 1 0 0 
C 0 0 0 0 
D 1 0 0 2 

由于A,B ,C和D是不同的变量,一个给定的ID可以有多个条目,我不知道如何做到这一点。谢谢你的帮助!

回答

5

矩阵乘法是否可以解决您的问题?

test <- data.frame(rbind(c(0,0,0,1),c(1,0,0,1),c(1,0,0,0),c(1,1,0,0))) 
colnames(test) <- c("A","B","C","D") 
rownames(test) <- c("ID1","ID2","ID3","ID4"); test 

test1 <- as.matrix(test) 
t(test1) %*% test1 
+0

哦,我的天哪,这是完美的!非常感谢你 - 我永远不会想出这个! – statgirl259