这是 R 中（中值 - 均值）差异的引导置信区间：

z = function() {s = sample(women$weight, replace=TRUE); median(s)-mean(s)}
k = replicate(10000, z())
c(quantile(k, c(.025, .5, .975)), mean=mean(k), sd=sd(k), qgte0=mean(k>=0))

     2.5%       50%     97.5%      mean        sd     qgte0 
-7.933333 -1.333333  5.800000 -1.218007  3.513462  0.362100 

我仍在思考是否k可以在 Wald(-like) 测试中使用差异重采样的平均值和 SD，或者是否可以将大于或等于 0 的分位数视为在某些情况下的单边 p 值假设——欢迎对此发表评论。

可以轻松设置排列测试以使用（均值 - 模式差）作为其测试统计量。这将为您提供差异的精确 P 值。