I am trying to define a model that looks like nvars, ncats = 3, 5 with pm.Model(): α = pm.Exponential('α', lam=λ, shape=(nvars, ncats)) b = pm.Dirichlet('b', a=α) trace = pm.sample() This fails, because the Dirichlet distribution expects an array, not a matrix, in PyMC3 (and in numpy,…

It should work, have you try specifying the shape of b?

Dynamic number of Dirichlet distributions?

junpenglao February 25, 2018, 12:44pm 4

Shape issue is one of our daily struggle