I have generated quantities and the y_pred has around 2000 draws. I’m using the pystan, is there any way that I can plot a specific draw of my posterior in pystan. For example in R we have a useful package that can use a specific draw given the fit : yrep= posterior_predict(fit, draws=50))

Now I’m using this function : fit.plot(‘y_rep’); which returns all the 2000 draws. I only need the mean of all the draws in my plot.