Regression predicted values in pymc

[cluhmann]

I think you can just do something like this?

az.plot_posterior(ppcs_at_190, var_names=["y", "μ"], ax=ax)

You should be able to add the labels afterwards.