Calculate univariate or multivariate (Mardia's test) skew and kurtosis for a vector, matrix, or data.frame