如何在R中将bootstrap 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

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

代码

import bin ascii
n = int('01100001011000100110001101100100', 2)
binascii.unhexlify('%x' % n)

用于启动创建统计信息

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