Regression predicted values in pymc

[cluhmann]

I think you can use arviz’s plot_posterior() function instead of plot_dist() and you get HDIs by default.