Posterior fit retrieval?


I have stan_glmer object. I know how to extract the distribution of each parameters .

But I would like to to be able to obtain the posterior (predictive) distribution of the fit for each data point (so, combination of parameters and covariates). If I use the ‘fitted.values’ method I only get point estimates, but I would like to be able to retrieve the full (conditional) posterior for each data point.

Is there a straightforward way to do this?