R/kurtosis.ade.R

In epade: Easy Plots

kurtosis.ade <-
function (x, na.rm = FALSE)
{
    if (is.matrix(x))
        apply(x, 2, kurtosis.ade, na.rm = na.rm)
    else if (is.vector(x)) {
        if (na.rm)
            x <- x[!is.na(x)]
        n <- length(x)
        n * sum((x - mean(x))^4)/(sum((x - mean(x))^2)^2)
    }
    else if (is.data.frame(x))
        sapply(x, kurtosis.ade, na.rm = na.rm)
    else kurtosis.ade(as.vector(x), na.rm = na.rm)
}