Good evening my colleagues,
I need to know how to calculate the probability of coverage for the empirical mean (Xbar), I used the bootstrap like this:
library(boot)
n=1000
x=rnorm(n,0,1)
a<-boot(x, function(x,i){mean(x[i])}, R=1000)
bci<-boot.ci(a, type = "norm")