library(modeest) x<-c(numeric sample data separated by commas) boxplot(x, horizontal = TRUE, col = " lightgrey") hist(x, breaks = "freedman-diaconis", col = " lightgrey", border = "black") m = mean(x) mdn = median(x) s = sd(x) skpII = 3*(m-mdn)/s cat("Pearson's skewness coefficient II = ",skpII,"\n") library(boot) b<-boot(data = x, function(x, i) {3*(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) |