如何在 R 中将引导 CI 添加到此函数中

Question

我想计算这两个估计的 90% 置信区间 (Bootstrap)。 Harrell-Davis 无分布分位数估计器。我认为这个函数是估计均值和标准差的非参数分位数的自举版本。现在我想知道如何计算 90% CI？

library(Hmisc)
x <- runif(100)
hdquantile(x, probs =  seq(0.025, 0.975,0.95), se=TRUE,names = TRUE, weights=FALSE)

Answer 1

如果您设置 se=TRUE，hdquantile 函数会计算估计量并给出标准误差。要找到置信区间，您可以通过 boot 函数引导它并通过 boot.ci 获取置信区间

代码

library(Hmisc)
x <- runif(100)
hdquantile(x, probs =  seq(0.025, 0.975,0.95), se=TRUE,names = TRUE, 
weights=FALSE)

为引导创建一个统计数据

library(boot)
hq <- function(x,i) {
hdquantile(x[i], probs =  seq(0.025, 0.975,0.95), se=TRUE,names = TRUE, 
weights=FALSE)
}
bootx <- boot(x,hq,1000)
boot.ci(bootx, conf = 0.90)