3

假设我们有 2 个类(A 和 B),每个类有 2 个样本,如下所示:

A A B B
1 2 3 4

我想生成所有可能的唯一排列,其中新类组是旧类的混合,以计算新的统计数据。在这种情况下,我们会得到:

A A B B
1 3 2 4
1 4 2 3

任何其他排列,在 A 类和 B 类之间进行 t 检验(或 ANOVA 等)时,将给出与原始排序相同的结果,我不想要它。有没有一种简单的方法可以做到这一点R

我可以想象生成原始索引向量的所有split排列,对每个排列的样本类执行 s,然后检查样本类的组合是否已经存在,如果不更新排列列表,然后转到下一个排列,等等等等。我只是想检查是否有某种方法可以使用permute包或R.

4

1 回答 1

2

找到了答案。这个partitions包有一些很好的工具来处理这些集合或分区的排列

对于上面的示例,我们在两个分区中有 2 个项目,一个简单的做法是:

library(partitions)
listParts(c(2,2))

这返回

listParts(c(2,2))
[[1]]
[1] (1,4)(2,3)

[[2]]
[1] (1,2)(3,4)

[[3]]
[1] (1,3)(2,4)
于 2012-10-01T14:38:15.440 回答