obtain_gamma function

Obtain the transition probabilities gamma for a fitted multilevel HMM