如何将自举CI添加到此函数R

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

hdquantile函数计算估计量，如果设置se = TRUE，也会给出标准错误。要找到通过启动功能的置信区间，你可以引导它通过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)