OK will do, thanks.
Using mean_func=tt.zeros_like seems to have gotten the dimensions to be correct but they still aren’t broadcasting correctly. It is still trying to dot product the two matrices rather than elementwise. I do not get an error until I start sampling as opposed to declaring the model as before. edit: Not sure if this is related to the column vector broadcasting issue you mention. Can I reshape the vector in some way?