Pymc Regression Of Many Regressions?

I haven't been using PyMC for long, but I was pleased at how quickly I was able to get a linear regression off the ground (this code should run without modification in IPython): im

Solution 1:

Have you tried using a simple list instead of numpy arrays to store the PyMC objects?