Hello,

I have defined N random variables that are multivariate Gaussian of dimension 3.

```
import pymc3 as pm
N=10
with pm.Model() as model:
cov = np.array([[1., 0.5,0.5], [0.5, 2,0.5], [0.5,0.5,3]])
mu = np.zeros(3)
theta=pm.MvNormal('Theta', mu=mu,cov=cov, shape=N)
```

What is the correct way to access the i-th RV?

Either:

`theta[i,:]`

Or:

`theta[:,i]`

I am not very good in debugging with PyMC3 so I do not know how to check this.

Thanks for your help and best regards!