[R] Percentile bootstrap for the median : error message
varin sacha
v@r|n@@ch@ @end|ng |rom y@hoo@|r
Sat Jan 8 18:04:21 CET 2022
Dear R-experts,
Here below my R code for the percentile bootstrap confidence intervals with an error message.
Is there a way to make my R code work ?
Many thanks for your help and time.
############################################
library(boot)
s=rnorm(100000,0,1)
(m<-median(s))
N <- 100
n<-5
out <- replicate(N, {
dat<-data.frame(sample(s,size=n))
med<-function(d,i) {
median(d[i, ])
}
boot.out <- boot(data = dat, statistic = med, R = 10000)
boot.ci(boot.out, type = "per")$per[, 4:5]
})
mean(out[1, ] < m & m < out[2, ])
############################################
More information about the R-help
mailing list