joint_probabilities function

Computes the joint probabilities of an ordinal time series