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?
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!