Multivariate normal distribution

[lucianopaz]

Ah, sorry. You have to pass in numpy.array's instead of plain list's. Try this:

mu0 = pm.Normal(‘mu0’, np.array([0.05, 1, 50]), np.array([0.06, 0.1, 0.5]), shape=3)