给定一个向量,1:4,和序列长度,2,我想把向量分成“子向量”,每个子向量的长度为2,并生成这些子向量的所有可能组合的矩阵。
输出将如下所示:
[,1] [,2] [,3] [,4]
[1,] 1 2 3 4
[2,] 3 4 1 2另一个例子。如果向量1:8和子向量长度为4,输出如下所示:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,] 1 2 3 4 5 6 7 8
[2,] 5 6 7 8 1 2 3 4如果向量1:9和子向量长度为3,输出如下所示:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
[1,] 1 2 3 4 5 6 7 8 9
[2,] 1 2 3 7 8 9 4 5 6
[3,] 4 5 6 1 2 3 7 8 9
[4,] 4 5 6 7 8 9 1 2 3
[5,] 7 8 9 4 5 6 1 2 3
[6,] 7 8 9 1 2 3 4 5 6给定的向量长度必须可被子向量长度整除。
发布于 2020-04-28 03:45:36
我可以回答整个问题,但要花点时间。这应该会给你答案的味道。
包combinat有一个名为permn的函数,它为您提供了向量的所有排列。你想要这个,但不太想。你需要的是所有区块的排列。因此,在第一个例子中,有两个长度为2的块,在第二个示例中,有三个长度为3的块。如果我们看第一个,并考虑排序块
> library(combinat)
> numBlocks = 2
> permn(1:numBlocks)
[[1]]
[1] 1 2
[[2]]
[1] 2 1所以我希望你能看到,第一个排列将采用块b1 = c(1,2)和b2 = c(3,4),并给它们排序c(b1,b2),第二个排列将对它们进行c(b2,b1)排序。
同样,如果你有三个街区,那么b1 = 1:3; b2 = 4:6; b3 = 7:9
permn(1:3)
[[1]]
[1] 1 2 3
[[2]]
[1] 1 3 2
[[3]]
[1] 3 1 2
[[4]]
[1] 3 2 1
[[5]]
[1] 2 3 1
[[6]]
[1] 2 1 3给你这些块的顺序。更普遍的解决方案是找出如何移动块,但这并不难。
更新:使用我的multicool包的。Note共同词法排序(coolex)不是你自己想出来的顺序。
library(multicool)
combs = function(v, blockLength){
if(length(v) %% blockLength != 0){
stop("vector length must be divisible by blockLength")
}
numBlocks = length(v) / blockLength
blockWise = matrix(v, nc = blockLength, byrow = TRUE)
m = initMC(1:numBlocks)
Perms = allPerm(m)
t(apply(Perms, 1, function(p)as.vector(t(blockWise[p,]))))
}> combs(1:4, 2)
[,1] [,2] [,3] [,4]
[1,] 3 4 1 2
[2,] 1 2 3 4> combs(1:9, 3)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
[1,] 7 8 9 4 5 6 1 2 3
[2,] 1 2 3 7 8 9 4 5 6
[3,] 7 8 9 1 2 3 4 5 6
[4,] 4 5 6 7 8 9 1 2 3
[5,] 1 2 3 4 5 6 7 8 9
[6,] 4 5 6 1 2 3 7 8 9https://stackoverflow.com/questions/61471725
复制相似问题