Adapting Hawkes process Stan code to PyMC

What do the pair plots look like? If there’s loads of divergences but good r_hats, it can mean there’s some kind of demented geometry in the posterior that the sampler can’t explore (so the samples are biased, because a long tail is missing for example). Sometimes you’ll see it in the pair plots.