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)