Comment by mvanaltvorst

In my experience, PyMC leads to models that are orders of magnitudes slower than equivalent models written in JAGS. Profiling is also extremely tedious, and there is no section in the PyMC docs that touches upon model performance.

I really like PyMC's API, but as soon as you move towards bigger datasets JAGS or Stan seem to be the only practical options.