Prediction interval for lmer() mixed effects model in R – Cross Validated