If you use pm.Normal, the resulting array will contain random variables which will all be univariate, independent draws (but they don’t necessarily need to be identically distributed - you can use arrays in mu and sd to set different means and variances). If you want draws which are correlated across entries, you have to use pm.MvNormal.