# [R] Repeating a R code and counting how many repetitions working

varin sacha v@r|n@@ch@ @end|ng |rom y@hoo@|r
Sat Dec 11 16:37:10 CET 2021

```Dear Rui,

I really thank you a lot for your response. Even if I don't fully understand your R code, there is something strange happening while running the code.
Indeed, I have run this R code here below 20 times and I always got the same answer : [1] 1

########################################
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)
}

N <- 100
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,])
########################################

Le samedi 11 décembre 2021, 12:55:43 UTC+1, Rui Barradas <ruipbarradas using sapo.pt> a écrit :

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,

À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 using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help