probabilité de couverture ?

Bonsoir mes collègues,
J'ai besoin connais comment calculer la probabilité de couverture pour la moyenne empirique (Xbar), moi j'ai utilisé le bootstrap comme ça

Code:

---

1 2 3 4 5 6

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")

---

Merci d'avance.