我正在寻找一种方法来计算数值变量的四分位距 IQR 的置信区间。当然，它们可以通过引导程序找到，但我明确地在寻找一种不同的、仍然无分发的方式。由于 IQR 是分位数差异，因此该参考可能会朝着正确的方向发展。

伪算法（或 R/Python 代码）会是什么样子？

这是R中的“便宜”引导方式：

library(boot)
set.seed(1)
x <- rnorm(100)
S <- boot(x, function(x, ix) IQR(x[ix]), R = 10000)
boot.ci(S, type = "bca")

# BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS
# Based on 10000 bootstrap replicates
# 
# CALL : 
#   boot.ci(boot.out = S, type = "bca")
# 
# Intervals : 
#   Level       BCa          
# 95%   ( 0.958,  1.448 )  
# Calculations and Intervals on Original Scale

是否有类似@whuber 在他的回答中使用的二项式方法？如何获得百分位数的置信区间？