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)