我正在从事一个研究项目，我想确定两个分布的等价性。我目前正在使用 Mann-Whitney Test for Equivalence 并且我正在运行的代码（下面）是由 Stefan Wellek (2010) 所著的 Test Statistical Hypotheses of Equivalence and Noninferiority 一书提供的。在运行我的数据之前，我正在使用具有相同均值和标准差的随机正态分布测试此代码。我的问题是有三个嵌套的 for 循环，当运行更大的分布时（如下例所示），代码需要永远运行。如果我只需要运行一次就不会出现这样的问题，但我正在做一个模拟测试并创建功率曲线，所以我需要运行此代码的多次迭代（大约 10,000）。目前，根据我如何改变分布大小，运行 10 个需要几天的时间，
任何有助于提高性能的帮助将不胜感激。
x <- rnorm(n=125, m=3, sd=1)
y <- rnorm(n=500, m=3, sd=1)
alpha <- 0.05
m <- length(x)
n <- length(y)
eps1_ <- 0.2 #0.1382 default
eps2_ <- 0.2 #0.2602 default
eqctr <- 0.5 + (eps2_-eps1_)/2
eqleng <- eps1_ + eps2_
wxy <- 0
pihxxy <- 0
pihxyy <- 0
for (i in 1:m)
for (j in 1:n) …