将矩阵排列成唯一的行和列

1 回答

• 1

  首先，有一种更好的方法来获取矢量 . 这将减少打字并使您的代码更清晰：

  vec <- rep(1:8, each=4)
  mat <- matrix(vec, ncol=6, byrow=T)
  

  接下来，请注意您要按行排序元素的唯一排序 . 也就是说，你想要你的元素的permutations . N 元素的唯一排列数为 N! . 我们可以在 R 中计算出这个数字：

  factorial(length(unique(vec)))
  # [1] 40320
  

  这比您需要的数字排序更多 . （看起来你只想要6.）假设你想绘制一些随机排序，让我们从可能的40320中抽取6个数字：

  set.seed(3931)  # this makes the example exactly reproducible
  rows <- sample.int(40320, 6)
  rows
  # [1] 36601 24136  5713 23405 25328 32973
  

  从这里开始，为了得到实际的排列，我们将使用 R 中的permute库 . 因为排列的总数很大，我们必须将 maxperm 参数设置为大于其默认值（ 9999 ，我收集）：

  library(permute)
  out <- allPerms(unique(vec), control=how(maxperm=41000))[rows,]
  out
  #      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
  # [1,]    8    2    7    1    3    4    6    5
  # [2,]    5    7    4    1    6    8    2    3
  # [3,]    2    1    8    6    3    4    7    5
  # [4,]    5    6    4    1    2    8    7    3
  # [5,]    6    1    3    2    5    7    4    8
  # [6,]    7    4    6    5    8    2    3    1
  

  ---

  我为误读这个问题而道歉 . N 元素的排列数，当你只使用 r 时，它们是 N!/r! . 这是一个简单的计算：

  factorial(length(unique(vec))) / factorial(length(unique(vec))-6)
  # [1] 20160
  

  我不清楚如何让 allPerms() 使用少于所有元素 . 这是一个简单的解决方案：

  library(permute)
  out <- allPerms(unique(vec), control=how(maxperm=41000))[rows,1:6]
  out
  #      [,1] [,2] [,3] [,4] [,5] [,6]
  # [1,]    8    2    7    1    3    4
  # [2,]    5    7    4    1    6    8
  # [3,]    2    1    8    6    3    4
  # [4,]    5    6    4    1    2    8
  # [5,]    6    1    3    2    5    7
  # [6,]    7    4    6    5    8    2
  

  回复于 2024-04-29T14:13:00+08:00