Any suggestions on speeding up the model?

Thank you for the reply. Does PyMC provide a function, like profiler in MALTAB, to visualize the time consumption on each line and help to locate the bottleneck?