library(modeest) x<- c(93, 60, 55, 91, 107, 86, 113, 90, 97, 92, 74, 144, 93, 116, 90, 87, 125, 89, 84, 104, 92, 109, 77, 77, 81, 124, 112, 125, 105, 128, 103, 96, 98, 91) boxplot(x, horizontal = TRUE, col = "lightgrey") hist(x, breaks = "freedman-diaconis", col = "lightgrey", border = "black") m = mean(x) mdn = median(x) s = sd(x) skY = (m-mdn)/s cat("Yule's skewness coefficient = ",skY,"\n") library(boot) b<-boot(data = x, function(x, i) { (mean(x[i])-median(x[i]))/sd(x[i])}, R = 1000) b plot(b) boot.ci(b, conf = 0.90, type = c("bca", "norm"), digits = 4) |