Question

我有一个数据框，每个人都有唯一的个人ID和测量值，例如

> personDat <- data.frame(personID = c("A","B","C","D"),value = rnorm(4))
> personDat
  personID      value
1        A -0.9246883
2        B  0.5175514
3        C -1.0109688
4        D  1.1614124

现在我需要创建所有独特的个体对，我可以使用combn函数：

> perCombs <- t(combn(personDat$personID,2))
> perCombs
     [,1] [,2]
[1,] A    B   
[2,] A    C   
[3,] A    D   
[4,] B    C   
[5,] B    D   
[6,] C    D   
Levels: A B C D

现在我想向perCombs添加两个额外的列，一个用于第一个personID的personDat中第二列的相应度量/值，然后是第二个personID的相应值。

基本上我需要对personID的唯一对的所有值进行二元运算，如果我有列，我可以用快速矢量化的方式进行。

编辑：这种天真的做法是：

perCombs <- data.frame(per1 = perCombs[,1],
                       per2 = perCombs[,2],
                       val1 = matrix(0,6,1),
                       val2 = matrix(0,6,1))
for(i in 1:6){
  print(i)
  perCombs[i,3] <- personDat[as.character(personDat$personID)==as.character(perCombs[i,1]),2]
  perCombs[i,4] <- personDat[as.character(personDat$personID)==as.character(perCombs[i,2]),2]
}

Answer 1

我们可以使用match

cbind(perCombs,  `dim<-`(personDat$value[match(perCombs, 
                         personDat$personID)], dim(perCombs)))

生成所有唯一的因子对并以相同的方式对齐其他变量R.

1 个答案: