将矢量转换为逻辑矩阵

Question

是否有一个本地R函数，它将采用输入向量并返回相应的二进制矩阵，其中矩阵与输入向量中的唯一值具有相同的列数？

例如，给定x <- 1:3，我想回到下面的矩阵：

 [,1] [,2] [,3] 
[1,] 1 0 0 
[2,] 0 1 0 
[3,] 0 0 1 

功能contrasts靠拢，但我似乎无法避开第n-1列返回：

> contrasts(as.factor(x)) 
    2 3 
1 0 0 
2 1 0 
3 0 1 

Answer 1

model.matrix()可能有助于在这里，但你需要压制截距：

> model.matrix(~ factor(1:3) - 1) 
    factor(1:3)1 factor(1:3)2 factor(1:3)3 
1   1   0   0 
2   0   1   0 
3   0   0   1 
attr(,"assign") 
[1] 1 1 1 
attr(,"contrasts") 
attr(,"contrasts")$`factor(1:3)` 
[1] "contr.treatment" 

东西稍微复杂一些：

> set.seed(1) 
> fac <- factor(sample(1:3, 10, replace = TRUE)) 
> model.matrix(~ fac - 1) 
    fac1 fac2 fac3 
1  1 0 0 
2  0 1 0 
3  0 1 0 
4  0 0 1 
5  1 0 0 
6  0 0 1 
7  0 0 1 
8  0 1 0 
9  0 1 0 
10 1 0 0 
attr(,"assign") 
[1] 1 1 1 
attr(,"contrasts") 
attr(,"contrasts")$fac 
[1] "contr.treatment" 

Answer 2

其实，contrasts是你想要的。

contrasts(as.factor(1:3), contrasts=FALSE) 

    1 2 3 
1 1 0 0 
2 0 1 0 
3 0 0 1