calculate_subfeatures function

Computes several subfeatures associated with a categorical time series