Any reason why the HPD in Pymc3 differs from what the code in Statistical Rethinking gets?

No problem here

code-4

PyMC3 version 3.3