r - 使用矩阵值作为索引

Question

我正在进行 Fisher 排列检验，其中我必须生成治疗状态的所有组合。

我们有 4 名受试者，其中 2 名接受了治疗。使用 combn，我可以生成治疗对象的所有组合。例如，第一行表示治疗了第一和第二受试者。

t(combn(4, 2))

     [,1] [,2]
[1,]    1    2
[2,]    1    3
[3,]    1    4
[4,]    2    3
[5,]    2    4
[6,]    3    4

如何从这个矩阵转到治疗状态矩阵，如下所示:

      [,1] [,2] [,3] [,4]
[1,]    1    1    0   0
[2,]    1    0    1   0
...

Answer 1

使用 base-R:

res <- t(apply(t(combn(4,2)),MARGIN=1,FUN=function(x){
  return(as.numeric(1:4 %in% x))
}))
> res
     [,1] [,2] [,3] [,4]
[1,]    1    1    0    0
[2,]    1    0    1    0
[3,]    1    0    0    1
[4,]    0    1    1    0
[5,]    0    1    0    1
[6,]    0    0    1    1