The KURTOSIS function is part of the IMLMLIB library . The KURTOSIS function returns the sample kurtosis for each column of a matrix. The sample kurtosis measures the heaviness of the tails of a data distribution. The KURTOSIS function returns an estimate for the excess kurtosis, which is 3 less than the standardized fourth central moment.
The KURTOSIS function returns the same sample kurtosis as the UNIVARIATE procedure. For a formula, see the section "Descriptive Statistics" in the chapter "The UNIVARIATE Procedure" in Base SAS Procedures Guide: Statistical Procedures.
The following example computes the kurtosis for each column of a matrix:
x = {1 0,
2 1,
4 2,
8 3,
16 . };
kurt = kurtosis(x);
print kurt;