Regression predicted values in pymc

Questions version agnostic
,
6

Does this work?

with model_adelie_flipper_regression:
    pm.set_data({"data": [190]})
    ppcs_at_190 = pm.sample_posterior_predictive(
        idata_r.posterior, var_names=["y", "μ"]
    )

fig, ax = plt.subplots()
az.plot_dist(
    ppcs_at_190.posterior_predictive["y"],
    label="Posterior Predictive of \nIndividual Penguin Mass",
    ax=ax,
)
az.plot_dist(
    ppcs_at_190.posterior_predictive["μ"],
    label="Posterior Predictive of μ",
    color="C4",
    ax=ax,
)