Multivariate categorical with observed data

Thanks! There was an error in the first version of the question that I just edited - I had the “shape” of the prior_probs variable as (1,2), but it should be (2,2) as is shown now