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)