如何將自舉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)