Which version of pymc are you on? There were some fixes in the latest versions that fixed some sampling issues with Multivariate variables. I don’t know if that applies to you. The code where I used with Multinomial noise sampled fine, did not try it yet with a normal one though.