It sounds like you are looking for the posterior distribution of your control_beta parameters? You can view their estimates by calling az.summary(idata, var_names = “control_beta”), and you can also get their estimates by using idata[“posterior”].mean(dim=["chain", "draw"])
1 Like