kurtosis <- function (x) {
n <- length(x)
K <- (1/n)*sum((x-mean(x))^4)/(((1/n)*sum((x-mean(x))^2))^2)-3
K <- ((n - 1)*((n+1)*K-3*(n-1))/((n-2)*(n-3)))+3
return(K)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.