Hello,
The problem can be solved with ?replicate.
in the code below I only repeat N <- 1e3, not 1e4.
set.seed(2021)
N <- 1e3
out <- replicate(N, {
boot.out <- boot(data = dat, statistic = med, R = 10000)
boot.ci(boot.out, type = "bca")$bca[, 4:5]
})
mean(out[1,] < mean(s) & mean(s) < out[2,])
Hope this helps,
Rui Barradas
Às 10:47 de 11/12/21, varin sacha via R-help escreveu:
Dear R-experts,
Here below my R code. I am trying to do 2 things :
1) I would like to repeat this R code 10000 times
2) Out of the 10000 repetitions I would have liked R to tell me how many times the "true"
mean value of the population called "s" in my R example here below is included in the
10000 BCa confidence intervals constructed.
Many thanks for your precious help.
########################################
library(boot)
s= sample(178:798, 10000, replace=TRUE)
mean(s)
a=sample(s,size=5)
mean(a)
dat<-data.frame(a)
med<-function(d,i) {
temp<-d[i,]
mean(temp)
}
boot.out<-boot(data=dat,statistic=med,R=10000)
boot.ci(boot.out,type="bca")
########################################
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
