ordinal_skewness function

Computes the estimated skewness of an ordinal time series