# NOT RUN {
## Skewness = 0
x <- rnorm(1000)
skewness(x, method = "bickel", M = shorth(x))
## Skewness > 0 (left skewed case)
x <- rbeta(1000, 2, 5)
skewness(x, method = "bickel", M = betaMode(2, 5))
## Skewness < 0 (right skewed case)
x <- rbeta(1000, 7, 2)
skewness(x, method = "bickel", M = hsm(x, bw = 1/3))
# }
Run the code above in your browser using DataLab