x10<-c(−3.5, −1.219685581, −0.76427577, −0.430514044, −0.139644873, 0.139644873, 0.430514044, 0.76427577, 1.219685581, 3.5) m = mean(x10) mdn = median(x10) s = sd(x10) AY = (m-mdn)/s cat("Coefficient de Yule = ",AY,"\n") length(x10) library(boot) b<-boot(data = x10, function(x10, i) {(mean(x10[i])-median(x10[i]))/sd(x10[i])}, R = 10000) b boot.ci(b, conf = c(0.90, 0.95, 0.99), type = "norm", digits = 4) plot(b) |