我想使用蒙特卡洛模拟对存在于三个组中的数据集执行置换测试(顺序无关紧要)。在零假设下,数据在三组中的位置无关紧要,因此我可以为每种可能的组合计算检验统计量,并找到统计量的分布函数。
我尝试在for循环中使用combin函数,但是问题是我得到了很多重复项,由于我的数据位于三个向量中,因此很难删除,而第一个向量可能是重复项,而后两个不是,例如,不是获得必需的(3!)^ 7 = 279.936唯一排列,我将计算116280 * 3432(= 399.072.960)测试统计信息。
我没有提供测试统计信息,因为它看起来似乎并不相关,我想自己做一次练习,但是我需要一些帮助,以找到可以将数据矢量划分成的所有组合。
#dataset
A <- c(2.90,2.56,2.88,2.73,2.50,3.18,2.83)
B <- c(2.97,2.45,2.76,2.20,2.16,2.89,2.87)
C <- c(2.67,2.62,1.84,2.33,1.27,2.39,2.39)
z <- c(A,B,C)