如何在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

if分位数函数计算估计量，如果设置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)