[英]Calculate / Translate vector of numbers in binary matrix / vector in R
我如何自动构建一个矩阵,将vector <- c(1,2,3)
的排列转换为某种二进制格式? 像这样:
x <- matrix(c(1,1,0,0,0,1,0,1,1,0,1,1,0,0,0,1,1,1), ncol = 3)
rownames(x) <- c("1", "1,2", "2", "3", "2,3", "1,2,3")
colnames(x) <- c("1", "2", "3")
x
1 2 3
1 1 0 0
1,2 1 1 0
2 0 1 0
3 0 0 1
2,3 0 1 1
1,2,3 1 1 1
虽然我想要的不是 3 个,而是 7 个值。
一种方法是
x <- c(2, 4, 5)
combs <- sapply(1:length(x), combn, x = x)
M <- do.call(rbind, sapply(combs, function(u)
t(apply(u, 2, function(v) 1 * x %in% v))))
dimnames(M) <- list(unlist(sapply(combs, apply, 2, paste, collapse = ",")), x)
M
# 2 4 5
# 2 1 0 0
# 4 0 1 0
# 5 0 0 1
# 2,4 1 1 0
# 2,5 1 0 1
# 4,5 0 1 1
# 2,4,5 1 1 1
这是一个将任何向量转换为适当的二进制矩阵的函数,
get_binary <- function(x){
v1 <- unlist(sapply(seq_along(x), function(i) combn(x, i, toString)))
mat <- t(sapply(v1, function(i)sapply(x, function(j) as.integer(grepl(j, i)))))
colnames(mat) <- x
return(mat)
}
get_binary(c(2, 8, 9))
这使,
2 8 9 2 1 0 0 8 0 1 0 9 0 0 1 2, 8 1 1 0 2, 9 1 0 1 8, 9 0 1 1 2, 8, 9 1 1 1
这是一个非常快的单行代码:
vector <- c(1,2,3)
library(RcppAlgos)
toBinary <- function(v) permuteGeneral(0:1, length(v), TRUE)[-1,]
toBinary(vector)
[,1] [,2] [,3]
[1,] 0 0 1
[2,] 0 1 0
[3,] 0 1 1
[4,] 1 0 0
[5,] 1 0 1
[6,] 1 1 0
[7,] 1 1 1
[-1, ]
是删除所有零的行。 这一行代表幂集中的空集。 实际上,从技术上讲,您要求的是从向量的幂集(当然减去空集)到二进制矩阵的映射。
如果您真的希望row.names
成为实际的排列,您可以使用rje
包中的powerSet
函数。 观察:
library(rje)
nameTest <- toBinary(vector)
row.names(nameTest) <- lapply(powerSet(rev(vector))[-1], sort)
nameTest
[,1] [,2] [,3]
3 0 0 1
2 0 1 0
c(2, 3) 0 1 1
1 1 0 0
c(1, 3) 1 0 1
c(1, 2) 1 1 0
c(1, 2, 3) 1 1 1
*免责声明:我是RcppAlgos
的作者
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.