[R] Speed up studentized confidence intervals ?

varin sacha v@r|n@@ch@ @end|ng |rom y@hoo@|r
Sun Dec 19 12:45:38 CET 2021


Dear R-experts,

Here below my R code working but really really slowly ! I need 2 hours with my computer to finally get an answer ! Is there a way to improve my R code to speed it up ? At least to win 1 hour ;=)

Many thanks

########################################################
library(boot)

s<- sample(178:798, 100000, replace=TRUE)
mean(s)

N <- 1000
out <- replicate(N, {
a<- sample(s,size=5)
mean(a)
dat<-data.frame(a)

med<-function(d,i) {
temp<-d[i,]
f<-mean(temp)
g<-var(replicate(50,mean(sample(temp,replace=T))))
return(c(f,g))

}

  boot.out <- boot(data = dat, statistic = med, R = 10000)
  boot.ci(boot.out, type = "stud")$stud[, 4:5]
})
mean(out[1,] < mean(s) & mean(s) < out[2,]) 
########################################################



More information about the R-help mailing list