In pymc3, a stochastic variable of array shape say 3 can be generated as follows
y = Normal('y', mu, sigma, shape=3, observed=some_data)
Now suppose that y depends on an array of parameters mu = [1,2,3] and sigma = [4,5,6] instead of single values, how would I specify it?