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?