I still don’t understand your problem. If I use pm.Normal('RV1', shape=3) or even pm.Normal('RV1', shape=(3,3)) I only get a single trace, not one for each dimension, which makes me think you have a different variable for each subject and are plotting the whole trace. Can you post your code?