我有一个矢量,说vec1,另一个矢量vec2如下:
vec1 = c(4,1) # [1] 4 1 vec2 = c(5,3,2) # [1] 5 3 2
我正在寻找的是vec1和vec2的所有可能的组合,同时保持向量元素的顺序.也就是说,得到的矩阵应该是这样的:@H_301_5@
> res [,1] [,2] [,3] [,4] [,5] [1,] 4 1 5 3 2 [2,] 4 5 1 3 2 [3,] 4 5 3 1 2 [4,] 4 5 3 2 1 [5,] 5 4 1 3 2 [6,] 5 4 3 1 2 [7,] 5 4 3 2 1 [8,] 5 3 4 1 2 [9,] 5 3 4 2 1 [10,] 5 3 2 4 1 # res=structure(c(4,4,5,1,# 3,2,# 2,1),.Dim = c(10L,5L))
两个向量不允许重复.也就是说,所得到的矩阵的所有行都具有唯一的元素.@H_301_5@
我实际上正在寻找最有效的方式.解决这个问题的一种方法是生成长度为n的所有可能的排列(n = 5),然后应用过滤.但是,随着n的增长,它是耗时的.@H_301_5@
有没有办法做到这一点?@H_301_5@
解决方法
尝试这个:
nv1 <- length(vec1) nv2 <- length(vec2) n <- nv1 + nv2 result <- combn(n,nv1,function(v) {z=integer(n);z[v]=vec1;z[-v]=vec2;z})
这个想法是产生所有组合的索引,以便放置vec1的元素.@H_301_5@