!date import numpy as np, pymc as pm %matplotlib inline # what pymc version is this? pm.__version__ np.random.seed(123456) a = pm.Normal('a', 0, 1) pm.MCMC([a]).sample(1000) np.random.seed(123456) b = pm.Normal('b', 0, 1) pm.MCMC([b]).sample(1000) np.sum((a.trace() - b.trace())**2)