Unfortunately it might be an issue with Multinomial, see eg this issue: https://github.com/pymc-devs/pymc3/issues/2614