Adapting Hawkes process Stan code to PyMC

Pass divergences=True to az.pair_plot so you can see if they are all clustered in one place