数据框：用匹配单元格的相应值填充矩阵

Question

如果数据帧列中的单元格值匹配，我试图用"1"s填充相邻矩阵。“1”的位置基于同一行中的相应值。

更准确地说：数据帧pat1

ID PATID SUB
 1     2  2A
 2     2  2B 
 3     3  2C
 4     3  2D

我想填充矩阵单元格[2A，2B]，[2B，2A]，[2C，2D]，和[2D，2C]在空矩阵patmat1与相应的行/cole. name与"1"起PATID[1]=PATID[2]，和PATID[3]=PATID[4]，分别。

所需的输出将是带有列名的矩阵（数据=c（0,1,0,0,1,0,0,0,0,0,0,1,0,0,0,1,0），nrow=4，byrow=T）

我找过了，但还没有找到方法。

Answer 1

这有帮助吗

library(reshape2)
 if(length(unique(pat1$PATID))==1) pat1$SUB2 <- rev(pat1$SUB)
acast(pat1, SUB~SUB2, value.var='PATID', length)
#   2A 2B
#2A  0  1
#2B  1  0

对于新的数据集，我们可以按“PATID”进行split，对单个list元素执行acast，并使用bdiag折叠它们。

library(Matrix)
patmat1[] <- as.matrix(bdiag(lapply(lst, function(x) 
        acast(transform(x, SUB2=rev(SUB)), SUB~SUB2, 
        value.var='PATID', length))))
patmat1
#   2A 2B 2C 2D
#2A  0  1  0  0
#2B  1  0  0  0
#2C  0  0  0  1
#2D  0  0  1  0

或者我们可以使用数据来实现这一点。表。我们转换数据。帧“到”数据。表'（setDT（pat1）），按“PATID”分组，创建“SUB2”作为“SUB”的rev，然后使用acastfromrestrape2将“long”格式更改为“wide”格式，并指定fun。聚合为长度。

library(data.table)
acast(setDT(pat1)[, SUB2:= rev(SUB), PATID], SUB~SUB2, 
                    value.var='PATID', length)
#   2A 2B 2C 2D
#2A  0  1  0  0
#2B  1  0  0  0
#2C  0  0  0  1
#2D  0  0  1  0