function(x, xname = "Skewness Analysis") { x <- sort(x) l <- length(x) ind <- 1:floor(l/2) dni <- l + 1 - ind sk <- (x[ind] + x[dni])/2 - median(x) sp <- (x[dni] - x[ind])/2 sk <- c(0, sk) sp <- c(0, sp) plot(sp, sk, xlab = "SPREAD", ylab = "SKEWNESS", sub = xname) lines(c(0, sp[2]), c(0, 0)) }