I ran a regression model and then ran a posterior sample check with new input data using the below code.
posterior_samples = pm.sample_ppc(trace, model = model, samples = 1500)
It looks like the result is a dictionary of posterior samples of the y value.
That said, how do I pull the mean of each observation out?
If I wanted to know the probability that the correct prediction is greater that say 1.2, how would I pull that out?