joint_probabilities function

Computes the joint probabilities of a categorical time series