marginal_probabilities function

Computes the marginal probabilities of a categorical time series